[問題] X(ij)=M+b(i)+e(ij) i=1...I,j=1...J
我要做的是單因子隨機效應的模型
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
推
09/16 00:47, , 1F
09/16 00:47, 1F
→
09/16 00:49, , 2F
09/16 00:49, 2F
→
09/16 01:37, , 3F
09/16 01:37, 3F
討論串 (同標題文章)
以下文章回應了本文：
完整討論串 (本文為第 1 之 2 篇)：
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
