Hi, I'm trying to set up AR(1) as a correlation structure in modeling some
data (attached file data.txt in text format) with lme, but have trouble
getting it to work.
Incent, Correctness, and Oppor are 3 categorical variables, Beta is a
response variable, and Time is an equally-spaced variable with 6 time points
(treated as a categorical variable as well). I'm not so sure about the form
in corAR1 here. Since the correlation is supposed to be along the Time
points, is form=~1 | Subject/Time a correct setup?
library(nlme)
Data <- read.table("data.txt", header=TRUE)
fm <- lme(Beta~Incent*Correctness+Oppor+(Time-1),
random = ~1|Subj, correlation=corAR1(0.5, form=~1|Subj/Time), Data)
The above lines do work, but my intuition seems pointing to a form like form
= ~time | Subj/Incent/Correctness/Oppor. However
fm <- lme(Beta~Incent*Correctness+Oppor+(Time-1),
random = ~1|Subj, correlation=corAR1(0.5,
form=~Time|/Incent/Correctness/Oppor), Data)
gives me the following error:
Error in attributes(.Data) <- c(attributes(.Data), attrib) :
'names' attribute [270] must be the same length as the vector [0]
It is probably due to the fact some combinations in the data are missing,
but I could not figure out how to group them so that the values on the time
variable would be unique within each subject. So how can I do this?
Any help is highly appreciated,
Gang
Subj Incent Oppor Correctness Time Beta
esp005 high noopp correct Time1 1.156180
esp005 high noopp correct Time2 -0.1524906
esp005 high noopp correct Time3 -0.1374741
esp005 high noopp correct Time4 0.1932847
esp005 high noopp correct Time5 -0.02334221
esp005 high noopp correct Time6 -0.4367026
esp005 high noopp incorr Time1 -1.046768
esp005 high noopp incorr Time2 0.7754745
esp005 high noopp incorr Time3 -0.4571449
esp005 high noopp incorr Time4 0.2027827
esp005 high noopp incorr Time5 0.1213075
esp005 high noopp incorr Time6 -0.9145339
esp005 high opp correct Time1 -0.09131422
esp005 high opp correct Time2 1.341441
esp005 high opp correct Time3 0.7011359
esp005 high opp correct Time4 -0.1799809
esp005 high opp correct Time5 0.465437
esp005 high opp correct Time6 -0.3150032
esp005 high opp incorr Time1 -0.8010447
esp005 high opp incorr Time2 0.7445193
esp005 high opp incorr Time3 -0.8163298
esp005 high opp incorr Time4 -0.7792278
esp005 high opp incorr Time5 -0.8525334
esp005 high opp incorr Time6 -0.9935786
esp005 low opp correct Time1 -1.135091
esp005 low opp correct Time2 2.080874
esp005 low opp correct Time3 0.9655223
esp005 low opp correct Time4 -0.4225955
esp005 low opp correct Time5 0.3518235
esp005 low opp correct Time6 0.7110413
esp005 low opp incorr Time1 -0.7773728
esp005 low opp incorr Time2 0.2896505
esp005 low opp incorr Time3 0.2957908
esp005 low opp incorr Time4 0.8173949
esp005 low opp incorr Time5 -0.8015972
esp005 low opp incorr Time6 -0.7611155
esp006 high noopp correct Time1 -0.2736552
esp006 high noopp correct Time2 0.0553967
esp006 high noopp correct Time3 0.06530985
esp006 high noopp correct Time4 -0.03726416
esp006 high noopp correct Time5 -0.4627595
esp006 high noopp correct Time6 -0.1016363
esp006 high noopp incorr Time1 -0.1492384
esp006 high noopp incorr Time2 0.02676032
esp006 high noopp incorr Time3 0.01891842
esp006 high noopp incorr Time4 0.2098106
esp006 high noopp incorr Time5 -0.07722018
esp006 high noopp incorr Time6 -0.0820671
esp006 high opp correct Time1 0.2699376
esp006 high opp correct Time2 0.4860634
esp006 high opp correct Time3 0.09422317
esp006 high opp correct Time4 0.1057343
esp006 high opp correct Time5 -0.09392627
esp006 high opp correct Time6 0.1700330
esp006 high opp incorr Time1 -0.2369925
esp006 high opp incorr Time2 -0.1076743
esp006 high opp incorr Time3 -0.1159856
esp006 high opp incorr Time4 0.07569676
esp006 high opp incorr Time5 -0.3941537
esp006 high opp incorr Time6 0.1117587
esp006 low opp correct Time1 0.05986514
esp006 low opp correct Time2 0.3611727
esp006 low opp correct Time3 0.1718891
esp006 low opp correct Time4 0.1112275
esp006 low opp correct Time5 0.08494787
esp006 low opp correct Time6 -0.08901607
esp006 low opp incorr Time1 0.1884523
esp006 low opp incorr Time2 0.2772646
esp006 low opp incorr Time3 0.01767265
esp006 low opp incorr Time4 -0.1886573
esp006 low opp incorr Time5 0.02343909
esp006 low opp incorr Time6 0.003311216
esp007 high noopp correct Time1 0.4830995
esp007 high noopp correct Time2 0.2181880
esp007 high noopp correct Time3 0.02222541
esp007 high noopp correct Time4 -0.06614286
esp007 high noopp correct Time5 -0.1391774
esp007 high noopp correct Time6 -0.1380095
esp007 high noopp incorr Time1 0.1973208
esp007 high noopp incorr Time2 0.2575948
esp007 high noopp incorr Time3 0.2585605
esp007 high noopp incorr Time4 0.2238915
esp007 high noopp incorr Time5 -0.3470621
esp007 high noopp incorr Time6 -0.3726538
esp007 high opp correct Time1 -0.08118145
esp007 high opp correct Time2 0.1611395
esp007 high opp correct Time3 -0.1699746
esp007 high opp correct Time4 -0.1834270
esp007 high opp correct Time5 -0.2152628
esp007 high opp correct Time6 -0.2895307
esp007 high opp incorr Time1 0.1846039
esp007 high opp incorr Time2 -0.02866489
esp007 high opp incorr Time3 0.5852502
esp007 high opp incorr Time4 -0.1285622
esp007 high opp incorr Time5 -0.1600584
esp007 high opp incorr Time6 -0.006998633
esp007 low opp correct Time1 -0.1103905
esp007 low opp correct Time2 -0.1153231
esp007 low opp correct Time3 -0.07760047
esp007 low opp correct Time4 -0.2746906
esp007 low opp correct Time5 -0.404328
esp007 low opp correct Time6 -0.1420324
esp007 low opp incorr Time1 -0.6065424
esp007 low opp incorr Time2 0.2828629
esp007 low opp incorr Time3 -0.009856104
esp007 low opp incorr Time4 -0.2008021
esp007 low opp incorr Time5 -0.09869743
esp007 low opp incorr Time6 -0.2824990
esp009 high noopp correct Time1 -1.505547
esp009 high noopp correct Time2 1.113987
esp009 high noopp correct Time3 0.4789144
esp009 high noopp correct Time4 -0.2463109
esp009 high noopp correct Time5 -0.2902244
esp009 high noopp correct Time6 0.7486134
esp009 high noopp incorr Time1 -0.1119943
esp009 high noopp incorr Time2 -0.8805183
esp009 high noopp incorr Time3 1.172378
esp009 high noopp incorr Time4 -0.005949792
esp009 high noopp incorr Time5 0.4770271
esp009 high noopp incorr Time6 -0.01373086
esp009 high opp correct Time1 -1.184459
esp009 high opp correct Time2 0.4254365
esp009 high opp correct Time3 0.7268725
esp009 high opp correct Time4 0.515629
esp009 high opp correct Time5 0.4637271
esp009 high opp correct Time6 0.4923486
esp009 high opp incorr Time1 -0.8852926
esp009 high opp incorr Time2 0.3829796
esp009 high opp incorr Time3 1.538067
esp009 high opp incorr Time4 3.054195
esp009 high opp incorr Time5 0.2178051
esp009 high opp incorr Time6 -0.4955013
esp009 low opp correct Time1 -0.6364805
esp009 low opp correct Time2 1.457549
esp009 low opp correct Time3 0.8776126
esp009 low opp correct Time4 1.062046
esp009 low opp correct Time5 0.201792
esp009 low opp correct Time6 1.594627
esp009 low opp incorr Time1 -3.359623
esp009 low opp incorr Time2 -0.3306261
esp009 low opp incorr Time3 2.607086
esp009 low opp incorr Time4 1.28359
esp009 low opp incorr Time5 1.469030
esp009 low opp incorr Time6 3.171498
esp011 high noopp correct Time1 0.3124615
esp011 high noopp correct Time2 -0.8804874
esp011 high noopp correct Time3 -0.2247847
esp011 high noopp correct Time4 -0.2671322
esp011 high noopp correct Time5 -0.7128651
esp011 high noopp correct Time6 0.01883687
esp011 high noopp incorr Time1 1.864378
esp011 high noopp incorr Time2 -0.750266
esp011 high noopp incorr Time3 -0.586842
esp011 high noopp incorr Time4 -0.6704107
esp011 high noopp incorr Time5 -0.3298926
esp011 high noopp incorr Time6 0.2421573
esp011 high opp correct Time1 -0.3177913
esp011 high opp correct Time2 -0.2981364
esp011 high opp correct Time3 -0.2146267
esp011 high opp correct Time4 -0.1804309
esp011 high opp correct Time5 0.08446378
esp011 high opp correct Time6 0.3543909
esp011 high opp incorr Time1 0.08388337
esp011 high opp incorr Time2 1.218138
esp011 high opp incorr Time3 -0.5016338
esp011 high opp incorr Time4 -0.1076762
esp011 high opp incorr Time5 0.5537836
esp011 high opp incorr Time6 -0.2860440
esp011 low opp correct Time1 -0.1256249
esp011 low opp correct Time2 -0.2261296
esp011 low opp correct Time3 -0.5473143
esp011 low opp correct Time4 -0.09333198
esp011 low opp correct Time5 -0.1071038
esp011 low opp correct Time6 0.02499052
esp011 low opp incorr Time1 -0.07223781
esp011 low opp incorr Time2 -0.149547
esp011 low opp incorr Time3 -0.1453610
esp011 low opp incorr Time4 -0.08941686
esp011 low opp incorr Time5 0.306246
esp011 low opp incorr Time6 0.2320469
esp012 high noopp correct Time1 -0.3615256
esp012 high noopp correct Time2 -0.3920005
esp012 high noopp correct Time3 -0.3296697
esp012 high noopp correct Time4 -0.01549745
esp012 high noopp correct Time5 -0.02441736
esp012 high noopp correct Time6 -0.2132724
esp012 high noopp incorr Time1 -0.373453
esp012 high noopp incorr Time2 -0.1007196
esp012 high noopp incorr Time3 0.1708301
esp012 high noopp incorr Time4 -0.01953903
esp012 high noopp incorr Time5 0.1978392
esp012 high noopp incorr Time6 -0.2226293
esp012 high opp correct Time1 -0.5573242
esp012 high opp correct Time2 -0.1668768
esp012 high opp correct Time3 -0.2498853
esp012 high opp correct Time4 0.02236604
esp012 high opp correct Time5 -0.001752809
esp012 high opp correct Time6 0.1219067
esp012 high opp incorr Time1 -0.1536199
esp012 high opp incorr Time2 -0.02398676
esp012 high opp incorr Time3 0.03128558
esp012 high opp incorr Time4 -0.5200846
esp012 high opp incorr Time5 -0.1566125
esp012 high opp incorr Time6 0.2605046
esp012 low opp correct Time1 -0.4457764
esp012 low opp correct Time2 -0.1480055
esp012 low opp correct Time3 0.0892735
esp012 low opp correct Time4 0.1253353
esp012 low opp correct Time5 0.160977
esp012 low opp correct Time6 0.06410025
esp012 low opp incorr Time1 0.01733564
esp012 low opp incorr Time2 -0.08083068
esp012 low opp incorr Time3 0.0009101268
esp012 low opp incorr Time4 0.1943673
esp012 low opp incorr Time5 0.1502581
esp012 low opp incorr Time6 0.2886754
esp014 high noopp correct Time1 0.3897176
esp014 high noopp correct Time2 0.2732953
esp014 high noopp correct Time3 -0.05613749
esp014 high noopp correct Time4 -0.08294953
esp014 high noopp correct Time5 0.03737959
esp014 high noopp correct Time6 0.1232754
esp014 high noopp incorr Time1 0.2681278
esp014 high noopp incorr Time2 0.3216812
esp014 high noopp incorr Time3 0.1409114
esp014 high noopp incorr Time4 -0.1725718
esp014 high noopp incorr Time5 -0.2724791
esp014 high noopp incorr Time6 -0.1426651
esp014 high opp correct Time1 0.4739094
esp014 high opp correct Time2 0.1792730
esp014 high opp correct Time3 0.176359
esp014 high opp correct Time4 0.3439640
esp014 high opp correct Time5 0.1924525
esp014 high opp correct Time6 0.1717151
esp014 high opp incorr Time1 0.3144462
esp014 high opp incorr Time2 0.2140982
esp014 high opp incorr Time3 0.09776947
esp014 high opp incorr Time4 0.08766148
esp014 high opp incorr Time5 -0.1158084
esp014 high opp incorr Time6 -0.03388017
esp014 low opp correct Time1 0.1512980
esp014 low opp correct Time2 0.4414964
esp014 low opp correct Time3 0.1622302
esp014 low opp correct Time4 0.1959245
esp014 low opp correct Time5 0.1135534
esp014 low opp correct Time6 -0.09733506
esp014 low opp incorr Time1 0.09931078
esp014 low opp incorr Time2 0.3388142
esp014 low opp incorr Time3 -0.2200165
esp014 low opp incorr Time4 -0.1145839
esp014 low opp incorr Time5 0.02053473
esp014 low opp incorr Time6 0.2068589
esp015 high noopp correct Time1 -0.3215265
esp015 high noopp correct Time2 0.4085156
esp015 high noopp correct Time3 0.2293433
esp015 high noopp correct Time4 -0.4682902
esp015 high noopp correct Time5 -1.071371
esp015 high noopp correct Time6 -0.7403038
esp015 high noopp incorr Time1 0.05301031
esp015 high noopp incorr Time2 -0.1548478
esp015 high noopp incorr Time3 -0.2736539
esp015 high noopp incorr Time4 -0.9517341
esp015 high noopp incorr Time5 -1.162436
esp015 high noopp incorr Time6 -0.6422004
esp015 high opp correct Time1 -0.2324042
esp015 high opp correct Time2 -0.3513996
esp015 high opp correct Time3 0.1431629
esp015 high opp correct Time4 -0.1607171
esp015 high opp correct Time5 -0.3969238
esp015 high opp correct Time6 -0.1688131
esp015 high opp incorr Time1 0.9312166
esp015 high opp incorr Time2 -1.121095
esp015 high opp incorr Time3 -0.1096821
esp015 high opp incorr Time4 0.04758765
esp015 high opp incorr Time5 -0.2127341
esp015 high opp incorr Time6 -0.5414439
esp015 low opp correct Time1 -0.2240375
esp015 low opp correct Time2 -0.08276297
esp015 low opp correct Time3 0.5523984
esp015 low opp correct Time4 -0.2775372
esp015 low opp correct Time5 -0.5579528
esp015 low opp correct Time6 -0.4111857
esp015 low opp incorr Time1 -0.2082284
esp015 low opp incorr Time2 0.1614645
esp015 low opp incorr Time3 -0.03589149
esp015 low opp incorr Time4 -0.768636
esp015 low opp incorr Time5 -0.8961188
esp015 low opp incorr Time6 -0.807752
esp019 high noopp correct Time1 -1.107367
esp019 high noopp correct Time2 0.6417076
esp019 high noopp correct Time3 1.025184
esp019 high noopp correct Time4 0.110842
esp019 high noopp correct Time5 -0.001747158
esp019 high noopp correct Time6 -0.2598766
esp019 high noopp incorr Time1 -0.2029920
esp019 high noopp incorr Time2 1.450233
esp019 high noopp incorr Time3 0.6189903
esp019 high noopp incorr Time4 0.6240945
esp019 high noopp incorr Time5 0.09637464
esp019 high noopp incorr Time6 -0.6620437
esp019 high opp correct Time1 -0.3997139
esp019 high opp correct Time2 -0.09644873
esp019 high opp correct Time3 0.6453347
esp019 high opp correct Time4 0.3290444
esp019 high opp correct Time5 0.003295310
esp019 high opp correct Time6 -0.1055445
esp019 high opp incorr Time1 -0.66997
esp019 high opp incorr Time2 1.186162
esp019 high opp incorr Time3 1.067688
esp019 high opp incorr Time4 0.7437585
esp019 high opp incorr Time5 0.3545609
esp019 high opp incorr Time6 0.1111318
esp019 low opp correct Time1 -1.113044
esp019 low opp correct Time2 -0.7529748
esp019 low opp correct Time3 -0.07419997
esp019 low opp correct Time4 -0.4987385
esp019 low opp correct Time5 -0.6503204
esp019 low opp correct Time6 -1.021259
esp019 low opp incorr Time1 -0.7887998
esp019 low opp incorr Time2 -0.48535
esp019 low opp incorr Time3 0.5645717
esp019 low opp incorr Time4 0.440514
esp019 low opp incorr Time5 -0.4387817
esp019 low opp incorr Time6 0.3583665
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
