Я хочу использовать «pic» для филогенетически независимых контрастов с 196 переменными в моем фрейме данных.Я могу сделать это один за другим, используя следующий код.Как я могу сделать это в меньшем количестве строк кода (применимо, для цикла или других методов)?Ниже приведен код для одной переменной.Заранее спасибо.С наилучшими пожеланиями Zillur
prot=pic(f1[,197],t1) # f1=dataframe, t1=phylogenetic tree
p1=pic(f1[,1],t1)
cor(prot,p1)