r.tmp<-me.test.lin1(1000)
r.tmp[2]/r.tmp[1]
r.tmp<-me.test.lin1(10)
r.tmp[2]/r.tmp[1]
r.tmp<-me.test.lin1(10000)
r.tmp[2]/r.tmp[1]
r.tmp[2]/r.tmp[1]
r.tmp<-me.test.lin1(10000)
r.tmp[2]/r.tmp[1]
r.tmp<-me.test.lin1(10000)
r.tmp[2]/r.tmp[1]
sqrt9.9)
sqrt(.9)
1/.81
fix(me.test.lin1)
me.test.lin1(10000)
me.test.lin1(1000)
me.test.lin1(100)
fix(me.test.lin1)
me.test.lin1(100)
fix(me.test.lin1)
me.test.lin1(100)
round(me.test.lin1(100),3)
round(me.test.lin1(1000),3)
round(me.test.lin1(100),5)
round(me.test.lin1(1000),5)
fix(me.test.lin1)
round(me.test.lin1(1000),5)
round(me.test.lin1(100),5)
q()
install.packages("MplusAutomation")
library(MplusAutomation)
runModels_Interactive()
extractModelParameters("d:/jouni/projects/lcat/computing/mplus/tmp/ex3.1.out")
help(runModels)
runModels_Interactive()
runModels("d:/jouni/projects/lcat/computing/mplus/tmp/")
extractModelParameters("d:/jouni/projects/lcat/computing/mplus/tmp/lcatm.3.1a.out")
extractModelParameters("d:/jouni/projects/lcat/computing/mplus/tmp/ex3.1.out")
runModels_Interactive()
runModels_Interactive()
runModels("d:/jouni/projects/lcat/computing/mplus/tmp/")
extractModelParameters("d:/jouni/projects/lcat/computing/mplus/tmp/lcatm.2.8m.out")
help(getSavedata_Data)
help(MplusAutomation)
help(read.csv)
fix(f.tmp0
fix(f.tmp)
f.tmp<-edit()
f.tmp()
ls()
fix(ess4.3c.dat)
prepareMplusData(ess.3c.dat,"d:/jouni/projects/lcat/r_mplus/tests/ess4_3c.dat")
prepareMplusData(ess4.3c.dat,"d:/jouni/projects/lcat/r_mplus/tests/ess4_3c.dat")
prepareMplusData(ess4.3c.dat,"d:/jouni/projects/lcat/computing/r_mplus/tests/ess4_3c.dat")
runModels("d:/jouni/projects/lcat/computing/r_mplus/tests",replaceOutfile="always")
readModels("d:/jouni/projects/lcat/computing/r_mplus/tests")
runModels("d:/jouni/projects/lcat/computing/r_mplus/tests",replaceOutfile="always")
runModels("d:/jouni/projects/lcat/computing/r_mplus/tests",replaceOutfile="always",showOutput=T)
runModels("d:/jouni/projects/lcat/computing/r_mplus/tests",replaceOutfile="always",showOutput=T)
r.tmp<-readModels("d:/jouni/projects/lcat/computing/r_mplus/tests")
r.tmp
extractModelSummaries
extractModelSummaries_1file
extractSummaries_1file
getSavedata_Fileinfo("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a.dat")
getSavedata_Fileinfo("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a")
getSavedata_Fileinfo("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a.out")
r.tmp
getSavedata_Fileinfo("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a.out")
extractModelParameters("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a")
extractModelParameters("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a.out")
extractModelParameters("d:/jouni/projects/lcat/computing/r_mplus/tests")
extractModelParameters("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a.out")
extractModelSummary("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a.out")
extractModelSummaries("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a.out")
getSavedata_Data("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a.out")
getSavedata_Fileinfo("d:/jouni/projects/lcat/computing/r_mplus/tests/lcatm.3.1a.out")
showSummaryTable(r.tmp)
extractModelSummaries("d:/jouni/projects/lcat/computing/r_mplus/tests")
r.tmp<-extractModelSummaries("d:/jouni/projects/lcat/computing/r_mplus/tests")
showSummaryTable(r.tmp)
r.tmp
showSummaryTable(r.tmp,keepCols=c("Parameters","LL","AIC","BIC"))
r.tmp
class(r.tmp)
r.tmp<-extractModelSummaries("d:/jouni/projects/lcat/computing/r_mplus/tests")
r.tmp
showSummaryTable(r.tmp,keepCols=c("Parameters","LL"))
showSummaryTable(r.tmp,keepCols=c("Parameters","LL"),sortBy=NULL)
showSummaryTable(r.tmp,keepCols=c("Parameters","LL"),sortBy=NA)
showSummaryTable(r.tmp,keepCols=c("Parameters","LL"),sortBy="Parameters")
HTMLSummaryTable(r.tmp,keepCols=c("Parameters","LL"),sortBy="Parameters")
getwd()
LatexSummaryTable(r.tmp,keepCols=c("Parameters","LL"),sortBy="Parameters")
q()
ls()
i
help(dump)
dump("ess4.3c.dat","dumpdata.R")
q()
duplication.matrix
library(matrixcalc)
install.package(matrixcalc)
install.packages(matrixcalc)
install.packages("matrixcalc")
library(matrixcalc)
duplication.matrix
q()
rmultinomial
getwd()
rmultinomial2
.51/2.25
 -.5102237/.2263979
2*(1-qnorm( -.5102237/.2263979))
2*(1-pnorm( -.5102237/.2263979))
qnorm( -.5102237/.2263979)
help(qnorm)
pnorm( -.5102237/.2263979)
2*pnorm( -.5102237/.2263979)
exp(-.5102237)
exp(-.5102237)*.2263979
-.5102237+c(-1,1)*qnorm(.975)*.2263979
exp(-.0165245)
exp(c(-.9539554,-.0664920))
-.0053501/2.05
-.0053501/.0026098
2*qnorm(-.0053501/.0026098)
2*pnorm(-.0053501/.0026098)
-.0053501+c(-1,1)*qnorm(.975)*.0026098
.9946642*.0026098
exp(c(-.0104652,-.000235))
q()
x.tmp<-seq(1,5,len=100)
f.tmp<-3*x.tmp^2-x.tmp^3/3-5*x.tmp
plot(x.tmp, f.tmp, type='l')
plot(x.tmp[-1], x.tmp[-1]-x.tmp[-100], type='l')
summ(x.tmp[-1]-x.tmp[-100])
summary(x.tmp[-1]-x.tmp[-100])
f.tmp<-6*x.tmp-x.tmp^2-5
plot(x.tmp, f.tmp, type='l')
help(chisq.test)
x.tmp<-rpoisson(100,2)
x.tmp<-rpois(100,2)
x.tmp
chisq.test(x.tmp,simulate.p.value=T)
chisq.test(x.tmp,simulate.p.value=T,B=10000)
x.tmp<-rpois(100,5)
x.tmp
chisq.test(x.tmp,simulate.p.value=T,B=10000)
dbinom(0:6,n=1,p=1/93)
help(dbinom)
dbinom(0:6,1,p=1/93)
100*dbinom(0:6,1,p=1/93)
q()
ls()
q()
1/(1+exp(-1.677)+exp(-1.893))
help(sample)
sample(1:26,26)
q()
ls()
fix(f.tmp)
source("d:/jouni/work/r/firstlast.r")
.First()
q()
fix(f.tmp)
fix(gssmac.p)
gssmac.p()
fix(gssmac.p)
gssmac.p()
fix(gssmac.p)
gssmac.p()
fix(gssmac.p)
gssmac.p(.25)
fix(gssmac.p)
gssmac.p(.25)
fix(gssmac.p)
gssmac.p(.25)
fix(gssmac.p)
gssmac.p(.25)
fix(gssmac.p)
q()
ls()
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
gssmac.p(.2)
fix(gssmac.p)
help(min)
fix(gssmac.p)
gssmac.p(.2)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
colours()
fix(gssmac.p)
colours()
fix(gssmac.p)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
fix(gssmac.p)
gssmac.p(.2,T)
ls()
fix(mi452.ex)
mi452.ex(30,30)
mi452.ex(25,25)
mi452.ex(30,25)
mi452.ex(30,100)
mi452.ex(0,0)
mi452.ex(20,20)
mi452.ex(40,2)
mi452.ex(40,30)
mi452.ex(20,40)
mi452.ex(30,1)
mi452.ex(30,10)
mi452.ex(20,1)
mi452.ex(25,2)
mi452.ex(30,3)
mi452.ex(20,1)
mi452.ex(1,1)
mi452.ex(20,2)
mi452.ex(25,3)
mi452.ex(30,1000)
mi452.ex(30,1)
mi452.ex(20,1)
mi452.ex(10,10)
mi452.ex(1,1)
mi452.ex(0,0)
mi452.ex(40,3)
mi452.ex(20,10)
mi452.ex(30,30)
mi452.ex(20,2)
mi452.ex(0,0)
mi452.ex(30,2)
mi452.ex(40,30)
mi452.ex(30,5)
mi452.ex(20,10)
mi452.ex(30,1)
mi452.ex(20,1)
mi452.ex(10,100)
ls()
update.packages(checkBuilt=TRUE, ask=FALSE)
q()
local({pkg <- select.list(sort(.packages(all.available = TRUE)),graphics=TRUE)
if(nchar(pkg)) library(pkg, character.only=TRUE)})
q()
ls()
ls()
getwd()
q()
update.packages(checkBuilt=TRUE, ask=FALSE)
help(update.packages)
help(update.packages)
help(help)
options()
options(browser = as.vector(Sys.getenv("R_BROWSER")))
options(browser)
options()$browser
Sys.getenv("R_BROWSER")
help(help)
q()
options()$browser
help(options)
getoptions(browser)
getOptions(browser)
getOption(browser)
getOption("browser")
help(help.start)
help(install.packages)
update.packages(checkBuilt=TRUE, ask=FALSE)
q()
help(help)
help(help)
getwd()
c(289,44,40,10,221,74,63,8)
c(289,44,40,10,221,74,63,8)/749
round(c(289,44,40,10,221,74,63,8)/749,2)
round(c(289,44,40,10,221,74,63,8)/749,3)
sum(round(c(289,44,40,10,221,74,63,8)/749,3))
sum(round(c(289,44,40,10,221,74,63,8)/383,3))
round(c(289,44,40,10,221,74,63,8)/383,3)
sum(round(c(289,44,40,10,221,74,63,8)/383,3)[1:4])
round(c(289,44,40,10,221,74,63,8)/366,3)
sum(round(c(289,44,40,10,221,74,63,8)/366,3)[5:8])
0.202*0.755 + 0.172*(0.755 + 0.115)+0.022*(0.755 + 0.115+ 0.104)
ppoisson(0,7)
help(poisson)
help(dpois)
dpois(0,7)
dpois(1,7)
dpois(2,7)
dpois(3,7)
exp(-7)
exp(-7)*7
exp(-7)*49/2
exp(-7)*(7^3)/6
sum(dpois(0:3,7))
pnorm(.4,.7,sqrt(.07))
pnorm(.35,.7,sqrt(.07))
pnorm(-1.13)
pnorm(-1.1338)
pnorm(-1.1338934)
qnorm(.95)
qnorm(.837,.7,sqrt(.007))
pnorm(.837,.7,sqrt(.007))
p.logit
p.logit(-0.12+.46*c(2,4))
2.28-0.06*35-.035*20
p.logit(-0.52+.46*c(2,4))
p.logit(-0.52-1.1+(.46-.2)*c(2,4))
help(library)
p1.tmp<-p.logit(1.56-.0039*46-.467)
p1.tmp
p1.tmp<-p.logit(1.56-(-.0039*46-.467))
p1.tmp
p2.tmp<-p.logit(.047-(-.0039*46-.467))
p2.tmp
p1.tmp-p2.tmp
1/(1+exp(-.508+.75+.0077*46)+exp(.997+.878-.0012*46)+exp(.855+.575-.0058*46))
dpois(0:9,2)
100*dpois(0:9,2)
round(100*dpois(0:9,2),2)
sum(round(100*dpois(0:9,2),2))
sum(round(100*dpois(0:9,2),2)[-(1:5))
sum(round(100*dpois(0:9,2),2)[-(1:5)]
)
e.tmp<-round(100*dpois(0:9,2),2)[(1:5)]
e.tmp
e.tmp<-c(e.tmp,5.26)
sum(e.tmp)
o.tmp<-c(0,5,16,30,17,32)
sum((o.tmp-e.tmp)^2/e.tmp)
p.logit
p.logit(-1.879+.0445*33.91)
p.logit(-1.879003+.0445*33.91)
p.logit(-1.879003+.0444766*33.91)
p.logit(-1.879003+.0444766*33.91317)
library(foreign)
help(foreign)
help(read.dta)
read.csv("e:\data.csv")
read.csv("e:\\data.csv")
read.csv("f:\\data.csv")
data.aleks<-read.csv("f:\\data.csv")
names(data.aleks)
data.aleks<-data.aleks[,-(1:3)]
names(data.aleks)
data.aleks[,1]==""
i.tmp<-data.aleks[,1]==""
i.tmp<-data.aleks[,1]!=""
data.aleks<-data.aleks[i.tmp,]
fix(data.aleks)
table(data.aleks[,1])
table(as.numeric(data.aleks[,1]))
fix(f.tmp)
f.tmp()
table(ind)
decision
decision %in% ind
Q
fix(f.tmp)
d.tmp<-f.tmp()
fix(d.tmp)
fix(f.tmp)
d.tmp<-f.tmp()
fix(f.tmp)
fix(d.tmp)
fix(f.tmp)
d.tmp<-f.tmp()
fix(d.tmp)
table(d.tmp$decision)
fix(f.tmp)
help(sort)
fix(f.tmp)
names(d.tmp)
fix(f.tmp)
aleks.f1<-f.tmp
d.tmp<-aleks.f1()
fix(d.tmp)
help(sample)
fix(aleks.f1)
sample(10,2)
fix(aleks.f1)
d.tmp<-aleks.f1()
fix(data)
sampled
n
d.tmp<-aleks.f1()
ls()
print(n)
Q
fix(d.tmp)
fix(aleks.f1)
d.tmp<-aleks.f1()
fix(data)
sampled
cbind(data,sampled)[1:30,]
Q
fix(aleks.f1)
aleks.data.final<-aleks.f1()
fix(aleks.data.final)
fix(aleks.f1)
aleks.data.final<-aleks.f1()
fix(aleks.f1)
fix(aleks.data.final)
fix(aleks.f1)
aleks.data.final<-aleks.f1()
fix(aleks.data.final)
aleks.data.final[aleks.data.final$sampled==1,]
nrow(aleks.data.final[aleks.data.final$sampled==1,])
help(write.csv)
write(aleks.data.final,"f:\sampling.csv")
write(aleks.data.final,"f:\\sampling.csv")
class(aleks.data.final)
write("aleks.data.final","f:\\sampling.csv")
write(aleks.data.final,file="f:\\sampling.csv")
class(aleks.data.final)
aleks.data.final[1:10,]
names(aleks.data.final)
ncol(aleks.data.final)
write(aleks.data.final,file="f:\\sampling.csv",row.names=F)
write.csv(aleks.data.final,file="f:\\sampling.csv")
aleks.data.final<-aleks.f1()
nrow(aleks.data.final[aleks.data.final$sampled==1,])
aleks.data.final[aleks.data.final$sampled==1,]
write.csv(aleks.data.final,file="f:\\sampling.csv")
q()
help(pnorm)
1-pnorm(34.1,20.5,7.6)
qnorm(.166,495,94)
.39+1.96*c(-1,1)*.39*(1-.39)/sqrt(1000)
.39+1.96*c(-1,1)*sqrt(.39*(1-.39)/1000)
(45/700-.084)/(sqrt(.084*(1-.084)/700))
1-pnorm(120,100,15)
.75+1.96*c(-1,1)*sqrt(.75*(1-.75)/1000)
(190-200)/(25/sqrt(40))
install.packages("rmeta")
install.packages("calibrate")
q()
library(jgjh)
library(matcalc)
library()
help(library)
q()
ls()
setwd("z:/research_projects/lcat/lcatworkshop/r")
q()
getwd()
ls()
help(rm)
ls()
remove(list=ls())
ls()
fix(.First)
q()
lcat
print.lcat.list)
print.lcat.list
fix(.First)
q()
q()
