作者sjgau (sjgau)
看板Fortran
標題Re: [問題] X(ij)=M+b(i)+e(ij) i=1...I,j=1...J
時間Tue Sep 16 06:48:45 2008
http://tinyurl.com/5td6ee
RNNOR
Generates pseudorandom numbers from a standard normal distribution using an
inverse CDF method.
Required Arguments
R — Vector of length NR containing the random standard normal deviates.
(Output)
Optional Arguments
NR — Number of random numbers to generate. (Input)
Default: NR = size (R,1).
FORTRAN 90 Interface
Generic: CALL RNNOR (R [,…])
Specific: The specific interface names are
S_RNNOR and D_RNNOR.
FORTRAN 77 Interface
Single: CALL RNNOR (NR, R)
Double: The double precision name is DRNNOR.
Description
※ 引述《hyekyo0608 (魚)》之銘言：
: 我要做的是單因子隨機效應的模型
: 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: 118.231.129.227