作者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://webptt.com/m.aspx?n=bbs/R_Language/M.1543289917.A.5EB.html
2F：推 ksherry: fn 裡面 input 沒有 n 11/29 11:39
3F：→ AmigoSafin: 我後來做的跟一樓差不多 11/29 21:08
5F：→ AmigoSafin: 看看一樓大大的跑出來圖是不是一樣 還有我不懂二樓大 11/29 21:10
6F：→ andrew43: 沒定義x啊原PO。 11/30 00:00
7F：→ AmigoSafin: x <-rlaplace(n,scale=theta) 這個是嗎？ 11/30 00:11
8F：→ andrew43: 好吧。儘量寫在一起寫完整吧。 11/30 01:46
9F：→ andrew43: 我和ksherry看到的毛病應該是一樣的，養成習慣把functio 11/30 02:10
10F：→ andrew43: 自定義的function需用到的參數都寫出來，包括theta,x,n 11/30 02:11
11F：→ andrew43: 而不是讓R自動到外部去幫你找內部未定義的變數在哪裡 11/30 02:13
12F：→ andrew43: 這不是錯誤，但養成好習慣以後比較不會出錯，也方便延用 11/30 02:22
13F：→ AmigoSafin: 我知道了！！謝謝A大～～我完食這題了 謝謝大家！ 11/30 10:30
※ 編輯: AmigoSafin (129.21.70.200), 11/30/2018 10:32:02
14F：→ AmigoSafin: 已修正 抱歉太粗心了 11/30 10:32