作者AmigoSafin ()
看板R_Language
標題[問題] 使用R plot likelihood function
時間Tue Nov 27 11:38:34 2018
[問題類型]:
程式諮詢(我想用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
推 ksherry: fn 裡面 input 沒有 n 11/29 11:39
→ AmigoSafin: 我後來做的跟一樓差不多 11/29 21:08
→ 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