r语言参数值意味着不同的行数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

data1=NULL
x1<- rnorm(1,0,1)
x2<- rlnorm(1,0,0.5)
x3<- rnorm(1,0,10)
x4<- rbinom(1,1,0.5)
x5<- rbinom(1,1,0.2)
x6<-c(rep(1,6),rep(2,6),rep(3,))
x7<-vector(length=12)
x8<-vector(length=12)
x9<-vector(length=12)
x10<-vector(length=12)
pai1<-vector(length=12)
pai2<-vector(length=12)
treat<-vector(length=12)
outcome<-vector(length=12)

for  (j in 1:12)
{



  x1[j]<- rnorm(1,0,1)
  x2[j]<- rlnorm(1,0,0.5)
  x3[j]<- rnorm(1,0,10)
  x4[j]<- rbinom(1,1,0.5)
  x5[j]<- rbinom(1,1,0.2)

  x7[j] <-sin(x1[j])
  x8[j] <- x2[j]^2
  x9[j] <- x3[j]*x4[j]
  x10[j] <- x4[j]*x5[j]


pai1[j]<-1/(1+exp(3.5-x1[j]-x2[j]-0.1*x3[j]-2.0*x4[j]-2.0*x5[j]-0.4*x6[j]))
treat[j]<-rbinom(1,1,pai1[j])

pai2[j]<-1/(1+exp(5-0.5*x1[j]-0.5*x2[j]-0.05*x3[j]-x4[j]-x5[j]-0.2*x6[j]))
outcome[j]<-rbinom(1,1,pai2[j])



data1=data.frame(x1, x2, x3, x4, x5, x6, treat, outcome)  

}

data1 这是所有代码 一直提示参数值意味着不同的行数: 1, 13, 12
> data1

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

data1=NULL

x7<-vector(length=12)
x8<-vector(length=12)
x9<-vector(length=12)
x10<-vector(length=12)
pai1<-vector(length=12)
pai2<-vector(length=12)
treat<-vector(length=12)
outcome<-vector(length=12)

for  (j in 1:12)
{
  

  
  x1[j]<- rnorm(1,0,1)
  x2[j]<- rlnorm(1,0,0.5)
  x3[j]<- rnorm(1,0,10)
  x4[j]<- rbinom(1,1,0.5)
  x5[j]<- rbinom(1,1,0.2)
  
  x7[j] <-sin(x1[j])
  x8[j] <- x2[j]^2
  x9[j] <- x3[j]*x4[j]
  x10[j] <- x4[j]*x5[j]
  
  
pai1[j]<-1/(1+exp(3.5-x1[j]-x2[j]-0.1*x3[j]-2.0*x4[j]-2.0*x5[j]-0.4*x6[j]))
treat[j]<-rbinom(1,1,pai1[j])

pai2[j]<-1/(1+exp(5-0.5*x1[j]-0.5*x2[j]-0.05*x3[j]-x4[j]-x5[j]-0.2*x6[j]))
outcome[j]<-rbinom(1,1,pai2[j])
  


data1=data.frame(x1, x2, x3, x4, x5, x6, treat, outcome)  

}

data1  本来的代码是这个 但是这个代码会提示找不到x1，所以我加了x1<- rnorm(1,0,1)
x2<- rlnorm(1,0,0.5)
x3<- rnorm(1,0,10)
x4<- rbinom(1,1,0.5)
x5<- rbinom(1,1,0.2)  但是又提示参数不对 请各位老师指出我的问题 谢谢

data1=NULL

x7<-vector(length=12)
x8<-vector(length=12)
x9<-vector(length=12)
x10<-vector(length=12)
pai1<-vector(length=12)
pai2<-vector(length=12)
treat<-vector(length=12)
outcome<-vector(length=12)

for  (j in 1:12)
{
  

  
  x1[j]<- rnorm(1,0,1)
  x2[j]<- rlnorm(1,0,0.5)
  x3[j]<- rnorm(1,0,10)
  x4[j]<- rbinom(1,1,0.5)
  x5[j]<- rbinom(1,1,0.2)
  
  x7[j] <-sin(x1[j])
  x8[j] <- x2[j]^2
  x9[j] <- x3[j]*x4[j]
  x10[j] <- x4[j]*x5[j]
  
  
pai1[j]<-1/(1+exp(3.5-x1[j]-x2[j]-0.1*x3[j]-2.0*x4[j]-2.0*x5[j]-0.4*x6[j]))
treat[j]<-rbinom(1,1,pai1[j])

pai2[j]<-1/(1+exp(5-0.5*x1[j]-0.5*x2[j]-0.05*x3[j]-x4[j]-x5[j]-0.2*x6[j]))
outcome[j]<-rbinom(1,1,pai2[j])
  


data1=data.frame(x1, x2, x3, x4, x5, x6, treat, outcome)  

}

data1  本来的代码是这个 但是这个代码会提示找不到x1，所以我加了x1<- rnorm(1,0,1)
x2<- rlnorm(1,0,0.5)
x3<- rnorm(1,0,10)
x4<- rbinom(1,1,0.5)
x5<- rbinom(1,1,0.2)  但是又提示参数不对 请各位老师指出我的问题 谢谢