data1=data.frame(x1, x2, x3, x4, x5, x6, treat, outcome)
}
data1 这是所有代码 一直提示参数值意味着不同的行数: 1, 13, 12
data1=data.frame(x1, x2, x3, x4, x5, x6, treat, outcome)
}
data1 这是所有代码 一直提示参数值意味着不同的行数: 1, 13, 12
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.1x3[j]-2.0x4[j]-2.0x5[j]-0.4x6[j]))
treat[j]<-rbinom(1,1,pai1[j])
pai2[j]<-1/(1+exp(5-0.5x1[j]-0.5x2[j]-0.05x3[j]-x4[j]-x5[j]-0.2x6[j]))
outcome[j]<-rbinom(1,1,pai2[j])
data1=data.frame(x1, x2, x3, x4, x5, x6, treat, outcome)
}
data1 这是所有代码