看板 R_Language 關於我們 聯絡資訊
[問題類型]: 程式諮詢(我想用R 做某件事情,但是我不知道要怎麼用R 寫出來) [軟體熟悉度]: 入門(寫過其他程式,只是對語法不熟悉) [問題敘述]: 請簡略描述你所要做的事情,或是這個程式的目的 各位大大好 想請教一下如何寫R code畫出laplace distribution的likelihood function? p.d.f: f_x(x)=(1/2\theta)exp{(-1/\theta)*|x|} 現在generate a sample with n=100, theta=2 我把function寫出來為: library(VGAM) theta <-2; n <- 100 x <-rlaplace(n,scale=theta) fn <- function(theta){(1/2*theta)^n*exp((-1/theta)*sum(abs(X)))} 這部分沒問題 但接下來該如何用這個function畫圖就困擾了我 plot(fn)是一條直線 有如我的心電圖(眼神死) 還請各位高手幫忙給小魯一點方向 謝謝~~ -- -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 8.41.66.212 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1543289917.A.5EB.html
HumuHumu: 是這樣嗎 https://pastebin.com/eScbMQck 11/27 15:59
ksherry: fn 裡面 input 沒有 n 11/29 11:39
AmigoSafin: 我後來做的跟一樓差不多 11/29 21:08
AmigoSafin: https://reurl.cc/AKn93 像這樣 我沒有用sort我晚點 11/29 21:09
AmigoSafin: 看看一樓大大的跑出來圖是不是一樣 還有我不懂二樓大 11/29 21:10
andrew43: 沒定義x啊原PO。 11/30 00:00
AmigoSafin: x <-rlaplace(n,scale=theta) 這個是嗎? 11/30 00:11
andrew43: 好吧。儘量寫在一起寫完整吧。 11/30 01:46
andrew43: 我和ksherry看到的毛病應該是一樣的,養成習慣把functio 11/30 02:10
andrew43: 自定義的function需用到的參數都寫出來,包括theta,x,n 11/30 02:11
andrew43: 而不是讓R自動到外部去幫你找內部未定義的變數在哪裡 11/30 02:13
andrew43: 這不是錯誤,但養成好習慣以後比較不會出錯,也方便延用 11/30 02:22
AmigoSafin: 我知道了!!謝謝A大~~我完食這題了 謝謝大家! 11/30 10:30
※ 編輯: AmigoSafin (129.21.70.200), 11/30/2018 10:32:02
AmigoSafin: 已修正 抱歉太粗心了 11/30 10:32