R: Step function not producing complete result output
0
0
Entering edit mode
6.0 years ago
wkmustahs21 ▴ 30

Hi

I am having issues with the step() in R I am running this lmer model:

lmer(TSW ~  Entry +(1|Block), data=inputData)
step(lmer(TSW ~  Entry +(1|Block), data=inputData))

I should be getting some results on Least Sq Means Differences but it's not showing up at all. I don't understand why. This is the result I get:

Backward reduced random-effect table:

            Eliminated npar   logLik    AIC    LRT Df Pr(>Chisq)   
<none>                   55  -8.8476 127.69                        
(1 | Block)          0   54 -12.1910 132.38 6.6868  1   0.009713 **
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Backward reduced fixed-effect table:
Degrees of freedom method: Satterthwaite 

      Eliminated Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
Entry          0 27.243 0.52391    52 10.181  7.5053 0.0007433 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Model found:
TSW ~ Entry + (1 | Block)
> st <- step(inputData.entry,  reduce.fixed=FALSE);
> print(st);
Backward reduced random-effect table:

            Eliminated npar   logLik    AIC    LRT Df Pr(>Chisq)   
<none>                   55  -8.8476 127.69                        
(1 | Block)          0   54 -12.1910 132.38 6.6868  1   0.009713 **
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Backward reduced fixed-effect table:
Degrees of freedom method: Satterthwaite 

      Eliminated Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
Entry          0 27.243 0.52391    52 10.181  7.5053 0.0007433 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Model found:
TSW ~ Entry + (1 | Block)

Can anyone tell me what is going on? I have tried to restart R and everything to run this script like 4 times but I am not getting the expected results.

**For additional reference please go here: https://passel.unl.edu/pages/informationmodule.php?idinformationmodule=1130447276&topicorder=6&maxto=7&minto=1

There is a video please jump to 6:09 and you will see the results that should come up.

Thanks again.

R least means • 1.2k views
ADD COMMENT

Login before adding your answer.

Traffic: 2505 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6