SAS/IML : RANDNORMAL에서 개별 분산 구성 요소를 사용하는 방법

Question

이것은 프로그래밍 관련 질문이지만 통계 배경을 먼저 설명해 드리겠습니다. 이 질문은 혼합 효과 위치 스케일 모델 (즉, 이질적 분산)에 대한 데이터 시뮬레이션의 일부를 지칭합니다. IML의 RANDNORMAL 함수를 사용하여 두 MVN 분산 구성 요소를 시뮬레이션하려고합니다. 두 분산 구성 요소가 이기종이므로 ​​RANDNORMAL에서 사용하는 분산은 사람마다 다릅니다. 따라서 IML을 사용하여 특정 행 (예 : 행 1 = 사람 1)을 선택하고 다음 행으로 이동하기 전에 RANDNORMAL 함수를 사용하는 등의 작업이 필요합니다.

아래 예제 코드는 2 인분입니다. DO를 사용하여 각 개인의 특정 분산 구성 요소 (VC1 및 VC2)를 반복합니다. 오류가 발생합니다 : "모듈 RANDNORMAL 이전 호출에서 종료하기 전에 다시 전화." DO 루프에서 BREAK 또는 EXIT 함수가 필요하다고 가정하고 있지만 아무 것도 시도하지 않았습니다.

PROC IML; 
    ColNames = {"ID" "VC1" "VC2"}; 
    A = {1 2 3, 
     2 8 9}; 
    PRINT A[COLNAME=ColNames]; 
    /*Set men of each variance component to 0*/ 
    MeanVector = {0, 0}; 
    /*Loop through each person's data using THEIR OWN variances*/ 
    DO i = 1 TO 2; 
     VC1 = A[i,2]; 
     VC2 = A[i,3]; 
     CovMatrix = {VC1 0, 
        0 VC2}; 
     CALL RANDSEED(1); 
     U = RANDNORMAL(2, MeanVector, CovMatrix); 
    END; 
QUIT; 

도움을 주시면 감사하겠습니다. 아, SAS 9.4를 사용하고 있습니다.

Answer 1

당신은 주위를 움직이기를 원하지만 대개는 U를 두 번 다시 쓰고 싶지 않습니다. 당신이하려고하는 것을 이해한다면 U의 첫 번째 행을 쓰고 U의 두 번째 행을 써야합니다. 다음은 루프를 통해 매번 새로운 것을 생성하는 것보다 (그리고 느리다) U와 _Cv 행렬을 j() 이후에 좀 더 효율적입니다.

proc iml; 

    a = {1 2 3,2 8 9}; 
    print(a); 
    _mv = {0,0}; 
    U = J(2,2); 
    _cv = J(2,2,0); 
    CALL RANDSEED(1); 

    do i = 1 to 2; 
    _cv[1,1] = a[i,2]; 
    _cv[2,2] = a[i,3]; 
    U[i,] = randnormal(1,_mv, _cv); 
    end; 
    print(u); 
quit; 

Answer 2

귀하의 실수는 유효한 SAS/IML 구문 아니다 라인

CovMatrix = {VC1 0, 0 VC2}; /* wrong */ 

입니다. "Constructing block matrices with applications to mixed models." 대신, 기사 "How to build matrices from expressions."

을 당신은 또한 블록 대각 행렬과 시뮬레이션을 수행하는 방법을 설명하는이 문서에 관심이있을 수 있습니다 참조 조의 접근 @ 사용하거나 자세한 내용은

CovMatrix = (VC1 || 0) // (0 || VC2); 

를 사용