Displaying 3 results from an estimated 3 matches for "pval1".
Did you mean:
val1
2012 Mar 04
1
p-value from GLM
Dear all,
I am fitting a GLM similar to
library(MASS)
anorex.1 <- glm(Treat~Postwt+Prewt,family = binomial, data = anorexia)
I have found two ways of computing the p-value of the fitted model:
pval1 <- 1-pchisq(anorex.1$deviance,anorex.1$df.residual)
pval2 <- 1-pchisq(anorex.1$null.deviance - anorex.1$deviance,
anorex.1$df.null - anorex.1$df.residual)
pval2 is testing LR chi2 from the null model, but what does pval1 says?
Many thanks in advance,
Dunia
[[alternative...
2011 Dec 03
1
partial mantel tests in ecodist with intential NA values.
...00,20)
> y<-sample(1:1000,20)
> z<-sample(1:1000,20)
> M1<-as.matrix( distance(x))
> M2 <-as.matrix( distance(y ))
> M3<-as.matrix( distance(z ))
> D1<-(lower(M1))
> D2<-(lower(M2))
> D3<-(lower(M3))
> mantel(D1 ~ D2+D3, nperm=1000)
mantelr pval1 pval2 pval3 llim.2.5% ulim.97.5%
0.09014696 0.10300000 0.89800000 0.18400000 0.01857311 0.18468621
> partial.mantel.test(M1,M2,M3,quiet=T)
$MantelR
r12 r13 r23 r12.3 r13.2
0.08977575 0.02170997 -0.01561346 0.09014696 0.02320821
$p
[1] 0.095...
2011 Jan 20
1
Problems with ecodist
Dear Dr.Goslee and anyone may intrested in matrix manipulate,
I am using your ecodist to do mantel and partial mantel test, I have
locality data and shape variation data, and the two distance matrixs are
given as belowings. When I run the analysis, it is always report that the
matrix is not square, but I didn't know what's wrong with my data. Would you
please help me on this. I am quite