1

で反復測定のロジスティック回帰を実行したいと思いRます。落ち込んでいる子供と落ち込んでいない子供で、その日の最後の食事の時間が異なるかどうかを確認したいと思います。私はすべての被験者（dep / non-dep）について14日間の食事時間を記録しています。私が抱えている唯一の問題は、これらの時間変数を適切な形式に変換して分析を実行することです。それらを10進数（例：15.5）に変換するのは良い考えではないようです。助けてください！

id  depressed mealtime     
B8         1 17:30:00       
B8         1 17:00:00      
B8         1 12:30:00      
B8         1     <NA>         
B8         1 19:45:00       
B8         1 19:30:00       
A1         0 19:30:00       
A1         0 18:45:00      
A1         0 19:30:00      
A1         0 18:30:00       
A1         0 20:30:00
4

1 に答える 1

1

参照点を使用して、食事時間変数を時差に変換してみませんか? たとえばstrptime、文字列を強制するために使用するPOSIXltと、次のdifftimeようなことができます。

dat$mealtime <- strptime(dat$mealtime,'%H:%M:%S')
dat$difference <- difftime(dat$mealtime,time2=strptime('00:00:00','%H:%M:%S'))

これで、新しく作成された変数を回帰に使用できます（glmロジットにはあなたを想定しています）

fit <- glm(depressed ~ difference,data=dat, family=binomial("logit"))

PS：ここにデータがあります：

dat <- read.table(text='id  depressed mealtime     
B8         1 17:30:00       
B8         1 17:00:00      
B8         1 12:30:00      
B8         1     <NA>         
B8         1 19:45:00       
B8         1 19:30:00       
A1         0 19:30:00       
A1         0 18:45:00      
A1         0 19:30:00      
A1         0 18:30:00       
A1         0 20:30:00',header=TRUE)
于 2013-03-16T08:37:58.303 に答える