作者hyekyo0608 (魚)
看板Fortran
標題[問題] X(ij)=M+b(i)+e(ij) i=1...I,j=1...J
時間Mon Sep 15 23:34:20 2008
我要做的是單因子隨機效應的模型
X(ij)=M+b(i)+e(ij)
where X(ij)is the jth observation from ith batch,
the b(i) and e(ij) are mutually independent normally distributed random
variables with E(bi)=E(eij)=0 and variances b^2 and w^2
我要用fortran產生X(ij)
OUTTER:DO NI=1,I
CALL DRNNOR(NR,B)
M(NI)=B*(DSORT(R))
INNTER1:DO NJ=1,Jn
CALL DRNNOR(NR,S)
E(NI,NJ)=S
X(NI,NJ)=M(NI)+E(NI,NJ)
END DO ONNTER1
END DO OUTTER
我想請問一下這個寫法是在產生X(ij)嗎?
可是我看不懂CALL DRNNOR(NR,B)
是說叫出NR個N(0,1)? 那B的意思是什麼?
CALL DRNNOR(NR,S)
E(NI,NJ)=S
一樣S的部份不懂?
然後為啥E(NI,NJ)=S?
最後最後X(NI,NJ)=M(NI)+E(NI,NJ)是不是就是我要的模型X(ij)=M+b(i)+e(ij)??
拜託拜託
好心的FORTRAN 高手們
幫我解決我的疑難,我困擾好久了
拜託拜託!!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.31.105.145
1F：推 andachi:call 是呼叫函式或是副程式阿 這是很基本的指令喔 09/16 00:47
2F：→ andachi:DRNNOR(NR,S(B?))是說依次序呼叫這兩個變數近副程式或函式 09/16 00:49
3F：→ DrStein:這程式不完整吧...請po出全部 謝謝 09/16 01:37