atial(~ logIKA04 + nbLitre0 +
nbLitre1, cnm05g))#complete model
>sp05_nt <- update(sp05_01, trend = trend.spatial("cte", cnm05g))
#null model
>1-(sp07_nt$loglik/sp07_l0$loglik)^(2/590)#Maddala (1983)
[1] 0.0001510074
Thanks,
Regards,
Marion
--
Marion Jacquot
Laboratoire d
cients associated to each regressor : we
can get them with
**
>bivmod$coefficients
NBLITRE0 NBLITRE1
-0.001131662 -0.001083650
But I do not know what means the " \ " in the formula. I suppose it is
different from "/" which means "divide by ".
Any hint appreciated
ds,
Marion.
--
Marion Jacquot
Laboratoire de Chrono-environnement
UMR UFC/CNRS 6249 USC INRA
Université de Franche-Comté
Place Leclerc
F-25030 Besançon cedex
FRANCE
Tel. : +33 (0)381 665 829
Fax : +33 (0)381 665 797
http://chrono-environnement.univ-fcomte.fr/
__
3 matches
Mail list logo