2011-11-19 8 views
11

Wiem, że zostało to zadane wiele razy tutaj pod rubryką "od długiego do szerokiego", ale natknąłem się na sytuację, w której mam dwie zmienne wartości, które są powtarzanymi pomiarami.Długie do szerokich w/dwa powtórzone pomiary

  id sex time  score1 score2 
1 subject 1 m Time1 -0.20926263 0.2499310 
2 subject 2 m Time1 0.17147511 3.2708905 
3 subject 3 m Time1 -0.82619584 0.5993917 
4 subject 4 f Time1 -0.95568823 4.4729726 
5 subject 5 f Time1 -2.29939525 8.0101254 
6 subject 1 m Time2 -0.37914702 3.6387589 
7 subject 2 m Time2 0.26759909 4.9027533 
8 subject 3 m Time2 0.07727621 2.1848642 
9 subject 4 f Time2 -0.08613439 5.8747074 
10 subject 5 f Time2 -0.02743044 4.3963938 
11 subject 1 m Time3 0.07176053 3.7959496 
12 subject 2 m Time3 0.46463917 5.2494579 
13 subject 3 m Time3 -0.68764512 2.2639503 
14 subject 4 f Time3 -0.56670061 2.3361909 
15 subject 5 f Time3 1.70731774 5.8345116 

Szybki sposób na odtworzenie ramki danych (DF).

DF<-data.frame(id=rep(paste("subject", 1:5, sep=" "), 3),    
       sex=rep(c("m","m","m","f","f"), 3),      
       time=c(rep("Time1",5), rep("Time2",5), rep("Time3",5)), 
       score1=rnorm(15), score2=abs(rnorm(15)*4))    

mogę rozwiązać kwestię długo szerokości dwóch zmierzonych powtarzanych zmiennych pomiarowych wykorzystujących funkcję reshape od podstawy ale miałem nadzieję na plyr lub reshape2/1 odpowiedzi, jak te pakiety są na ogół znacznie bardziej intuicyjny do mnie. Jeśli masz jakieś inne rozwiązania, przekaż je, ponieważ nauka będzie świetna.

Rozwiązanie z bazy:

wide <- reshape(DF, v.names=c("score1", "score2"), idvar="id",   
       timevar="time", direction="wide")      
wide 
+1

? Nie widzę, jakie było zamknięcie ("poza tematem") dla ...? –

+1

Ben przybił go ze zmianą kształtu i tak małym kodem. Jestem bardzo otwarty na inne metody (nie jestem pewien, czy/jak to zrobić z data.table itp.) –

Odpowiedz

12

Myślę, że będzie to zrobić:

library(reshape) 
m <- melt(DF) 

najprostszy, ale czas i wynik są w odwrotnej kolejności od Twojego przykład (w przypadku jest to ważne)

cast(m,id+sex~...) 

Lub bardziej bezpośrednio:

cast(m,id+sex~variable+time) 

Można to ograniczyć do jednego-liner:

recast(DF,id+sex~...) 

Jeśli chcesz, możesz użyć nowszej reshape2 pakiet zamiast reshape, zastępując cast z dcast (od wersji recast zawartych w reshape2 robi nie daje pożądanego rezultatu.)

+0

Doskonałość. Dzięki za możliwość uczenia się. –

2

Odpowiedź Bena jest wspaniała, biorąc pod uwagę datę, w której zadano to pytanie. Należy jednak zauważyć, że funkcja dcast zawarta w nowszych wersjach "data.table" może poradzić sobie z takimi problemami, bez konieczności wcześniejszego przesyłania danych do postaci melt. Jako takie byłoby bardziej efektywne podejście do korzystania.

library(data.table) 
dcast(as.data.table(DF), id + sex ~ time, value.var = c("score1", "score2")) 
##   id sex score1_Time1 score1_Time2 score1_Time3 score2_Time1 score2_Time2 score2_Time3 
## 1: subject 1 m 0.78213630 -0.1557955 -0.10278773 1.5771598  1.013447 1.4583278 
## 2: subject 2 m 0.07456498 -1.4707524 0.38767161 0.2372536  2.787854 3.0741317 
## 3: subject 3 m -1.98935170 -0.4781501 -0.05380504 4.4001015  2.226653 0.4493848 
## 4: subject 4 f 0.61982575 0.4179416 -1.37705956 3.0527030  2.755023 3.5244309 
## 5: subject 5 f -0.05612874 1.3586796 -0.41499456 0.6580944  2.829981 1.5924235 
Powiązane problemy