Я работаю с примером PTW в R, и он работает хорошо. Однако, когда я пытаюсь заменить свои данные данными в этом примере, я получаю ошибки. Мои данные отформатированы так же, как и данные примера, поэтому я не знаю, что я делаю неправильно.
Это работает:
data(gaschrom)
ref <- gaschrom[1:8,]
samp <- gaschrom[9:16,]
gaschrom.ptw <- ptw(ref, samp, warp.type = "individual",
optim.crit = "RMS", init.coef = c(0, 1, 0, 0))
plot(gaschrom.ptw)
plot(gaschrom.ptw, what = "signal")
Однако, когда я загружаю свои собственные данные (ниже), онивыходит из строя. Мои данные 3 наблюдения с 123 переменными. R распознает, что когда я импортирую данные в RStudio b
-13, -20,637, -19, -11,2.604,10,44,57,327,44,10,2.313, -11, -20, -6.545,1,52,70.082,57,17,18.593,0,5,4.754,4, -3, -1,204, -1,10,26.964,25,43,57.435,40,17,12.874, -2,-16, -27,441, -23, -24, -31,418, -20, -11, -13,463, -6, -3, -4,643, -2, -2, -4,643, -3, -6,7.727,10,44,53.049,40, -3, -14,796, -24, -24, -31,708, -20, -11, -14,598, -7, -7, -10,321, -7, -7, -10,321,-7, -7,6.591,10,44,57.327,44,10,7.727, -6, -3,6.056,8,30,38.152,29,5, -2,474, -11, -19, -30,548,-23, -24, -31,418, -20, -11, -14,598, -7, -7,6.591,10,44,57.327,44,10,8.388, -5, -1, -3,134, -1,-5, -12.802, -11, -20, -27.431, -20, -11, -8.966, -1,10,12.209,10 -2,7.554, -2,40,49.062,40, -3, -14,796, -24, -24, -29,912, -18, -6,18.559,24,73,93.119,69,12,23.458,3,41,51.912,41,2, -0,353, -11,1,15.803, 18,41,47.903,33, -5, -12,677, -20, -11, -18,876, -11, -20, -28,646, -21, -15, -23,738, -15, -21, -25,161, -16, -1,2.826,7,7,12.736,9,14,11.417,6, -14, -24,368, -22, -19, -28,863, -19, -21, -33,117, -23,-24, -35,696, -24, -24, -26,076, -14,6,3.654,6, -14, -9,964, -8,24,28.462,24, -8, -14,461, -19, -8, -8,92, -2, -2, -7,212, -6, -13, -23,74, -19, -22, -18,166, -8,24,38.372,34,22,25.179,11,2, -2,715, -6, -13, -2.55,2,41,66.592,56,46,49.651,24, -8, -9,964, -14,6,3.654,6, -14, -20,663, -19 -19, -27,154, -17, -17, -22,032, -13, -6, -10,917, -6, -13, -22,032, -17, -17, -22,032, -13, -6, -10,917, -6, -13, -22,032, -17, -17, -22,032, -13, -6, -10,917, -6, -13, -2.702,2,41,55.285,45,13,16.614,2,11,14.481,11,2, -6,924, -10, -26, -23,842, -15,24,36.204,34,16,27.907,15,32,46.134,34,22, 21.475,7, -9, -17,535, -17, -17, -22,032, -13, -6, -4.863,0,5,5.889,5,0, -4,863, -6, -13, -22,032, -17, -17, -25,736, -17, -17, -25,736, -17, -17, -25,736, -17, -17, -25,736, -17, -17, -25,736, -17, -17, -25,736, -17, -17, -25,736, -17, -17, -25,736, -17, -17, -25,736, -17, -17, -25,736, -17, -17, -25,736,-17, -17, -25,736, -17, -17, -25,736, -17, -17, -6,407,2,41,50,162,39, -2, -13,788, -23
Myданные 3 наблюдения с 123 переменными. R распознает, что когда я импортирую данные в RStudio, но когда я пытаюсь запустить анализ следующим образом, он не работает.
train <- read.csv(file="threesamps2.csv", header=FALSE, sep=",")#this works, the data is brought into "train" and shows to have 3 observations with 123 variables.
ref <- train[1,]# This works and shows 1 obs, 123 variables
samp <- train[2:3,]# This works and shows 2 obs, 123 variables
train.ptw <- ptw(ref, samp, warp.type = "individual",
optim.crit = "RMS", init.coef = c(0, 1, 0, 0))#This seems to be what is failing. Why? It works fine with the sample data.
plot(train.ptw)
plot(train.ptw, what = "signal")
И ref, и sample загружаются правильно, но я получаю ошибку «Ошибка вplot (train.ptw, what = "signal"): объект train.ptw не найден. Если я правильно читаю код, необходимо создать файл train.ptw, но он не найден. Любые идеи будут оценены.