作者查詢 / phil5566

總覽項目：發文 | 留言 | 暱稱

作者 phil5566 在 PTT 全部看板的留言(推文), 共200則

限定看板：全部

看板排序：

全部R_Language91Statistics86Tech_Job12home-sale3AllTogether2Math2Salary2ShuangHe1SongShan1<< 收起看板(9)

首頁

尾頁

Re: [請益] 擎天峰和都更案

[ SongShan ]19 留言, 推噓總分: +12

作者: shchho - 發表於 2022/03/11 22:24(2年前)

14^F推phil5566: 推推，已成為鄰居03/12 20:14

Re: [請益] 擎天峰和都更案

[ ShuangHe ]19 留言, 推噓總分: +12

作者: shchho - 發表於 2022/03/11 22:24(2年前)

14^F推phil5566: 推推，已成為鄰居03/12 20:14

[心得] 板橋畢業心得文

[ home-sale ]850 留言, 推噓總分: +366

作者: Darkside18 - 發表於 2022/03/09 18:14(2年前)

359^F推phil5566: 年收300買2500都擔心自己要吃土03/09 21:57

360^F→phil5566: 只能說有爸媽靠真好03/09 21:58

[討論] 外商學歷普遍不高？

[ Tech_Job ]146 留言, 推噓總分: +45

作者: AlainDelon - 發表於 2021/12/08 17:07(2年前)

123^F→phil5566: 114純血在IC外商算還可以12/09 08:49

[請益] 大台北有準時下班的公司/產業嗎？

[ Tech_Job ]77 留言, 推噓總分: +41

作者: xinye0622 - 發表於 2021/10/03 17:29(2年前)

62^F→phil5566: 推外商，薪水快3M，每天工時六小時10/04 08:21

[請益] offer請益及規劃

[ Tech_Job ]12 留言, 推噓總分: +7

作者: TRF777 - 發表於 2021/09/15 14:59(2年前)

6^F→phil5566: 新北有外商IC很缺AE呀09/15 15:42

[問題] 模擬資料算估計值的MSE問題?

[ Statistics ]3 留言, 推噓總分: 0

作者: phil5566 - 發表於 2017/06/18 21:36(7年前)

3^F→phil5566: 感謝回答,不知道還有沒有其它說法?06/19 17:16

[問題] 設限資料作圖問題?

[ R_Language ]7 留言, 推噓總分: 0

作者: phil5566 - 發表於 2017/05/31 15:24(7年前)

7^F→phil5566: 太酷了~大家都code我都試過了~感謝,又學到東西了,謝謝06/01 15:37

[問題] 信賴區間和覆蓋率的問題?

[ Statistics ]50 留言, 推噓總分: +1

作者: phil5566 - 發表於 2017/02/21 02:10(7年前)

3^F→phil5566: 首先先用模擬的方式產生一組我有興趣研究的樣本,用這組02/21 13:52

4^F→phil5566: 樣本求出θ^,因為θ^很複雜,所以沒辦法用手算求出Var,只02/21 13:54

5^F→phil5566: 好用bootstrap的方式求stand error代替,假設bootstrap02/21 13:57

6^F→phil5566: 1000次,再每bootstrap一次下也可求出θ^,假設叫θb^,共02/21 14:00

7^F→phil5566: 有1000個θb^,這1000個θb^的s.d.,即bootstrap下的02/21 14:02

8^F→phil5566: stand error,接這建立C.I.,我建立C.I.的方式有四種,02/21 14:12

9^F→phil5566: 以下用R的語法說明較快,先說明一些代號,th0=θ^,02/21 14:17

10^F→phil5566: th=那1000組的θb^,alpha<-c(0.025,0.975),02/21 14:20

11^F→phil5566: 1.normalC.I.:th0 + qnorm(alpha) * sd(th)02/21 14:23

12^F→phil5566: 2.basic C.I.:2*th0-quantile(th,rev(alpha),type=1)02/21 14:25

13^F→phil5566: 3.percentile C.I.:quantile(th, alpha, type=6)02/21 14:27

14^F→phil5566: 剛剛忘了說我bootstrap的方式是採抽後放回的抽樣方式02/21 14:30

15^F→phil5566: 第4種 BCa C.I.最麻煩,zalpha <- qnorm(alpha)02/21 14:48

16^F→phil5566: z0<- qnorm(sum(th0 < th) / length(th))02/21 14:49

17^F→phil5566: L<- mean(th.jack) - th.jack02/21 14:49

18^F→phil5566: a <- sum(L^3)/(6 * sum(L^2)^1.5)02/21 14:49

19^F→phil5566: adj.alpha <- pnorm(z0 + (z0+zalpha)/(1-a*(z0+zalpha)02/21 14:50

20^F→phil5566: ))02/21 14:50

21^F→phil5566: quantile(th, adj.alpha, type=6)--->BCa C.I.02/21 14:50

22^F→phil5566: 這裡面的th.jack就是把我產生的那組資料02/21 14:53

23^F→phil5566: 假設叫T(n筆data),逐一的剔掉一筆data後,在剩餘的data下02/21 14:58

24^F→phil5566: 所算出的估計量,假設叫θj^,所以從第一筆開始剔除到最後02/21 15:00

25^F→phil5566: 一筆,會得到n筆的θj^,這n筆θj^就叫作th.jack02/21 15:02

26^F→phil5566: 更正:是在剩餘的data下所算出的"估計值",不是"估計量"02/21 15:04

27^F→phil5566: ---------------------------------------------------02/21 15:05

28^F→phil5566: 以上是產生一組data下所算出的4種C.I.,那我重複產生10002/21 15:06

29^F→phil5566: 組data下可算出100組4種C.I.,在分別去計算4種C.I.有沒有02/21 15:08

30^F→phil5566: 包含我所設定的參數θ的比率,即為文中的覆蓋率02/21 15:10

31^F→phil5566: 我比較疑惑的是我的4種C.I.都是建立在95%信心水準下02/21 15:13

32^F→phil5566: 照理說做100組應該要有95組左右的C.I都會包到設定的參數02/21 15:15

33^F→phil5566: 怎麼會掉到80多組?所以我想是不是估計量是有偏的,所以02/21 15:16

34^F→phil5566: 覆蓋率才會有落差?02/21 15:17

35^F→phil5566: 不好意思~我說了一大堆~我也不知道有沒有講到重點,實際02/21 15:28

36^F→phil5566: 上我的θ是一個向量,不過就只是每個分量的元素都做和我02/21 15:30

37^F→phil5566: 上面所說的一樣的事而已02/21 15:31

38^F→phil5566: 附上原始的codehttp://pastebin.com/jgUBxtvj02/21 15:32

43^F→phil5566: 推文裡前3種C.I.只要用到th0和th,th就是您說的一組數字02/21 17:41

44^F→phil5566: x就是輸入產生的data,而stat是算th.jack會用到的,簡單來02/21 17:45

46^F→phil5566: x和stat是算BCa C.I.會用到的東西所以我要輸入02/21 17:46

49^F→phil5566: boot.CI回傳的會是4種C.I.的值,我把4種C.I.的算法都寫02/21 17:48

50^F→phil5566: 在一boot.CI這個函數了02/21 17:49

[問題] 實務上如何用模擬來證明不偏性和有效性?

[ Statistics ]34 留言, 推噓總分: +1

作者: phil5566 - 發表於 2017/02/18 02:30(7年前)

3^F→phil5566: 產生很多組樣本算出估計量平均,根據大數法則平均估計量02/18 11:44

4^F→phil5566: 機率收斂到估計量的期望值,假設叫a,看|a-α|是否為0?對02/18 11:48

5^F→phil5566: 吧?那有|a-α|有說要多小才是不偏嗎?另外好奇問一下,02/18 11:51

6^F→phil5566: bootstrap1000次裡也1000組的估計量,直接把這1000組估計02/18 11:54

7^F→phil5566: 量算平均,可以代表估計量的平均值嗎?bootstrap算出的02/18 11:56

8^F→phil5566: stand error去比大小,可以衡量f和g的相對有效性嗎?謝謝02/18 11:58

14^F→phil5566: 那麼比較f和g的不偏,重複產生n組的samples一定要共用嗎?02/19 14:55

15^F→phil5566: 還是各別去產生n組samples去算不偏也行?02/19 14:57

16^F→phil5566: 該不會是重複產生n組樣本給f和g共用分別去算這n組的平均02/19 15:03

17^F→phil5566: 數和標準差,同時比較不偏和有效,而不須要用到bootstrap?02/19 15:04

18^F→phil5566: 謝謝02/19 15:05

22^F→phil5566: 那麼f和g的樣本點個數相同只是函數不同驗證不偏這樣有不02/19 15:12

23^F→phil5566: 合理的地方嗎?02/19 15:12

26^F→phil5566: 會這樣做是因為f和g函數太複雜沒法手算E(X)和V(X),只好02/19 15:14

27^F→phil5566: 用模擬來驗證,即使沒法驗證不偏,驗證漸進不偏和有效也02/19 15:16

28^F→phil5566: ok?比較兩者的優劣02/19 15:16

34^F→phil5566: 我也是這麼認為~感謝大家的意見02/20 13:59

首頁

尾頁