diff --git a/R/calculatePower.R b/R/calculatePower.R
index 30d4c6e..90be49f 100644
--- a/R/calculatePower.R
+++ b/R/calculatePower.R
@@ -61,11 +61,13 @@ g
 data2$yearsOld <- data2$age / 365
 
 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)
 kmodel3 <- lm(up.fac.mean ~ formal.score, data=data1)
 summary(kmodel2)
 summary(kmodel3)
 summary(kmodel4)
+summary(kmodel5)
 
 #pilotM <- glm(up.fac.mean ~ ((mmt) / (milestones/age)),   # give the anticipated regression a try
 #                  family=gaussian(link='identity'), data=data1)
@@ -100,5 +102,9 @@ powerCheck(n, nSims)
 
 #Sample values
 powerCheck(50, 1000) 
-powerCheck(80, 1000) 
-powerCheck(200, 5000)
+powerCheck(200, 1000) 
+powerCheck(500, 1000)
+
+powerCheck2(50, 1000) 
+powerCheck2(200, 1000) 
+powerCheck2(500, 1000)