R程序求助

我用T=10没有问题

1. set.seed(100)
   
2. library(MASS)
   
3. library(binarySimCLF)
   
4. 
   
5. N=2100
   
6. T=10
   
7. 
   
8. x1=list()
   
9. for(i in 1:N){x1[[i]]=runif(T)}
   
10. 
    
11. x2=list()
    
12. for(i in 1:N){x2[[i]]=runif(T)}
    
13. 
    
14. x3=list()
    
15. for(i in 1:N){x3[[i]]=runif(T)}
    
16. 
    
17. x4=list()
    
18. for(i in 1:N){x4[[i]]=runif(T)}
    
19. 
    
20. x5=list()
    
21. for(i in 1:N){x5[[i]]=runif(T)}
    
22. 
    
23. X1=unlist(x1)
    
24. X2=unlist(x2)
    
25. X3=unlist(x3)
    
26. X4=unlist(x4)
    
27. X5=unlist(x5)
    
28. 
    
29. f1=1.5*X1^2-mean(1.5*X1^2)
    
30. f2=2*sin(2*pi * X2)-mean(2*sin(2*pi * X2))
    
31. f3=.8*cos(2*pi*X3)+X3^.5-mean(.8*cos(2*pi*X3)+X3^.5)
    
32. f4=-1.5*sin(2*pi*X4)+X4^2-mean(-1.5*sin(2*pi*X4)+X4^2)
    
33. f5=1.5*sqrt(X5)+cos(2*pi*X5)-mean(1.5*sqrt(X5)+cos(2*pi*X5))
    
34. f=f1+f2+f3+f4+f5
    
35. 
    
36. # generating correlated discrete data
    
37. library(binarySimCLF)
    
38. 
    
39. y = matrix(rep(NA,T*N),N,byrow=TRUE)
    
40. f0=matrix(f,N,byrow=TRUE)
    
41. mu =exp(f0)/(1+exp(f0))# mean matrix
    
42. for (i in 1:N)
    
43. {
    
44. rho=(rhoRange(mu[i,])$rhomax+rhoRange(mu[i,])$rhomin)/2
    
45. temp = ranXch(mu[i, ], rho)
    
46. if (temp$succeed)
    
47.     y[i,] = temp$y
    
48. }

复制代码

谢谢，非常谢谢 qioqpwqr 同学，希望以后有机会多多交流。