大家好:
毕业论文需要用到R模拟项目反应理论的反应数据。这个是我编写的程序,运行之后可以在数据编辑器里查看到生成的数据矩阵,但是导出为txt文件之后,行和列的数字就不对了。请问怎样可以导出像数据编辑器里面一样的数据表啊(行和列正确的)?谢谢谢谢~~~~
附上程序部分:
I <- 1000
J <- 60
response <- matrix(0, I, J)
a <- rnorm(J, 0.8, 0.04)
b <- rnorm(J, 0, 1)
gamma<-rnorm(15,0,0.5)
theta <- rnorm(I, 0,1)
sdata <- matrix(0, I, J)
{
simulate <- function(I, J)
sdata <- matrix(0, I, J)
for( i in 1:I ) {
for( j in 1:30 ) {
temp <- a[j]*(theta[i]-b[j])
pij <- exp(temp)/(1+exp(temp))
ptemp <- runif(1)
if( pij >= ptemp) response[i,j] <- 1
sdata [i,j] <- response[i,j]
}
中间可以省略了
sdata<-data.frame
cat(sdata,file="redata.txt")


雷达卡



京公网安备 11010802022788号







