小编Udd*_*din的帖子

我试图在 R markdown 中运行以下代码以打印代码。我不想打印 R 代码的结果。代码中的某些行确实很大。当我运行时，它工作得很好，但是有些行超出了页面？这是 R markdown 代码

title: R code
output: 
  pdf_document:
    toc: true
    number_sections: true
documentclass: article
classoption: 
  - portrait
  - a4paper
---



```{r,eval=F,echo=TRUE}
hrs_data_GS$rSLn==i,\]%*%beta\[k,\])^2/2*sigma2)+
                          sum(z_ij_str\[i\]*xdata\[hrs_data_GS$rSLn==i,\]%*%alpha\[k,\])-
                          0.5*sum(w_ij\[hrs_data_GS$rSLn==i\]*(xdata\[hrs_data_GS$rSLn==i,\]%*%alpha\[k,\])^2)+
                          ( (lamda/sigma2)*sum(hrs_data_GS$y_ij\[hrs_data_GS$z_ij==1 & hrs_data_GS$rSLn==i\]-
                                                 xdata\[hrs_data_GS$z_ij==1 & hrs_data_GS$rSLn==i,\]%*%beta\[k,\])
                            +sum_z_ij_str\[i\]- sum(w_ij\[hrs_data_GS$rSLn==i\]*
                                                    xdata\[hrs_data_GS$rSLn==i,\]%*%alpha\[k,\]))^2*
                          2*((n_i\[i\]*lamda^2/sigma2)+sum_w_ij_j\[i\]+1/sig2.theta)^-1)+
          (log(pi_cls\[i,k\]))  #(log(pi_k\[k\]))  
      }
    }
    nume_prod<-nume_prod-apply(nume_prod,1,max) #it works 
    nume_prod <- exp(nume_prod)
    prob_cls_rsum=apply(nume_prod, 1, sum) 
    prob_cls=nume_prod/ prob_cls_rsum 
    #cls=miscF::rMultinom(p=cbind(prob_cls)) #miscF is not available at CRC
    cls=Hmisc::rMultinom(p=cbind(prob_cls),1)  
    #Update sum_y.ij_j and sum_wxa_ij_j before updating the theta
    #these values involve 'cls'
    #******
    for (i in …

我想从具有三个值（例如 ）的多项分布生成大小为 20 的样本1,2 and 3。例如，样本可以是这样的sam=(1,2,2,2,2,3,1,1,1,3,3,3,2,1,2,3,...1)

下面的代码可以工作，但没有得到预期的结果

> rmultinom(20,3,c(0.4,0.3,0.3))+1
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20]
[1,]    1    1    3    2    2    1    1    2    3     2     3     2     1     2     2     3     1     2     2     2
[2,]    2    1    2    1    3    2    4    2    1     2     2     1     1     2     1     2     3     2     3     3
[3,]    3    4    1    3    1 …

我想创建18名为mat10ci, mat15ci,....,的相同维度的空矩阵mat95ci。我试过下面的代码

for( s in seq(10,95,by=5)){
  paste0("mat",s,"ci")=array(NA, dim = c(10, 20))
}

我有这个错误

target of assignment expands to non-language object

帮助表示赞赏。

对于以下数据

> data
   x1 x2 x3
1   1  0  0
2   0  1  0
3   0  0  1
4   1  0  1
5   1  0  0
6   0  1  1
7   1  0  0
8   0  1  0
9   1  1  1
10  1  0  0

如何识别所有列中具有相同值的行。对于data，行号4 ，并6没有列相同的值1 & 3，并2 & 3分别和行号9对所有列相同的值。如何在R中标识这些行号。不胜感激。