Displaying 2 results from an estimated 2 matches for "numerat".
Did you mean:
numeral
2009 Mar 28
1
Error in R??
...eciprocal condition number = 0
In addition: There were 50 or more warnings (use warnings() to see the first
50)
The R code is part of a bigger program.
##sample from full conditional distribution of Si
#Prob(Si = 1)
for (j in 1:N)
{
numerat =
currphi1*exp((-1/(2*currvar))*t(Y1-(t(XX1)%*%currbeta1))%*%(Y1-(t(XX1)%*%currbeta1)))
denomin =
currphi2*exp((-1/(2*currvar))*t(Y2-(t(XX2)%*%currbeta2))%*%(Y2-(t(XX2)%*%currbeta2)))
sum=denomin + numerat
ProbSi = numerat/sum
arunofSi[j]=rbinom(1,1,P...
2007 May 02
1
Degrees of freedom in repeated measures glmmPQL
...says Groups:
Subjects = 106, Treatment = 3.
Should I be worried or is this what to expect?!
I was expecting it to be more like an ANOVA table, where the error
degrees of freedom should reflect the number of subjects, not all the
data points.
I can't see the usual differentiation between the numerater and
denominator/error degrees of freedom, so am I right in thinking the DF
column shows the error degrees of freedom? Or do glmms not work like
this?
Thank you very much in advance,
Charlotte