プログラムの書き方についてです。

質問

ひまわりの種が２０００年のと２００６年ながあるとします。
その平均値の差tを求めるプログラムはどう書けばいいでしょうか？

t=|平均(!2006)ー平均(!2000)|/√(n(!2006)-1)s^2(!2006)+(n(!2000)-1)s^2(!2000)/n(!2006)+(n(!2000)-2(*1/n(!2006)+1/n(!2006))

分子は|平均(!2000)ー平均(!2006)|です。あとはルートの中に含まれています。
(*1/n(!2006)+1/n(!2006))の部分はルート中の分母にかけるのではないです。そのままの掛け算です。
不偏分散も平均のプログラムもわかる状態なのでサブルーチンで出ことはできますが、この式はどう書けばいいのかわかりません。

教えてください。わかりにくかったら補足していきたいので教えてください。

fortranでお願いします。

回答No.3

noname#227064 · ベストアンサー

> なぜそう解釈されてしまうのですか？

演算子には優先順位があります。
参考URLを見ると、+, -より*, /の優先順位が高くなっています。
したがって、
a+b/c+d
という式は、
a+(b/c)+d
と評価されます。

参考URL：

http://www.hiroshima-cu.ac.jp/japanese/IPC/hunet99/sun/WorkShop/ja/html_docs/fortran/f77rm/expressions.doc.html#1345

回答No.2

noname#227064

fortranはわからないので実行されない理由についてはアドバイスできませんが、それ以外で気になった点がありましたので回答します。

> t=abs(mean2-mean1)/sqrt(float(n-1)*unvar2+float(n-1)*unvar1/float(n)+float(n)-2)*(1/float(n)+1/float(n))

この式は、

t=(abs(mean2-mean1)/sqrt(float(n-1)*unvar2+(float(n-1)*unvar1/float(n))+float(n)-2))*(1/float(n)+1/float(n))

と解釈されるのではないでしょうか。

補足

なぜそう解釈されてしまうのですか？

教えてください

回答No.1

Wr5

>自分でいろいろやった結果がここでの質問となりました。

そのいろいろやった内容を掲示して、それに対するアドバイスなどをもらうようにしないと丸投げと言われても仕方ありません。
本当にあなたが「いろいろやった」のか、ネットの向こうに居る人間には判断できません。
ラクをする為に「いろいろやった」ことにしようとしているかも知れません。
そんな訳で「努力の結果」を掲示して質問してください。

補足

t=abs(mean2-mean1)/sqrt(float(n-1)*unvar2+float(n-1)*unvar1/float(n)+float(n)-2)*(1/float(n)+1/float(n))

私はこのように書いたのですが実行されません。

unvar1,2は2000年と2006年の不偏分散を。
mean1,2はそれぞれの平均です。これはサブルーチンでcallしています。
floatの部分がおかしいのだと思うのですが、どうしたらよいのかわかりません。
実行されないのはなにがいけないのでしょうか？