Hi all,
I have been trying to run a phylogenetic model using a recently updated
version of MCMCglmm, but for some reason the model will not run when using
cbind to form the response variable. I have two variables, one for number
of individuals uninfected (NoUninf) and one for number infected (NoInf),
and I have tried both separately to make sure the model works. However,
when I combine the two using cbind, the model cannot assign the R
structure, even if I make sure to assign the family using c("gaussian,
"gaussian") to account for the two variables in the response (error
repeatedly says r structure not specified properly). I'm not sure if the
length of the variable is confusing the program, as cbind works perfectly
well in a glm using lme4.
On a separate note, I'm not sure how to combine rows to make sure that the
number of records in the data file equals the number of records in the
phylogeny (1) for each species. Basically, I have a data file with multiple
rows for each species depending on sampling location, ie:
Species Site NoInf NoUninf
C. canutus WA 10 25
C. canutus NA 11 34
Since each line is important as it represents a unique sample location, I
would like to incorporate them all into the data file. However, my
phylogeny and covariance matrix will only have one record per species, so
should I try to shape my data so that each row is a matrix of multiple row
records? Any help would be greatly appreciated, sorry if these are rather
mundane problems but I'm quite new to using R.
Thanks,
Nick
--
Nick Clark
PhD Candidate, School of Environment
Griffith University, Gold Coast
Researchgate profile:
https://www.researchgate.net/profile/Nicholas_Clark4/?ev=hdr_xprf
[[alternative HTML version deleted]]