版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

本帖產(chǎn)生 1 個(gè) 程序強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

ian_zhangty

木蟲(chóng) (著名寫(xiě)手)

程序強(qiáng)帖: 1
應(yīng)助: 1 (幼兒園)
金幣: 4358.4
散金: 63
紅花: 3
帖子: 1043
在線: 86.6小時(shí)
蟲(chóng)號(hào): 296621
注冊(cè): 2006-11-12
專業(yè): 植物學(xué)研究的新技術(shù)、新方

[求助] 請(qǐng)問(wèn)，我想用fortran計(jì)算統(tǒng)計(jì)中的p值以及95%信度空間

請(qǐng)問(wèn)有沒(méi)有現(xiàn)成的程序，謝謝

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有9人回復(fù)
293求調(diào)劑已經(jīng)有14人回復(fù)
0703化學(xué)336分求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有9人回復(fù)
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請(qǐng) 已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

請(qǐng)問(wèn)這種“粘手”玩具的具體材料是什么？如何制備成大塊的軟膠墊？已經(jīng)有10人回復(fù)
化工計(jì)算基礎(chǔ) 已經(jīng)有468人回復(fù)
用MS中的Forcite計(jì)算冰，結(jié)果發(fā)現(xiàn)溫度波動(dòng)異常的大，感覺(jué)像做周期性的震蕩衰減已經(jīng)有6人回復(fù)
請(qǐng)問(wèn)對(duì)于生產(chǎn)制劑的原料藥產(chǎn)地變更的申請(qǐng)，其變更說(shuō)明該怎么寫(xiě) 已經(jīng)有8人回復(fù)
【求助】Fortran語(yǔ)言賦值問(wèn)題？已經(jīng)有3人回復(fù)
【求助】Arcgis與GS+在地統(tǒng)計(jì)分析中的區(qū)別已經(jīng)有39人回復(fù)

1樓 2011-08-05 16:56:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲(chóng) (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時(shí)
蟲(chóng)號(hào): 183750
注冊(cè): 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

ian_zhangty(金幣+4): 2011-08-06 10:36:46

那就上我說(shuō)的那兩個(gè)上面找找，也可以去看看 R (http://www.r-project.org) 或者 dataplot (http://www.itl.nist.gov/div898/software/dataplot/) 的源碼中關(guān)于 p 的計(jì)算。后者主要是 fortran 編的，前者有一部分是 C，有一部分是 fortran，還有大部分是 R....

贊一下

回復(fù)此樓

4樓2011-08-06 09:32:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 5 個(gè)回答

snoopyzhao

至尊木蟲(chóng) (職業(yè)作家)

程序強(qiáng)帖: 16
應(yīng)助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時(shí)
蟲(chóng)號(hào): 183750
注冊(cè): 2006-02-13
專業(yè): 污染生態(tài)化學(xué)

【答案】應(yīng)助回帖

★
jjdg(金幣+1): 感謝您的關(guān)注 2011-08-06 01:01:29
ian_zhangty(金幣+1): 2011-08-06 09:10:51

計(jì)算 p 值可不是簡(jiǎn)單的事情，呵呵……
去 statlib 和 netlib 上找找

其實(shí)，這種事情，找個(gè)統(tǒng)計(jì)軟件做了就結(jié)了……

贊一下(1人)

回復(fù)此樓

2樓2011-08-05 23:00:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ian_zhangty

木蟲(chóng) (著名寫(xiě)手)

程序強(qiáng)帖: 1
應(yīng)助: 1 (幼兒園)
金幣: 4358.4
散金: 63
紅花: 3
帖子: 1043
在線: 86.6小時(shí)
蟲(chóng)號(hào): 296621
注冊(cè): 2006-11-12
專業(yè): 植物學(xué)研究的新技術(shù)、新方

★
dubo(金幣+1): 歡迎常來(lái)程序語(yǔ)言版討論 2011-08-21 14:00:44

引用回帖:

2樓: Originally posted by snoopyzhao at 2011-08-05 23:00:50:
計(jì)算 p 值可不是簡(jiǎn)單的事情，呵呵……
去 statlib 和 netlib 上找找

其實(shí)，這種事情，找個(gè)統(tǒng)計(jì)軟件做了就結(jié)了……

我也知道的。但是由于我的整個(gè)程序都是用fortran編的，如果有直接的程序就非常省事了

贊一下(1人)

回復(fù)此樓

3樓2011-08-06 09:10:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ian_zhangty

木蟲(chóng) (著名寫(xiě)手)

程序強(qiáng)帖: 1
應(yīng)助: 1 (幼兒園)
金幣: 4358.4
散金: 63
紅花: 3
帖子: 1043
在線: 86.6小時(shí)
蟲(chóng)號(hào): 296621
注冊(cè): 2006-11-12
專業(yè): 植物學(xué)研究的新技術(shù)、新方

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★
jjdg(金幣+1): 感謝分享你的解決方案 2011-08-06 15:28:45
dubo(金幣+5, 程序強(qiáng)帖+1): good 2011-08-21 14:01:09

已經(jīng)解決，方案如下，謝謝

FUNCTION betai(a,b,x)
REAL betai,a,b,x
!USES betacf,gammln
REAL bt,betacf,gammln
if(x<0..or.x>1.) pause 'bad argument x in betai'
if(x==0..or.x==1.) then
  bt=0.
else
  bt=exp(gammln(a+b)-gammln(a)-gammln(b)+a*log(x)&
   +b*log(1.-x))
endif
if(x<(a+1.)/(a+b+2.)) then
  betai=bt*betacf(a,b,x)/a
  return
else
  betai=1.-bt*betacf(b,a,1.-x)/b
  return
endif
END FUNCTION betai

FUNCTION betacf(a,b,x)
INTEGER maxit
REAL betacf,a,b,x,EPS,fpmin
PARAMETER (maxit=100,EPS=3.e-7,fpmin=1.e-30)
INTEGER m,m2
REAL aa,c,d,del,h,qab,qam,qap
qab=a+b
qap=a+1.
qam=a-1.
c=1.
d=1.-qab*x/qap
if(abs(d) d=1./d
h=d
do m=1,maxit
  m2=2*m
  aa=m*(b-m)*x/((qam+m2)*(a+m2))
  d=1.+aa*d
  if(abs(d)   c=1.+aa/c
  if(abs(c)   d=1./d
  h=h*d*c
  aa=-(a+m)*(qab+m)*x/((a+m2)*(qap+m2))
  d=1.+aa*d
  if(abs(d)   c=1.+aa/c
  if(abs(c)   d=1./d
  del=d*c
  h=h*del
  if(abs(del-1.) betacf=h
return
  end if
end do
pause 'a or b too big, or maxit too small in betacf'
END FUNCTION betacf

FUNCTION gammln(xx)
REAL gammln,xx
INTEGER j
DOUBLE PRECISION ser,stp,tmp,x,y,cof(6)
SAVE cof,stp
DATA cof,stp/76.18009172947146d0,-86.50532032941677d0,&
   24.01409824083091d0,-1.231739572450155d0,&
   .1208650973866179d-2,-.5395239384953d-5,&
   2.5066282746310005d0/
x=xx
y=x
tmp=x+5.5d0
tmp=(x+0.5d0)*log(tmp)-tmp
ser=1.000000000190015d0
do j=1,6
  y=y+1.d0
  ser=ser+cof(j)/y
end do
gammln=tmp+log(stp*ser/x)
END FUNCTION gammln

program pvalue
REAL df, r, t, prob
df=17-2 ! 17 is sample size

r=0.188446645 ! pearson coefficient
t=(abs(r)*sqrt(df))/sqrt(1-r**2)
prob=betai(0.5*df,0.5,df/(df+t**2))
write (*,*) prob

end program pvalue

贊一下(2人)

回復(fù)此樓

5樓2011-08-06 10:41:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 5 個(gè)回答

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 268求調(diào)劑 +3	簡(jiǎn)單點(diǎn)0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 08工科 320總分求調(diào)劑 +4	梨花珞晚風(fēng) 2026-03-17	4/200	2026-03-17 13:38 by houyaoxu
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7	Ncdx123456 2026-03-13	8/400	2026-03-16 12:15 by karry wen
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +5	Becho359 2026-03-11	5/250	2026-03-14 11:35 by 哦哦123
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[基金申請(qǐng)] 有必要更換申報(bào)口嗎 20+3	fannyamoy 2026-03-11	3/150	2026-03-14 00:52 by zhanghaozhu
[考研] 0703，333分求調(diào)劑一志愿鄭州大學(xué)-物理化學(xué) +3	李魔女斗篷 2026-03-11	3/150	2026-03-13 22:24 by JourneyLucky
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 材料工程調(diào)劑 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant