Number of Groups in SpeciesMix
Dear friends, I am willing to apply the SAM analytical framework to a dataset of plant species in coastal Brazil using the SpeciesMix package. The SpeciesMix package fits Species Archtype Models, a special type of finite mixture of regression model motivated by the analysis of multi-species data. In appying function clusterSelect, which helps in defining the best number of species groups G, I would like to confirm if my understanding is correct: the formula reported there as "obs ~ 1 + x" in clusters <- clusterSelect(obs~1+x,dat1$pa,dat,G=2:5,em.refit=2) is a generic formulation not directly related to the species (dat1$pa) or environmental (dat) data, isn't it? So in principle I should use this same formulation as well, understanding that obs stands for the whole species data matrix, 1 for the presence of a constant, and x for the whole explanatory dataset? I tried to apply obs~1+x but it returns an error message, however. I am kind of blocked here so any thoughts could help... Sincerely, Alexandre
Dr. Alexandre F. Souza Professor Adjunto III Universidade Federal do Rio Grande do Norte CB, Departamento de Ecologia Campus Universit?rio - Lagoa Nova 59072-970 - Natal, RN - Brasil lattes: lattes.cnpq.br/7844758818522706 http://www.docente.ufrn.br/alexsouza [[alternative HTML version deleted]]