Группировка переменных между различными таблицами данных в R - PullRequest
0 голосов
/ 09 июля 2019

Я хотел бы выполнить регрессию с помощью этой формы:

Мой код выглядит следующим образом:

R_j_k <- read.table(text = "
    j     R_j_1  R_j_2     R_j_3 
     1      20    15         40 
     2      45    30         50 
     3      60    45         60 
     4      75    60         60 
     5      90    70         60 
     6     105    70         90 
     7     120    70        120 
     8     125    70        140 
     9     130    70        160 
    10     145    70        180 
           ",header = T)

R_i_k <-read.table(text = "
             i     R_i_1  R_i_2    R_i_3 
             1     150    70        190 
             2     155    70        200 
             3     150    75        195 
             4     160    80        190 
             5     150    75        180 
               ",header = T)

Я хотел бы оценить параметр b, и я не знаю, как группировать переменные из разных таблиц данных вместе.Итак, я хотел бы научиться назначать IV для моего DV.

Моя попытка:

Sum_R_j_k <- data.frame(
    j = seq(1:10),
    k_1 = rnorm(10,0,0),  
    k_2 = rnorm(10,0,0),
    k_3 = rnorm(10,0,0)
)

Sum_R_j_k[1,2] <-  R_j_k[1,2]
Sum_R_j_k[1,3] <-  R_j_k[1,3]
Sum_R_j_k[1,4] <-  R_j_k[1,4]

    for (k in 2:4) {
    for (j in 2:10) {
        Sum_R_j_k[j,k] <- Sum_R_j_k[j-1,k] + R_j_k[j,k]

    }}


  Sum_R_j_k <- read.table(text = " 
        j        k_1    k_2      k_3
    1   20  15  40
    2   65  45  90
    3   125 90  150
    4   200 150 210
    5   290 220 270
    6   395 290 360
    7   515 360 480
    8   640 430 620
    9   770 500 780
    10  915 570 960
        ",header = T)

for (i in 1:5) {
for (k in 2:4) {
for (j in 1:10) {

   nls(R_i_k[i,k] ~ R_j_k[j,k]*b^j)

    }}}

Это не работает.Я думаю, что мне нужно сгруппировать переменные, а не использовать для циклов.К сожалению, я понятия не имею, как.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...