作者查詢 / clickhere
作者 clickhere 在 PTT 全部看板的留言(推文), 共1142則
限定看板:全部
看板排序:
1F→:X=2 is 0.75 quantile if P(X=1) is 1/4.07/21 02:36
2F→:Otherwise, this is not a probability distribution.07/21 02:36
5F→:bootstrap要resampling 1000筆, 500是不對的.07/21 02:23
6F→:把你寫好的產生係數的函數丟給boot(),boot.ci()就行了.07/21 02:24
7F→:sample的問題, boot()會處理.07/21 02:25
8F→:用sample(..., replace=T)可以直接rexampling.07/21 02:26
1F→:1&2同時成立, a才是a-th quantile.07/21 02:29
2F→:因為P(x=a)可能不為0.07/21 02:31
3F→:2.應該是P(X≧F^-1(a) ) ≦ 1-a才對!!07/21 02:32
4F→:不加等號, 你的cdf可能達不到 a, 尤其當X是離散分佈時.07/21 02:34
2F→:這跟Gibb無關吧?!07/09 08:13
3F→:如果m是正整數,那這是cauchy的變型,在R用rcauchy就行了.07/09 08:21
4F→:上面是錯的, m=2才是cauchy07/09 08:23
5F→:sorry, m=1才是07/09 08:24
6F推:所以你這是t分佈的變型. m應該是a的函數吧?07/09 08:29
7F→:在R中,用 rt(), 再套用rejection method就行了.07/09 08:30
1F→:1.找統算的書,方法很多. inverse cdf, rejection, ...07/08 04:16
2F→:2. atan() or atan2(), 請看 help.07/08 04:16
3F→:3. 這跟 R 沒太大關係.07/08 04:17
7F→:統計計算. 找 Simulation by Sheldon Ross07/09 08:05
8F→:inv. cdf就是general了.其他方法都是這個的延伸07/09 08:06
9F→:af 是啥?07/09 08:07
3F推:gsub("A+", "X", String) 請參考 regular expression07/03 03:37
5F→:2F是Unix用的,Window至少還需要Rtools,minGW06/18 06:16
6F→:R CMD SHLIB 才能造出*.dll on Windows as *.so on Unix06/18 06:19
7F→:http://www.murdoch-sutherland.com/Rtools/ 全包了06/18 06:27
9F→:看你的R是幾版的,就抓相對的Rtools2x.exe06/19 04:48
1F推:Singular Value Decomposition.06/16 16:03
3F推:a<-ts(1:5);a;b<-lag(a,k=1);b;is.ts(a);is.ts(b)06/15 12:46
4F→:str(a);str(b);06/15 12:50
5F→:Use object, save life. View() is only worked on MS.06/15 12:54
9F推:load nlme&MASS. geepack亦可. 但均無multinominal.06/14 03:08
10F→:family請看help(family). Multinominal的link是啥?06/14 03:10
11F→:是multivariate logit?06/14 03:11
12F→:可能要merge成2個類別, 再針對某兩類跑兩次glmm.06/14 03:15