Блог

Итак, я создал дерево без указания длины ветвей, и у меня есть набор данных с точно такими же именами наконечников и 5 измерений для каждого из этих видов в общей сложности (450 измерений). Я пытаюсь выполнить MLE (оценку максимального правдоподобия), используя эти измерения, чтобы изменить длину ветви, но я не уверен, как это сделать. Мне нужно использовать все 5, и длины ветвей, которые мне нужны, относятся к прямым ветвям к образцу. Короткий путь состоит в том, чтобы определить значения MLE для каждого образца на основе 5 измерений, которые все они имеют, и ввести их в Mesquites.exe но еще раз код, который у меня есть, позволяет мне взглянуть только на одну из переменных, а не на все 5. 5 переменных, которые у меня есть, — это 5 различных измерений длины, но мне нужно использовать их все, а не только одну

 wingL <- geodata$wingL
MLreconstruction <- ace(wingL, geotree, type="continuous", method="ML")
MLreconstruction

 

Любая помощь будет признательна