search for: hmtest

...uot;, "T(70%)c-N(0%)f", "T(80%)a-N(0%)f", "T(90%)-N(0%)f", "T(70%)a-N(0%)i", "T(70%)c-N(0%)i", "T(90%)-N(0%)i", "T(70%)a-N(0%)c", "T(70%)c-N(0%)c", "T(80%)a-N(0%)c")) rownames(sig.data)<-sig.data[,2] my.hmtest <- structure(list( estimate = t(t(structure(sig.data[,"estimate"], .Names = rownames(sig.data)))), conf.int = sig.data[,4:5], ctype = "ABCC4-2007"), class = "hmtest") par(mex=0.5) #This helps to accomodate the margins when text is getting cut off plot(my.hmt...

...= "within", data = new.frame,index = c("id")) re=plm(gd ~ l+g+o+c+g1+h+n+r, model = "random", data = new.frame,index = c("id"),random.method="amemiya") then I wrote the following function for the Hausman Test using an auxiliary regression method: hmtest=function(re=0,fe=0){ y.re=pmodel.response(re) X.re=model.matrix(re) X.fe=model.matrix(fe) auxdata<-data.frame(cbind(y.re,X.re,X.fe)) colnames(auxdata)<-c("y", paste("x", 1:17, sep="")) auxmod<-lm(y~x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11...

Here it is again, hope this is more clear I am using the following data (only a small subset is given): Habitat Fungus.yield Birch 20.83829053 Birch 22.9718181 Birch 22.28216829 Birch 24.23136797 Birch 22.32147961 Birch 20.30783598 Oak 27.24047258 Oak 29.7730014 Oak 30.12608508 Oak 25.76088669 Oak 30.14750974 Hornbeam 17.05307949 Hornbeam 15.32805111 Hornbeam 18.26920177 Hornbeam 21.30987049