mgaughan/24_deb_pkg_gov
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

model without age, and a few more sim samples

Browse Source
This commit is contained in:
Kaylea Champion 2023-11-14 16:14:31 -08:00
parent d62f7a8b63
commit 68653a6470
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
R/calculatePower.R
View File

@ -61,11 +61,13 @@ g
data2$yearsOld <- data2$age / 365 data2$yearsOld <- data2$age / 365
kmodel2 <- lm(up.fac.mean ~ mmt + milestones + age, data=data1) kmodel2 <- lm(up.fac.mean ~ mmt + milestones + age, data=data1)
kmodel5 <- lm(up.fac.mean ~ mmt + milestones, data=data1)
kmodel4 <- lm(up.fac.mean ~ mmt + age, data=data1) kmodel4 <- lm(up.fac.mean ~ mmt + age, data=data1)
kmodel3 <- lm(up.fac.mean ~ formal.score, data=data1) kmodel3 <- lm(up.fac.mean ~ formal.score, data=data1)
summary(kmodel2) summary(kmodel2)
summary(kmodel3) summary(kmodel3)
summary(kmodel4) summary(kmodel4)
summary(kmodel5)
#pilotM <- glm(up.fac.mean ~ ((mmt) / (milestones/age)), # give the anticipated regression a try #pilotM <- glm(up.fac.mean ~ ((mmt) / (milestones/age)), # give the anticipated regression a try
# family=gaussian(link='identity'), data=data1) # family=gaussian(link='identity'), data=data1)
@ -100,5 +102,9 @@ powerCheck(n, nSims)
#Sample values #Sample values
powerCheck(50, 1000) powerCheck(50, 1000)
powerCheck(80, 1000) powerCheck(200, 1000)
powerCheck(200, 5000) powerCheck(500, 1000)
powerCheck2(50, 1000)
powerCheck2(200, 1000)
powerCheck2(500, 1000)