많은 anovas를 실행하고 각 행에 대해 F 값을 얻으려고합니다.

Question

아래와 같은 데이터 세트로 작업하고 있습니다. 이제 저는 이것이 R이 좋아하는 형식의 형식이 아니라는 것을 압니다. 나는 데이터를 정리하는 방법을 알고 있지만, 각 unique_id에 대한 F 통계량 을 얻으려면 어떻게해야할지 잘 모릅니다. 이것이 내 목표입니다. 쉽게 할 수있는 방법이 있습니까? 그렇지 않으면, 각 행을 독립적으로 정리하고, 수행하고, anova를 적용하고, 새로운 통계로 F 통계를 추가하는 몇 가지 유형의 적용 함수를 사용할 수있는 방법이 있습니까?

unique_id heart heart heart kidney kidney kidney cortex cortex cortex 
373020.8 1.39 1.18 1.30 2.71 2.96 2.52 1.97 1.67 1.44 
371588.9 1.93 2.35 2.50 2.54 1.63 2.23 2.68 2.89 1.86 
367772.8 0.42 0.51 0.97 1.02 0.03 0.82 0.01 0.90 1.01 

Answer 1

나는 data.tables의 일부분이므로 데이터를 녹인 후 DT로 쉽게 할 수 있습니다. DT가 제공 한 정보가 들어있는 data.table과 함께 내 취해진 조치입니다.

DT <- dt<- melt(dt, 
     id.vars = c("unique_id"), 
     measure.vars = c("heart","cortex","kidney")) 

DT[,fstat:=summary(aov(value~variable))[[1]][1,"F value"],by=unique_id] 

이것은 F_Stat를 unique_id로 계산하므로 작동해야합니다.