作者查詢 / phil5566

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

作者 phil5566 在 PTT [ R_Language ] 看板的留言(推文), 共91則

限定看板：R_Language

看板排序：

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

首頁

尾頁

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

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

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

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

[問題] 多變數函數中如何固定某些變數做圖/畫線?

[ R_Language ]37 留言, 推噓總分: +5

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

7^F→phil5566: 謝謝E大回答,您說的把資料改成array,是否就是02/11 00:47

8^F→phil5566: array(T,c(20,9))?但我不太懂為何可以省去02/11 00:48

9^F→phil5566: "t[t[,4]==i,2]這段"?還有我試了一下image好像是把data02/11 00:50

10^F→phil5566: http://imgur.com/a/G7Q33變成這種效果,我想做的效果是02/11 00:55

11^F→phil5566: 類似這種https://goo.gl/kwjQnV可以看出函數發生極值點02/11 01:01

12^F→phil5566: 的立體圖形,不知是否有誤解您的意思,能否給個code範例?02/11 01:03

17^F→phil5566: 可是我想畫的是函數圖形,和data結構有什麼關係?不懂?02/11 13:09

19^F→phil5566: 看了範例persp比較符合我的需求,但它是定出x,y和f02/12 01:50

20^F→phil5566: 假設x1代表x,x2代表y,那他的f就相當於我的02/12 01:53

21^F→phil5566: surlikelihood(a,T),a<-c(x1,x2,5)02/12 01:55

22^F→phil5566: x1 <- seq(-10, 10, length= 30)02/12 02:21

23^F→phil5566: x2 <- x102/12 02:22

24^F→phil5566: a<-c(x1,x2,5)02/12 02:22

25^F→phil5566: 前面有說x1,x2範圍,但先假設沒範圍限制下來作圖02/12 02:24

26^F→phil5566: f <-function(x,y) surlikelihood(a,T)02/12 02:25

27^F→phil5566: wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww02/12 02:28

28^F→phil5566: 上面這一行我只是想表示原本surlikelihood代入data T和02/12 02:30

29^F→phil5566: 固定x3=5之下,只剩x(x1)和y(x2)兩個變數的函數02/12 02:32

30^F→phil5566: 但這樣寫應該是錯的,不知道這一塊該怎麼處理才好?02/12 02:33

31^F→phil5566: 之後再接z <- outer(x, y, f)02/12 02:38

32^F→phil5566: z[is.na(z)] <- 102/12 02:38

33^F→phil5566: persp(x, y, z, theta = 30, phi = 30, expand = 0.5, c02/12 02:39

34^F→phil5566: ol = "lightblue")02/12 02:39

35^F→phil5566: 應該可以畫出立體圖......只是02/12 02:40

36^F→phil5566: f <-function(x,y) surlikelihood(a,T)這行該怎麼改才對02/12 02:41

37^F→phil5566: ?02/12 02:41

Re: [問題] 如何更有效率的整理data的問題?

[ R_Language ]12 留言, 推噓總分: +5

作者: celestialgod - 發表於 2017/01/31 17:42(7年前)

1^F推phil5566: 謝謝C大回答~這個方法蠻厲害的,完全沒用到廻圈01/31 22:21

2^F→phil5566: 但還是有幾點疑問?您說的outputMat排序是指01/31 22:23

3^F→phil5566: order(outputDF$Var1, outputDF$Var2)吧?order(X,Y)排序01/31 22:24

4^F→phil5566: 原則是否先比X再比Y,比完後越小的放越前面?01/31 22:26

5^F推phil5566: 最後一行我還可以理解,只是tapply那一行就看不太懂?01/31 22:29

6^F→phil5566: google"tapply"的用法還是有看沒有懂,能否麻煩C大解說一01/31 22:30

7^F推phil5566: 下?最後就是假設產生sample的那個步驟如果只有產生10個01/31 22:34

8^F→phil5566: 那您做g2那個動作,取min還會取到0.5嗎?謝謝01/31 22:35

9^F→phil5566: 可能是我沒說清楚,我指的是不論今天sample產生出來是什01/31 22:38

10^F→phil5566: 都一定要有0.5,1.5,......,19.5這20個數的存在01/31 22:40

11^F推phil5566: 讚讚讚~02/01 01:13

12^F推phil5566: 馬上來試試02/01 18:34

Re: [問題] 如何更有效率的整理data的問題?

[ R_Language ]20 留言, 推噓總分: +2

作者: clansoda - 發表於 2017/01/31 08:53(7年前)

1^F推phil5566: 謝謝clansoda大的回答,但是我用你的code跑過一遍之後發01/31 16:02

2^F→phil5566: 現如果一開始產生data時沒有出現1.5,分類一樣不會出現01/31 16:05

3^F→phil5566: 1.5,即產生出來的矩陣列數是小於120的矩陣01/31 16:06

6^F→phil5566: 例如:http://pastebin.com/5yzV0TyT01/31 16:08

7^F→phil5566: 不過再請教一下,您的code裡有一行是01/31 16:15

8^F→phil5566: dat[, .(Counting = .N), .(x, y)][order(x, y)]01/31 16:16

9^F→phil5566: dat是一個矩陣,為何dat[,,]裡會有兩個","01/31 16:18

10^F→phil5566: ?".N"和".(x, y)"代表什麼意思?01/31 16:24

11^F→phil5566: `:=`和:=又是什麼意思?謝謝01/31 16:28

19^F推phil5566: 受教了~謝謝01/31 17:12

[問題] 重新打code卻反而變慢,不知哪裡有問題?

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

作者: phil5566 - 發表於 2017/01/24 03:07(7年前)

7^F→phil5566: 回C大這個跟之前你改的那一個是不同的01/26 11:48

8^F→phil5566: 重貼舊版 http://pastebin.com/T4c4hwnn01/26 11:59

9^F→phil5566: 有喔~之前C大改的我有採用快超多的,只是這是不同的方法01/26 12:01

10^F→phil5566: 小弟資質駑鈍又是新手,簡化到新版就是我的極限了01/26 12:02

11^F→phil5566: 還是請大大再一次賜教01/26 12:04

12^F→phil5566: 重貼新版 http://pastebin.com/xqHHBHsx01/26 12:18

13^F→phil5566: 我在新版的第631行和第689行的確是用matrix存變數的01/26 12:21

14^F→phil5566: 不知C大指的是哪一部分?謝謝01/26 12:21

18^F→phil5566: 請教W大...google了benchmark還是有看沒有懂,能否給個例01/30 19:30

19^F→phil5566: 子,謝謝01/30 19:30

[問題] 在R上使用R2OpenBUGS如何使程式跑得更快?

[ R_Language ]8 留言, 推噓總分: +1

作者: phil5566 - 發表於 2017/01/22 02:14(7年前)

5^F→phil5566: 感謝C大回答,有問題我會再發問~01/22 19:49

Re: [問題] boot函數和optim函數搭配使用的問題

[ R_Language ]6 留言, 推噓總分: +3

作者: f496328mm - 發表於 2017/01/10 16:07(7年前)

3^F推phil5566: 讚讚讚~我再來試試看01/11 13:07

4^F推phil5566: 謝謝f大提醒01/11 17:23

5^F推phil5566: 也謝謝celestialgod的幫忙,真的學到蠻多的01/11 17:25

6^F→phil5566: ^大大01/11 17:26

[問題] boot函數和optim函數搭配使用的問題

[ R_Language ]32 留言, 推噓總分: +2

作者: phil5566 - 發表於 2017/01/10 00:14(7年前)

5^F→phil5566: 那如果repeat loop強迫 optim 直到收斂解? 要怎麼做呢?01/10 01:32

7^F→phil5566: 會用boot是因為我不寫bootstrap抽樣的code,恰好最複雜01/10 01:34

11^F→phil5566: 的部分已經有專家寫好function了01/10 01:35

16^F→phil5566: @x88776544能否舉個例子一下,沒用過loop,不清楚如何使用01/10 13:11

17^F→phil5566: @f496328mm 70%的意思是什麼?是只有70%的結果接受,剩下01/10 13:16

18^F→phil5566: 30%直接捨棄嗎?謝謝01/10 13:17

32^F→phil5566: 感謝各位大大的回答,我花個時間吸收一下,有問題再發問01/11 13:05

[問題] table()函數問題

[ R_Language ]23 留言, 推噓總分: +3

作者: phil5566 - 發表於 2016/12/27 02:11(7年前)

6^F→phil5566: 回andrew43大,您的方法取出來後,無法和其他向量作運算12/27 14:49

7^F→phil5566: 回koai大,你的方法,我有成功!!!12/27 15:01

8^F→phil5566: 回f496328mm大,我您的方法,作完後好像還table的型式12/27 15:02

9^F→phil5566: 不是矩陣的型式,謝謝版上高手指教,我學很到很多,感謝12/27 15:04

10^F→phil5566: 附上小弟後來的作法以及跑出來的結果12/27 16:18

11^F→phil5566: f<-as.double(names(table(T[,6])))12/27 16:19

12^F→phil5566: e<-data.frame(cumsum(table(T[,6]))/sum(table(T[,6]))12/27 16:20

13^F→phil5566: D<-rbind(f,e[,1])12/27 16:20

14^F→phil5566: 處理的data,以及跑的過程 http://codepad.org/Ju9SFS2L12/27 16:22

15^F→phil5566: rbind的結果好像怪怪的,列的部分為何不是顯示"[1] [2]"12/27 16:24

16^F→phil5566: 而是顯示"f "12/27 16:25

17^F→phil5566: 有人有更好的作法嗎?謝謝12/27 16:26

23^F→phil5566: 讚讚讚~ok了12/27 18:07

[問題] 函數boot.ci的問題?

[ R_Language ]9 留言, 推噓總分: +2

作者: phil5566 - 發表於 2016/12/21 18:03(7年前)

3^F→phil5566: 謝謝a大,試過您提供的方法,有成功,又學到一招,感謝12/23 19:33

4^F→phil5566: 不知道是不是boot.ci這個函數一次只能算第一列的C.I,還12/23 19:36

5^F→phil5566: 是要作什麼調整之後,第二列2~第5列的C.I也可以計算12/23 19:38

9^F→phil5566: 了解了,感謝12/23 19:51

首頁

尾頁