[R] R-square prob is not calculated by randomization in lmPerm::lmp

Jeff Newmiller jdnewmil at dcn.davis.ca.us
Tue Sep 6 16:35:32 CEST 2016


That is contributed code. It could do anything the author felt like. I recommend reading the source code. 
-- 
Sent from my phone. Please excuse my brevity.

On September 5, 2016 11:52:15 PM PDT, Agustin Lobo <aloboaleu at gmail.com> wrote:
>Any reason why the R-square prob is not calculated by randomization in
>lmPerm::lmp? The help pages states "Either permutation test p-values
>or the usual F-test p-values will be output", but I always get the F
>test for R-square as with lm():
>
>require(lmPerm)
>x <- 1:1000
>set.seed(1000)
>y1 <- x*2+runif(1000,-100,100)
>dat <- data.frame(x =x,y=y1)
>summary(lmp(y~x, data=dat,center=FALSE,perm="Prob"))
>
>[1] "Settings:  unique SS "
>
>Call:
>lmp(formula = y ~ x, data = dat, center = FALSE)
>
>Residuals:
>     Min       1Q   Median       3Q      Max
>-100.431  -48.645    2.843   48.640  101.800
>
>Coefficients:
>  Estimate Iter Pr(Prob)
>x    1.993 5000   <2e-16 ***
>---
>Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
>Residual standard error: 57.3 on 998 degrees of freedom
>Multiple R-Squared: 0.9902, Adjusted R-squared: 0.9902
>F-statistic: 1.009e+05 on 1 and 998 DF,  p-value: < 2.2e-16



More information about the R-help mailing list