Arhn - архитектура программирования

Объединение вложенных выходных данных цикла for в вектор

Я подогнал многоэтапную модель Маркова (MSM), используя 3 ковариаты. Я пытаюсь применить функцию, основными аргументами которой являются подобранная модель МСМ и значения ковариат. Функция:

transition<-qmatrix.msm(MSMmod,  ci="normal", covariates=list(grossTon= 10, activ=0, period=1))

Выходом этой функции являются оценочные значения:

              Owning<10 Owning10-40   Owning>40        left
Owning<10   -0.18037446  0.06140559  0.00000000 0.118968868
Owning10-40  0.01609661 -0.09027454  0.04502546 0.029152476
Owning>40    0.00000000  0.04801757 -0.05137377 0.003356196
left         0.00000000  0.00000000  0.00000000 0.000000000

а меня интересуют только некоторые из них, так что с transition$estimates[c(5,2,10,7,13,14,15)] получаю; например, значения, которые дает функция qmatrix.msm с комбинацией ковариат GrossTon=120, active=0 и period=1 (см. ниже).

[1] 0.06140559 0.01609661 0.04502546 0.04801757 0.118968868 0.029152476 0.003356196

С обоими циклами for в моей функции я ожидаю получить 48 векторов (по 7 значений в каждом, аналогично предыдущему) в результате всех комбинаций GrosTon (12 значений) и activ (4 значения), таким образом, наконец, чтобы объединить все из них в единый вектор (336 значений, 7*12*4).

Это моя функция:

transRate<-function(period){
  estim<-data.frame(matrix(rep(0,336),336,1))

  for(i in seq(10,120,by=10)){
    for(j in seq(0,3, by=1)){

    estim[c(i,j)]<-qmatrix.msm(msm.Mult4,  ci="normal", covariates=list(grossTon=i, activ=j, period=period))$estimates[c(5,2,10,7,13,14,15)]
    outp[c(i,j)]<-c(estim[c(i,j)])#Here I'm trying to get my 336 values vector
    }
  }

  grosTvect<-sort(rep(seq(10,120,by=10),28))
  rate<-rep(c("q12","q21","q23","q32","q14","q24","q34"),48)
  estimRate<-data.frame(grosTvect,rate,outp)

  return(estimRate)                     
} 
dataFrame<-transRate(period=1)

Я не знаю, как обращаться с i и j, чтобы получить свой вектор.

Я получаю следующую ошибку, когда пишу estim[c(i,j)] и outp[c(i,j)]<-c(estim[c(i,j)])

error in `[<-.data.frame`(`*tmp*`, c(i, j), value = c(0.0614055886960195,  : 
  new columns would leave holes after existing columns

И это другое, когда я пишу estim[i,j] и outp<-estim[i,j]

Error in `[<-.data.frame`(`*tmp*`, i, j, value = c(0.0614055886960195,  : 
  replacement has 7 rows, data has 1

Любая помощь будет оценена.

r for-loop
31.01.2013

  • Рафаэль, у вас есть хорошее начало для вопроса, но не могли бы вы добавить больше информации о векторе, который должна создавать ваша функция? И вы получаете какую-то ошибку с вашим текущим кодом? 31.01.2013
  • Надеюсь, теперь стало понятнее. 31.01.2013

Ответы:


1

Ваш фрейм данных "estim" имеет неправильные размеры для хранения вектора длины 7 из "оценок". Вы также не объявляете «outp», поэтому, если код интерпретатора когда-либо попадет туда после успешного выполнения назначения кадра данных, это вызовет еще одну ошибку. (В любом случае неясно, почему вы предпринимаете дополнительный шаг по назначению «outp».)

Было бы проще создать фрейм данных из 48 строк со столбцами с именами c("q12","q21","q23","q32","q14","q24","q34") и просто назначить правильную строку = i+12*j с использованием последовательных индексов?

 transRate<-function(period){
  estim<-data.frame(q12=1:48, q21=0, q23=0, q32=0, q14=0, q24=0, q34=0)))

  for(i in seq(1,12)){
    for(j in seq(0,3, by=1)){

    estim[ i+12*j , ] <- qmatrix.msm(cav.msm,  ci="normal", 
              covariates=list( grossTon = i*10, activ=j, period=period))$
                                                  estimates[c(5,2,10,7,13,14,15)]
       }
  }

    return(estim)                     
} 

dataFrame<-transRate(period=1)

(Не мог понять, как этот объект «grosTvect» с разными размерами должен был совпадать с результатами «оценок».)

31.01.2013

2

кадр данных estim имеет только один столбец. Этот столбец является матрицей. Когда вы делаете что-то вроде estim[i,j], вы получаете доступ к столбцам фрейма данных. Я полагаю, что вы действительно хотите получить доступ к матрице внутри фрейма данных. Итак, вам придется использовать что-то вроде estim[1][,c(i,j)], если вы хотите заполнить матрицу.

31.01.2013
Новые материалы
Коллекции публикаций по глубокому обучению
Последние пару месяцев я создавал коллекции последних академических публикаций по различным подполям глубокого обучения в моем блоге https://amundtveit.com - эта публикация дает обзор 25..
Представляем: Pepita
Фреймворк JavaScript с открытым исходным кодом Я знаю, что недостатка в фреймворках JavaScript нет. Но я просто не мог остановиться. Я хотел написать что-то сам, со своими собственными..
Советы по коду Laravel #2
1-) Найти // You can specify the columns you need // in when you use the find method on a model User::find(‘id’, [‘email’,’name’]); // You can increment or decrement // a field in..
Работа с временными рядами спутниковых изображений, часть 3 (аналитика данных)
Анализ временных рядов спутниковых изображений для данных наблюдений за большой Землей (arXiv) Автор: Рольф Симоэс , Жильберто Камара , Жильберто Кейрос , Фелипе Соуза , Педро Р. Андраде ,..
3 способа решить квадратное уравнение (3-й мой любимый) -
1. Методом факторизации — 2. Используя квадратичную формулу — 3. Заполнив квадрат — Давайте поймем это, решив это простое уравнение: Мы пытаемся сделать LHS,..
Создание VR-миров с A-Frame
Виртуальная реальность (и дополненная реальность) стали главными модными терминами в образовательных технологиях. С недорогими VR-гарнитурами, такими как Google Cardboard , и использованием..
Демистификация рекурсии
КОДЕКС Демистификация рекурсии Упрощенная концепция ошеломляющей О чем весь этот шум? Рекурсия, кажется, единственная тема, от которой у каждого начинающего студента-информатика..