Verwendung von Zeit-Tagebuch-Daten mit TraMineR

Question

Ich versuche, Sequenzanalyse mit Zeit-Tagebuch-Daten (American Time Use Survey) mit TraMineR in R zu tun. Ich habe die Daten als SPELL Daten (ID, Startzeit, Stoppzeit, Ereignis) aber ich die folgende Fehlermeldung, wenn es zu konvertieren versuchen, Daten STS oder SPS:

Error in as.matrix.data.frame(subset(data, , 2)) : dims [product 0] do not match the length of object [9]

ich glaube, es hat etwas damit zu tun, wie ich meine Zeit umwandeln (als Zeichen) Datum/Uhrzeit-Typen. Ich glaube, dass TraMineR ein POSIXlt-Format benötigt?

Hier ein Ausschnitt aus meinen Rohdaten (trcode ist das Ereignis)

head(atus.act.short)

 tucaseid tustarttim tustoptime trcode 
1 2.00701e+13 04:00:00 08:00:00 10101 
2 2.00701e+13 08:00:00 08:20:00 110101 
3 2.00701e+13 08:20:00 08:50:00 10201 
4 2.00701e+13 08:50:00 09:30:00 20102 
5 2.00701e+13 09:30:00 09:40:00 180201 
6 2.00701e+13 09:40:00 11:40:00 20102 

ich strptime verwenden, um die Zeichenkette zu POSIXlt konvertieren:

atus.act.short$starttime.new <- strptime(atus.act.short$tustarttim, format="%X") 
atus.act.short$stoptime.new <- strptime(atus.act.short$tustoptime, format="%X") 

Ich habe ich auch die ID bis zu nur zwei Ziffern

atus.act.short$id <- atus.act.short$tucaseid-20070101070000 

ich mit einem neuen Datenrahmen am Ende wie folgt:

id  starttime.new  stoptime.new trcode 
1 44 2012-08-03 04:00:00 2012-08-03 08:00:00 10101 
2 44 2012-08-03 08:00:00 2012-08-03 08:20:00 110101 
3 44 2012-08-03 08:20:00 2012-08-03 08:50:00 10201 
4 44 2012-08-03 08:50:00 2012-08-03 09:30:00 20102 
5 44 2012-08-03 09:30:00 2012-08-03 09:40:00 180201 
6 44 2012-08-03 09:40:00 2012-08-03 11:40:00 20102 
7 44 2012-08-03 11:40:00 2012-08-03 11:50:00 180201 
8 44 2012-08-03 11:50:00 2012-08-03 12:05:00 20102 
9 44 2012-08-03 12:05:00 2012-08-03 13:05:00 120303 
10 44 2012-08-03 13:05:00 2012-08-03 13:20:00 180704 
11 44 2012-08-03 13:20:00 2012-08-03 15:20:00 70104 
12 44 2012-08-03 15:20:00 2012-08-03 15:35:00 180704 
13 44 2012-08-03 15:35:00 2012-08-03 17:00:00 120303 
14 44 2012-08-03 17:00:00 2012-08-03 17:20:00 180701 
15 44 2012-08-03 17:20:00 2012-08-03 17:25:00 180701 
16 44 2012-08-03 17:25:00 2012-08-03 17:55:00 70101 
17 44 2012-08-03 17:55:00 2012-08-03 18:00:00 181203 
18 44 2012-08-03 18:00:00 2012-08-03 19:00:00 120303 
19 44 2012-08-03 19:00:00 2012-08-03 19:30:00 110101 
20 44 2012-08-03 19:30:00 2012-08-03 21:30:00 120303 
21 44 2012-08-03 21:30:00 2012-08-03 23:00:00 10101 
22 44 2012-08-03 23:00:00 2012-08-03 23:03:00 10201 
26 48 2012-08-03 06:45:00 2012-08-03 08:15:00 10201 
27 48 2012-08-03 08:15:00 2012-08-03 08:45:00 180209 
28 48 2012-08-03 08:45:00 2012-08-03 09:00:00 20902 
29 48 2012-08-03 09:00:00 2012-08-03 11:00:00 50101 
30 48 2012-08-03 11:00:00 2012-08-03 11:45:00 120312 

Dann versuche ich, ein Sequenzobjekt [mit Bibliothek (Traminer)]

atus.seq <- seqdef(atus.act.short, informat = "SPELL", id="id") 

Und ich folgende Fehlermeldung erhalten zu erstellen:

Error in as.matrix.data.frame(subset(data, , 2)) : dims [product 0] do not match the length of object [9]

Gedanken?

Answer 1

Ich habe es geschafft, dies zu umgehen, indem Sie die Zeit in Minuten umwandeln (nach weiteren Fragen zu stackoverflow), den Statuscode ein Zeichen (as.character), mit Seqformat, und zu einer Zeitachse zuweisen. Der neue Code lautet:

atus.seq2 <- seqformat(atus.act.short2, id="id", from="SPELL", to="STS", begin = "startmin", end = "stopmin", status="trcode", process = "FALSE")