c=(4,3,2,1)
w=c(4,3,2,1)
w2 = w / sum(w)
w2
c=(2,1.5,1,0.5)
w=c(2,1.5,1,0.5)
w2 = w / sum(w)
w2
w=c(2,1.75,1.50,1.25)
w2 = w / sum(w)
w2
W = c(.2,.35,.25,.01)
normW = W / sum(W)
normW
library(tmvtnorm)
env
Sys.env
System.env
library(tmvtnorm)
tmp = function(G){print(G)}
tmp(matrix(nrow=2,ncol=2))
tmp = function(G){G1=G;print(G)}
tmp(matrix(nrow=2,ncol=2))
a=c(log(1/6)-log(.18),log(1/6)-log(.14),rep(log(1/6)-log(.17),4))
a
theta=c(.18,.14,rep(.17,4))
theta
sum(a*theta)
atmp = a-sum(a*theta)
atmp
atmp^2
sum(theta*atmp^2)
libary(bas)
library(bas)
help bas
??bas
library(MCMCpack)
library(tmvtnorm)
libraruy(tmvtnorm)
library(tmvtnorm)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,1)
fakedata(6,1,.8,,10,10,1)
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
help(rbinom)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
runif(10,0,1)>0.8
as.numeric(runif(10,0,1)>0.8)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
fakedata(6,1,.8,10,10,1)
data = rbind(fakedata(6,1,.8,10,29,1),fakedata(4,1,.6,10,29,2),fakedata(2,1,.1,10,29,3))
data
data = valdata(c(2,4,6),c(1,1,1),c(.1,.6.9),50,p,c(1,2,3))
data = valdata(c(2,4,6),c(1,1,1),c(.1,.6,.9),50,p,c(1,2,3))
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
data = valdata(c(2,4,6),c(1,1,1),c(.1,.6,.9),50,p,c(1,2,3))
data = valdata(c(2,4,6),c(1,1,1),c(.1,.6,.9),50,29,c(1,2,3))
data
length(c(2,4,6))
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
data = valdata(c(2,4,6),c(1,1,1),c(.1,.6,.9),50,29,c(1,2,3))
data
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
data = valdata(c(2,4,6),c(1,1,1),c(.1,.6,.9),50,29,c(1,2,3))
data
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
data
data = valdata(c(2,4,6),c(1,1,1),c(.1,.6,.9),50,29,c(1,2,3))
data
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
data = valdata(c(2,4,6),c(1,1,1),c(.1,.6,.9),50,29,c(1,2,3))
data
outdata = valdata(c(2,4,6),c(1,1,1),c(.1,.6,.9),50,29,c(1,2,3))
outdata
means=c(2,2,2)
vars=c(1,1,1)
probs=c(.1,.6,.9)
n=50
p=29
names=c(1,2,3)
for(i in 1:length(means)){rbind(outdata,fakedata(means[i],vars[i],probs[i],n,p,names[i]))}
outdata
i=1
fakedata(means[i],vars[i],probs[i],n,p,names[i])
rbind(NULL,fakedata(means[i],vars[i],probs[i],n,p,names[i]))
i=2
rbind(NULL,fakedata(means[i],vars[i],probs[i],n,p,names[i]))
i=3
rbind(NULL,fakedata(means[i],vars[i],probs[i],n,p,names[i]))
source("/Users/andrewcron/Desktop/AnalysisSoftware/Validation/fakedata.R")
outdata = valdata(c(2,4,6),c(1,1,1),c(.1,.6,.9),50,29,c(1,2,3))
outdat
outdata
env()
Sys.env()
??env
Sys.getenv()
log(.1)-log(1-.9)
.1 - (1-.9)
.1-(1-.9)
log(.1)-log(1-.9)
help(get)
X=matrix(rep(1,9),3,3)
X
Y=matrix(rep(1,9),3,3)
X%*%Y
for(i in 1:1000){
for(j in 1:1000){
X%*%Y}}
sum(1:1000000)
sum(as.numeric(1:1000000))
seq(1,10,0.001)
library(tmvtnorm)
help(tmvtnorm)
??tmvtnorm
dtmvnorm
help(dmvnrom)
help(dmvnorm)
x=seq(1,7,0.01)
help(dmvtnorm)
help(dvtnorm)
help(mvtnorm)
help(tmvtnorm)
help(dtmvnorm)
plot(x,dtmvnorm(x,4,2,1,7))
plot(x,dtmvnorm(t(x),4,2,1,7))
??tmvtnorm
help(dtmvnrom.marginal)
help(dtmvnorm.marginal)
help.search("truncated normal")
help.search("truncated")
mtmvnorm(2,5,1,7)
plot(x,dnorm(x,2,5),type="l")
lines(c(2,2),c(0,.08))
lines(c(2,2),c(0,.08),col="r")
lines(c(2,2),c(0,.08),color="r")
lines(c(2,2),c(0,.08),col="red")
lines(c(3.11,3.11),c(0,0.08))
lines(c(3.11,3.11),c(0,0.08),col="red")
plot(x,dnorm(x,2,5),type="l",xlab="Survey Question",ylab="Posterior Density")
lines(c(1,1),c(0,dnorm(1,2,5)))
lines(c(7,7),c(0,dnorm(7,2,5)))
lines(c(2,2),c(0,0.085),col="red")
lines(c(3.11,3.11),c(0,0.085),col="blue")
legend("topright",lty=c(1,1),col=c("red","blue"),legend=c("Non-Truncated Mean",Truncated Mean))
help(legend)
legend("topright",lty=c(1,1),col=c("red","blue"),legend=c("Non-Truncated Mean","Truncated Mean"))
??format
format(pi,digits=3)
printf
fprint
sprintf
time
help(time)
date()
x=dat()
x=date()
x
y=date()
y-x
as.numeric(x)
time(x)
x
help(date)
Sys.time()
x=time()
help(date)
help(proc.time)
system.time
x=system.time
x
x=system.time()
help(system.time)
x=proc.time()
x-proc.time()
system.time(x=1)
proc.time()-x
proc.time()
(1/7)^3
sqrt(.5/(100*99))
s3*qrt(.5/(100*99))
3*sqrt(.5/(100*99))
3*sqrt(.5/(100))
3*sqrt(.5/(200))
3*sqrt(.5/(300))
3*sqrt(.5/(400))
3*sqrt(.5/(500))
3*sqrt(.5/(750))
3*sqrt(.5/(800))
3*sqrt(.5/(2000))
3*sqrt(.5/(20000))
3*sqrt(.5/(1000))
3*sqrt(0.005)
3*sqrt(.5/(50))
3*sqrt(.5/(20))
.5*.5/40
.5*.5/50
library(sn)
help(sn)
help(SN)
x=seq(-3,3,0.01)
plot(x,dsn(x),type="l")
plot(x,dsn(x,0,1,-1),type="l")
plot(x,dsn(x),type="l")
for i in seq(-1,1,0.01)
for (i in seq(-1,1,0.01))
plot(x,dsn(x,0,1,i))
for (i in seq(-1,1,0.01))
lines(x,dsn(x,0,1,i))
plot(x,dsn(x),type="l")
for (i in seq(-1,1,0.01))
lines(x,dsn(x,0,1,i),type="l")
library(sn)
x=seq(-3,3,0.001)
plot(x,dsn(x,0,1,-1))
plot(x,dsn(x,0,1,-1),type="l'")
plot(x,dsn(x,0,1,-1),type="l")
plot(x,dsn(x,0,1,10),type="l")
line(c(0,0),c(0,1),type="l")
lines(c(0,0),c(0,1),type="l")
lines(x,dsn(x,0,1,5))
lines(x,dsn(x,0,1,1))
lines(x,dsn(x,0,1,100))
lines(x,dsn(x,0,1,50))
help(dsn)
setwd("~/Desktop/CDP_BETA/Wrappers/RWrapper/")
source("/Users/andrewcron/Desktop/CDP_BETA/Wrappers/RWrapper/GPUmix.R")
