data(tortoise)
A<-tortoise$T + tortoise$F.med.high
Re(eigen(A)$values[1])
y<-numeric(4)
for (i in 2:5)
{
A<-tortoise$T+tortoise[[i]]
y[i-1]<-Re(eigen(A)$values[1])
}
data.frame(row.names=names(tortoise)[2:5], lambda=y)
Run the code above in your browser using DataLab