Ошибка при использовании для l oop в блестящей приборной панели R

Question

Я пытаюсь создать блестящий на основе ранее созданного кода R.

Проблема, с которой я столкнулся, заключается в том, что циклы for не работают.

Одним из примеров может быть :

    
    centros.acp=reactive({colMeans(a.acp())})
    
    lcov.acp=reactive({solve(cov(a.acp()))})
    
    dm.acp=reactive({rep(0,length(a.acp()[,1]))})
    
    na.acp=reactive({a.acp()[1]})
    
    dm.acp()=reactive({
      value=dm.acp()
      for(i in na.acp()){

      value()[i]<-reactive({round(t(a.acp()[i,]-centros.acp())%*%lcov.acp()%*%(a.acp()[i,]-centros.acp()),3)

    })}})
    T2<-reactive({dm.acp()})
    my_list=reactive({list("T2" = T2(),"cc.sw.UCL.99" = cc.sw.UCL.99(),"cc.sw.UCL.95" = cc.sw.UCL.95(), "num.com"=num.com())})
   
   output$plot_t2=renderPlot({plot(1:n(),my_list()$T2, ylim = c(0,max(my_list()$T2,my_list()$cc.sw.UCL.99)+3),main = "",
         xlab = 'OBS',ylab = 'T2') })

Я получаю сообщение об ошибке

Error in dm.acp() = reactive({ : invalid (NULL) left side of assignment

Код, который я показываю, - единственная часть, которая имеет значение для проблемы, все остальное работает. Я пробовал разные способы, но всегда получаю ошибку или построенный график основан на исходном dm.acp (), который состоит из группы 0. Когда я запускаю для l oop, я хочу заменить 0 значениями ( на каждой итерации меняется новый номер).

Исходный dm.acp () равен

> dm.acp
 [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[50] 0 0 0 0 0 0 0 0 0 0 0

, а после for l oop должно быть

[1]  6.210  3.663  4.246  4.279  3.766  3.583  2.957  4.796  4.578  3.225  4.232  2.865  6.976  3.454
[15]  3.411  3.530  3.610 22.747  9.936  2.659  3.080  6.412  9.183  2.781  3.253  5.585  3.183  5.171
[29]  3.898  5.868  2.498 12.901  2.772  2.795  1.644  1.393 21.078  2.489  3.198  2.412  4.213  4.251
[43]  8.144  3.201  3.295  4.378  4.746  9.481  4.281  4.402  5.329  2.637  2.627  4.507  5.521  3.545
[57]  2.623  5.835  2.438  3.230

Я ценю любую помощь. Спасибо!

Answer 1

Попробуйте заменить for l oop на lapply, как указано ниже:

lapply(1:na.acp(), function(i){
  value()[i]<-reactive({
    round(t(a.acp()[i,]-centros.acp())%*%lcov.acp()%*%(a.acp()[i,]-centros.acp()),3)
  })
})

Обновление: Попробуйте этот код.

dm.acp<-c(rep(0,60))
    
    n <- length(dm.acp)
    lapply(1:n, function(i) { 
      dm.acp[i] <<- i+2
    })

Вы должны иметь возможность реализовать это в своей программе. Если вы собираетесь использовать функцию в lapply, заключите его как lapply(1:(length(na.acp)),function(i){…})

Update2: См. Код ниже

###  replace this code  ###
dm.acp=reactive({rep(0,length(a.acp()[,1]))})

na.acp=reactive({a.acp()[1]})

dm.acp()=reactive({
  value=dm.acp()
  for(i in na.acp()){
    
    value()[i]<-reactive({round(t(a.acp()[i,]-centros.acp())%*%lcov.acp()%*%(a.acp()[i,]-centros.acp()),3)
      
    })}})
T2<-reactive({dm.acp()})


###  with this code   ###
dm.acpv=reactive({rep(0,length(a.acp()[,1]))})

na.acp=reactive({a.acp()[1]})

T2 <- reactive({
  value=dm.acpv()
  
  lapply(1:(na.acp()), function(i){
    value[i] <<- round(t(a.acp()[i,]-centros.acp())%*%lcov.acp()%*%(a.acp()[i,]-centros.acp()),3)
  })
  value
})
## no need to define T2 again as I renamed dm.acp as T2