作者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