Re: [問題] X(ij)=M+b(i)+e(ij) i=1...I,j=1...J

看板Fortran作者sjgau (sjgau)時間16年前 (2008/09/16 06:48)推噓0(0推 0噓 1→)

留言1則, 1人參與討論串2/2 (看更多)

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