Kay, This should not work if the function is correctly written. You say that the terms in your formula are in data=env, but there are no variables env$NFac and env$Fac in env -- but there are NFac and Fac.
Cheers, Jari Sent from my iPad On 23.2.2013, at 12.15, "Kay Cichini" <kay.cich...@gmail.com> wrote: > Hi all, > > Does anyone have a clue why this is not working: > > library(BiodiversityR) > library(vegan) > > data(mite) > mite <- mite[1:60,] > > env <- data.frame(Fac = rep(LETTERS[1:2], each = 30), NFac = > rep(letters[1:4], each = 15)) > > nested.npmanova(mite ~ env$Fac + env$NFac, data = env, method = "jac", > permutations = 1000) > > Regards, > Kay > > -- > > Kay Cichini, MSc Biol > > Grubenweg 22, 6071 Aldrans > > Tel.: 0650 9359101 > > E-Mail: kay.cich...@gmail.com > > Web: > www.theBioBucket.blogspot.co.at<http://www.thebiobucket.blogspot.co.at/><http://www.theBioBucket.blogspot.co.at> > -- > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-ecology mailing list > R-sig-ecology@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-ecology _______________________________________________ R-sig-ecology mailing list R-sig-ecology@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology