• Daum
  • |
  • 카페
  • |
  • 테이블
  • |
  • 메일
  • |
  • 카페앱 설치
 
카페정보
Biostatistics
 
 
 
 

친구 카페

 
 
카페 게시글
묻고 답하기 행(row)과 열(column) 출력하기 ( 수정 질문)
UUUUUU 추천 0 조회 164 15.03.10 14:15 댓글 2
게시글 본문내용
 
다음검색
댓글
  • 15.03.10 16:56

    첫댓글 많은 대안중에 가장 간단하고 속도도 손해보지 않는 방법은, theta1, theta2를 N*kitem matrix를 만드는 방법입니다.

    N<-1000 #row number
    kitem<-20 # col number

    theta1<-matrix(rnorm(N,0,1),nrow = N, ncol = kitem)
    theta2<-matrix(rnorm(N,-0.5,1),nrow = N, ncol = kitem)

    irf<-function(theta1,theta2, a1=.5, a2=1,b=0,cc=0.5)
    for(i in 1:N){
    for(i in 1:kitem){
    {
    d<- -b*sqrt(a1^2+a2^2)
    prob <- cc+(1-cc)/(1+exp(-((theta1*a1)+(theta2*a2)+d)))
    }
    return(prob)
    }
    }

    x=irf(theta1,theta2, a1=.5, a2=1,b=0,cc=0.5)
    dim(x)

  • 작성자 15.03.10 18:19

    죄송합니다. ,제가 우문을 했습니다. 행이 1000개, 열이 20개인 matrix 이고요, 각 값은 prob의 값이 발생되는 겁니다. 즉, 2만개의 값이 있는 겁니다.

최신목록