Hello Jari, Thank you so much for your help. Yes, I have used percent data. I just know tried with the absolute counts and got, what I guess are right, values of AIC and BIC.
Now coming back to the similar values, how can I discriminate in between the models with so similar AIC and BIC (and Deviance) values. I apologise if you wrote it in the paragraph below, but I did not understand. I copy here the new results: Cluster I RAD models, family poisson No. of species 24, total abundance 166 par1 par2 par3 Deviance AIC BIC Null 91.6634 161.7168 161.7168 Preemption 0.23969 53.4299 125.4832 126.6613 Lognormal 0.71216 1.6811 13.3666 87.4199 89.7761 Zipf 0.42497 -1.4264 2.4005 76.4538 78.8100 Mandelbrot 0.42497 -1.4264 2.9057e-06 2.4005 78.4538 81.9880 Cluster II RAD models, family poisson No. of species 35, total abundance 228 par1 par2 par3 Deviance AIC BIC Null 58.597 165.650 165.650 Preemption 0.13664 45.830 154.884 156.439 Lognormal 1.0417 1.3473 10.898 121.951 125.062 Zipf 0.27724 -1.0959 11.181 122.234 125.344 Mandelbrot 0.3957 -1.221 0.4 10.721 123.774 128.440 Thanks again! Sol On Dec 17, 2013, at 1:48 PM, Jari Oksanen <jari.oksa...@oulu.fi> wrote: > Hello, > On 17/12/2013, at 17:01 PM, Sol Noetinger wrote: > >> Hello, >> >> Cluster II >> RAD models, family poisson >> No. of species 35, total abundance 100 >> >> par1 par2 par3 Deviance AIC BIC >> Null 25.7004 Inf Inf >> Preemption 0.1 27.8760 Inf Inf >> Lognormal 0.21756 1.3473 4.7797 Inf Inf >> Zipf 0.27724 -1.0959 4.9038 Inf Inf >> Mandelbrot 0.64175 -1.3825 1 4.9181 Inf Inf >> >> I read from the manual that to see which models fits better you use the AIC >> values. >> What is the meaning of getting "infinite"? > > It seems you can get infinite AIC and BIC when the distribution you selected > is in conflict with the nature of your data. For instance, when you postulate > a Poisson model (like in your case), but your data are not integers (counts). > Was that the case with you? Distribution families that go along with > non-integer (real) data are gaussian and Gamma. You can neither use > quasipoisson nor other quasi models, because these do not have AIC. > >> Can I use the Deviance value to compare the models? >> And in case I can use the deviance, since there are very close values, >> should I run a test to see if the differences are significant?� in that >> case, which one?. > > The models have different numbers of estimated parameters and they are not > nested. Many people claim that you can use neither deviance nor AIC (and they > are right). At least you must take into account the number of estimated > parameters that varies from zero to three. > > Cheers, Jari Oksanen > [[alternative HTML version deleted]]
_______________________________________________ R-sig-ecology mailing list R-sig-ecology@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology