setwd("C:\\Users\\user\\Desktop")
kkk <- read.csv("kkk.csv")
library(plyr)
kkk2 <- ddply(kkk, .(dq1,sq2,job,age),transform, df_m = mean(diff_sec))
kkk2$rl_dv <- with(kkk2, diff_sec - df_m)
write.csv(kkk2, "kkk3.csv", na="")
kkk5 <- ddply(kkk, .(dq1,age),transform, df_m = mean(diff_sec))
kkk5$rl_dv <- with(kkk5, diff_sec - df_m)
write.csv(kkk5, "kkk5.csv", na="")
ddply(kkk, .(dq1,age),summarize, df_m = mean(diff_sec))
## dq1 age df_m
## 1 1 1 837.0000
## 2 1 2 427.8000
## 3 1 3 437.5588
## 4 1 4 432.0297
## 5 1 5 465.5026
## 6 2 1 494.4444
## 7 2 2 375.8148
## 8 2 3 401.7500
## 9 2 4 412.9408
## 10 2 5 442.7431
## 11 3 1 319.3518
## 12 3 2 346.0890
## 13 3 3 367.7421
## 14 3 4 394.8369
## 15 3 5 434.1065
## 16 4 1 312.8865
## 17 4 2 324.8392
## 18 4 3 350.0650
## 19 4 4 368.9863
## 20 4 5 427.7339
## 21 5 1 274.1304
## 22 5 2 312.4023
## 23 5 3 351.4059
## 24 5 4 370.6220
## 25 5 5 438.3710
## 26 NA 1 396.0000
## 27 NA 2 338.6667
## 28 NA 3 341.0000
## 29 NA 4 478.9286
## 30 NA 5 534.5714
head(kkk2)
## X_FIRST id wt wt2 survey sid type sq1 sq1_1 sq1_2 sq1_3
## 1 8427 5177 0.479779 552.2500 2 5272 1 3 14 NA 198
## 2 6183 7262 0.594457 684.2500 2 7389 1 3 11 NA 170
## 3 1 7563 0.331176 381.2000 2 7696 1 3 15 NA 210
## 4 257 3669 0.701471 807.4286 2 3729 1 3 10 NA 156
## 5 257 3829 0.605099 696.5000 2 3889 1 3 14 NA 202
## 6 1395 3095 0.419530 482.9000 2 3144 1 3 16 NA 229
## sq1_4 sq2 sq3 q1_1 q1_1ori q1_1re q1_2 q1_2ori q1_2re q1_3 q1_3ori
## 1 1 1 24 2 2 NA 2 2 NA 5 9999
## 2 2 1 39 2 2 NA 4 4 NA 2 2
## 3 2 1 49 4 4 NA 4 4 NA 4 4
## 4 2 1 48 3 9999 3 3 3 NA 3 3
## 5 1 1 48 2 2 NA 2 2 NA 3 9999
## 6 2 1 49 3 9999 3 4 9999 4 2 9999
## q1_3re q2_1 q2_1ori q2_1re q2_2 q2_2ori q2_2re q2_3 q2_3ori q2_3re q2_4
## 1 5 3 9999 3 4 4 NA 4 4 NA 4
## 2 NA 2 2 NA 2 2 NA 2 2 NA 4
## 3 NA 4 4 NA 5 5 NA 4 4 NA 4
## 4 NA 2 2 NA 1 1 NA 1 1 NA 2
## 5 3 3 3 NA 3 3 NA 2 2 NA 2
## 6 2 5 5 NA 1 1 NA 2 2 NA 2
## q2_4ori q2_4re q2_5 q2_5ori q2_5re q2_6 q2_6ori q2_6re q2_7 q2_7ori
## 1 4 NA 1 1 NA 2 2 NA 4 4
## 2 4 NA 2 2 NA 5 5 NA 5 5
## 3 4 NA 4 4 NA 2 2 NA 2 2
## 4 2 NA 3 3 NA 3 3 NA 2 2
## 5 2 NA 3 3 NA 3 3 NA 3 3
## 6 2 NA 5 5 NA 5 5 NA 5 5
## q2_7re q2_8 q2_8ori q2_8re q2_9 q2_9ori q2_9re q2_10 q2_10ori q2_10re
## 1 NA 1 1 NA 2 2 NA 1 1 NA
## 2 NA 2 2 NA 4 4 NA 2 2 NA
## 3 NA 4 4 NA 2 2 NA 4 4 NA
## 4 NA 1 1 NA 3 3 NA 3 3 NA
## 5 NA 2 2 NA 2 2 NA 2 2 NA
## 6 NA 2 2 NA 2 2 NA 2 2 NA
## q2_11 q2_11ori q2_11re q2_12 q2_12ori q2_12re q2_13 q2_13ori q2_13re
## 1 4 4 NA 4 4 NA 4 4 NA
## 2 2 2 NA 2 2 NA 2 2 NA
## 3 4 4 NA 5 5 NA 4 4 NA
## 4 2 2 NA 2 2 NA 3 3 NA
## 5 2 2 NA 2 2 NA 1 1 NA
## 6 2 2 NA 4 4 NA 4 4 NA
## q2_14 q2_14ori q2_14re q2_15 q2_15ori q2_15re q2_16 q2_16ori q2_16re
## 1 4 4 NA 2 2 NA 5 5 NA
## 2 2 2 NA 2 2 NA 2 2 NA
## 3 4 4 NA 4 4 NA 4 4 NA
## 4 1 1 NA 2 2 NA 3 3 NA
## 5 2 2 NA 1 1 NA 4 4 NA
## 6 2 2 NA 4 4 NA 5 5 NA
## q3_1 q3_1ori q3_1re q3_2 q3_2ori q3_2re q3_3 q3_3ori q3_3re q3_4 q3_4ori
## 1 2 2 NA 4 4 NA 2 2 NA 4 4
## 2 2 2 NA 4 4 NA 2 2 NA 2 2
## 3 4 4 NA 4 4 NA 4 4 NA 4 4
## 4 4 4 NA 2 2 NA 2 2 NA 2 2
## 5 2 2 NA 2 2 NA 2 2 NA 2 2
## 6 2 2 NA 2 2 NA 1 1 NA 2 2
## q3_4re q3_5 q3_5ori q3_5re q3_6 q3_6ori q3_6re q3_7 q3_7ori q3_7re q3_8
## 1 NA 2 2 NA 4 4 NA 2 2 NA 3
## 2 NA 2 2 NA 2 2 NA 2 2 NA 2
## 3 NA 4 4 NA 5 5 NA 4 4 NA 4
## 4 NA 3 3 NA 1 1 NA 3 3 NA 2
## 5 NA 2 2 NA 3 3 NA 2 2 NA 2
## 6 NA 2 2 NA 2 2 NA 4 4 NA 4
## q3_8ori q3_8re q4 q4ori q4re q5 q5ori q5re q6 dq1 dq2 dq3 dq4 dq5 dq6 x
## 1 3 NA 2 2 NA 1 1 NA NA 1 1 1 3 1 1 1
## 2 2 NA 2 2 NA 3 3 NA 2 1 2 1 1 2 1 1
## 3 4 NA 4 4 NA 4 4 NA 1 1 1 1 3 1 4 1
## 4 2 NA 2 2 NA 3 3 NA NA 1 NA 1 2 2 1 1
## 5 2 NA 2 2 NA 3 3 NA NA 1 1 1 2 2 3 1
## 6 4 NA 2 2 NA NA 9999 9999 NA 1 1 1 2 1 4 1
## sex age edc job inc mar fam rel tq4 tq5 qq1_1 qq1_2 qq1_3 qq2_1 qq2_2
## 1 1 1 1 1 1 3 1 1 1 1 4 4 1 3 2
## 2 1 2 1 1 1 1 2 1 1 2 4 2 4 4 4
## 3 1 3 1 1 1 3 1 5 3 3 2 2 2 2 1
## 4 1 3 1 1 NA 2 2 1 1 2 3 3 3 4 5
## 5 1 3 1 1 1 2 2 3 1 2 4 4 3 3 3
## 6 1 3 1 1 1 2 1 5 1 NA 3 2 4 1 5
## qq2_3 qq2_4 qq2_5 qq2_6 qq2_7 qq2_8 qq2_9 qq2_10 qq2_11 qq2_12 qq2_13
## 1 2 2 5 4 2 5 4 5 2 2 2
## 2 4 2 4 1 1 4 2 4 4 4 4
## 3 2 2 2 4 4 2 4 2 2 1 2
## 4 5 4 3 3 4 5 3 3 4 4 3
## 5 4 4 3 3 3 4 4 4 4 4 5
## 6 4 4 1 1 1 4 4 4 4 2 2
## qq2_14 qq2_15 qq2_16 qq3_1 qq3_2 qq3_3 qq3_4 qq3_5 qq3_6 qq3_7 qq3_8 qq4
## 1 2 4 1 4 2 4 2 4 2 4 3 4
## 2 4 4 4 4 2 4 4 4 4 4 4 4
## 3 2 2 2 2 2 2 2 2 1 2 2 2
## 4 5 4 3 2 4 4 4 3 5 3 4 4
## 5 4 5 2 4 4 4 4 4 3 4 4 4
## 6 4 2 1 4 4 5 4 4 4 2 2 4
## qq5 qq6 y_q1_1 y_q1_2 y_q1_3 y_q2_1 y_q2_2 y_q2_3 y_q2_4 y_q2_5 y_q2_6
## 1 1 NA 75 75 0 50 25 25 25 100 75
## 2 3 3 75 25 75 75 75 75 25 75 0
## 3 4 2 25 25 25 25 0 25 25 25 75
## 4 3 NA 50 50 50 75 100 100 75 50 50
## 5 3 NA 75 75 50 50 50 75 75 50 50
## 6 NA NA 50 25 75 0 100 75 75 0 0
## y_q2_7 y_q2_8 y_q2_9 y_q2_10 y_q2_11 y_q2_12 y_q2_13 y_q2_14 y_q2_15
## 1 25 100 75 100 25 25 25 25 75
## 2 0 75 25 75 75 75 75 75 75
## 3 75 25 75 25 25 0 25 25 25
## 4 75 100 50 50 75 75 50 100 75
## 5 50 75 75 75 75 75 100 75 100
## 6 0 75 75 75 75 25 25 75 25
## y_q2_16 y_q3_1 y_q3_2 y_q3_3 y_q3_4 y_q3_5 y_q3_6 y_q3_7 y_q3_8 y_q4
## 1 0 75 25 75 25 75 25 75 50 75
## 2 75 75 25 75 75 75 75 75 75 75
## 3 25 25 25 25 25 25 0 25 25 25
## 4 50 25 75 75 75 50 100 50 75 75
## 5 25 75 75 75 75 75 50 75 75 75
## 6 0 75 75 100 75 75 75 25 25 75
## y_q5 yq1_1 yq1_2 yq1_3 yq2_1 yq2_2 yq2_3 yq2_4 yq2_5 yq2_6 yq2_7 yq2_8
## 1 0 2.5 2.5 5.5 NA 5.5 5.5 5.5 2.5 2.5 5.5 2.5
## 2 50 2.5 5.5 2.5 2.5 2.5 2.5 5.5 2.5 5.5 5.5 2.5
## 3 75 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 2.5 2.5 5.5
## 4 50 NA NA NA 2.5 2.5 2.5 2.5 NA NA 2.5 2.5
## 5 50 2.5 2.5 NA NA NA 2.5 2.5 NA NA NA 2.5
## 6 NA NA 5.5 2.5 5.5 2.5 2.5 2.5 5.5 5.5 5.5 2.5
## yq2_9 yq2_10 yq2_11 yq2_12 yq2_13 yq2_14 yq2_15 yq2_16 yq3_1 yq3_2 yq3_3
## 1 2.5 2.5 5.5 5.5 5.5 5.5 2.5 5.5 2.5 5.5 2.5
## 2 5.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 5.5 2.5
## 3 2.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5 5.5
## 4 NA NA 2.5 2.5 NA 2.5 2.5 NA 5.5 2.5 2.5
## 5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 5.5 2.5 2.5 2.5
## 6 2.5 2.5 2.5 5.5 5.5 2.5 5.5 5.5 2.5 2.5 2.5
## yq3_4 yq3_5 yq3_6 yq3_7 yq3_8 yq4 yq5 phonetype intvwerID
## 1 5.5 2.5 5.5 2.5 NA 2.5 2.5 2 34709
## 2 2.5 2.5 2.5 2.5 2.5 2.5 NA 2 34531
## 3 5.5 5.5 5.5 5.5 5.5 5.5 5.5 2 24568
## 4 2.5 NA 2.5 NA 2.5 2.5 NA 2 25394
## 5 2.5 2.5 NA 2.5 2.5 2.5 NA 2 25394
## 6 2.5 2.5 2.5 5.5 5.5 2.5 NA 2 29248
## StartTime EndTime diff1
## 1 24NOV13:16:28:00 24NOV13:16:43:00 1900-01-00 00:14:40
## 2 29NOV13:17:03:00 29NOV13:17:08:00 1900-01-00 00:04:21
## 3 30NOV13:11:27:00 30NOV13:11:34:00 1900-01-00 00:06:55
## 4 17NOV13:13:27:00 17NOV13:13:34:00 1900-01-00 00:06:41
## 5 17NOV13:14:45:00 17NOV13:14:53:00 1900-01-00 00:07:38
## 6 15NOV13:16:42:00 15NOV13:16:52:00 1900-01-00 00:09:15
## diff2 year month day hour min sec diff_sec contact cont1
## 1 1900-01-00 00:14:40 1900 1 0 0 14 40 880 1 1
## 2 1900-01-00 00:04:21 1900 1 0 0 4 21 261 1 1
## 3 1900-01-00 00:06:55 1900 1 0 0 6 55 415 1 1
## 4 1900-01-00 00:06:41 1900 1 0 0 6 41 401 1 1
## 5 1900-01-00 00:07:38 1900 1 0 0 7 38 458 1 1
## 6 1900-01-00 00:09:15 1900 1 0 0 9 15 555 1 1
## cont2 cont3 cont4 cont5 TIME missing LAGID LAGTIME T SEQNUM
## 1 NA NA NA NA 900 1 34709 1700927580 8427 82
## 2 NA NA NA NA 300 0 34531 1701360600 6183 363
## 3 NA NA NA NA 420 0 24568 1701376260 1 178
## 4 NA NA NA NA 420 2 25394 1700313540 257 108
## 5 NA NA NA NA 480 1 25394 1700317020 257 112
## 6 NA NA NA NA 600 2 29248 1700151000 1395 151
## df_m rl_dv
## 1 880.0000 0.000000
## 2 261.0000 0.000000
## 3 451.2857 -36.285714
## 4 451.2857 -50.285714
## 5 451.2857 6.714286
## 6 451.2857 103.714286