People, many thanks for keeping this list so functional. Sorry for the
slow reply, had a small accident, but is ok now. I'll answer Liam
first and Ted and Brian in sequence.
Liam, yes, the error is fatal, and the BrownieREML worked. But seems
it did not converged:
These are the results:
Hope not being abusing of this list.
Possibly, I missed something. I tried both with matrix and dataframe and
reached to the same problem:
phy.manova(tree, c(table$option,table$difoption), table$morph,
data.names=NULL, nsim=1000, test=Wilks)Warning: no tip labels, order
assumed to be the
Hi Agus.
I think the problem is that c() concatenates your data for the
continuous dependent variables in your model into one long vector. Try
instead:
table-as.matrix(table)
Y-cbind(table[,2],table[,3])
x-as.factor(table[,1])
phy.manova(tree,Y,x)
Good luck. - Liam
Liam J. Revell,