2016-09-09 54 views
0

可有人請與下面的情景幫助嗎?我對SaS非常陌生,我不確定如何讓它起作用?的SaS情景生成

模擬從下面的線性模型200點意見:

Y = alpha + beta1 * X1 + beta2 * X2 + noise 

其中:

alpha=1, beta1=2, beta2=-1.5

X1 ~ N(1, 4), X2 ~ N(3,1), noise ~ N(0,1)

我曾嘗試這個代碼,但不知道它的完全準確:

DATA ONE; 
alpha = 1; 
beta1 = 2; 
beta2 = -1.5; 
RUN; 

DATA CALC; 
    SET ONE; 
    DO i = 1 to 200; 
Y=alpha+beta1*X1+beta2*X2+Noise; 
X1=Rannor(1); 
X2=rannor(3); 
Noise=ranuni(0); 
    OUTPUT; 
    END; 

RUN; 

PROC PRINT DATA=CALC; 
RUN; 

回答

1

您必須在SAS幫助下找了主題 「rannor」,「ranuni」,「生成隨機數」,......

rannor:生成標準正態分佈的房車。

ranuni:統一分佈式房車。

在rannor的參數是種子數量,而不是預期的值。

如果N(X,Y)在你的例子是指隨機變量通常與預期值x和標準偏差Y分佈(或做你指的是方差???)那麼代碼可以是(看看在改變順序的發言上; Y的定義必須是隨機數的定義後...):

DATA ONE; 
alpha = 1; 
beta1 = 2; 
beta2 = -1.5; 
RUN; 

DATA CALC; 
SET ONE; 
seed = 1234; 
DO i = 1 to 200; 
    X1=1+4*Rannor(seed); 
    X2=3+rannor(seed); 
    Noise=rannor(seed); 
    Y=alpha+beta1*X1+beta2*X2+Noise; 
    OUTPUT; 
END; 
RUN; 

PROC PRINT DATA=CALC; 
RUN; 

有用於生成隨機數也變體,例如「打電話給rannor」。在SAS中處理種子數量有不同的概念。有關這些主題,請參閱SAS幫助。 here