[問題] 使用R plot likelihood function

作者: AmigoSafin   2018-11-27 11:38:34
[問題類型]:
程式諮詢(我想用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)是一條直線
有如我的心電圖(眼神死)
還請各位高手幫忙給小魯一點方向
謝謝~~
作者: HumuHumu (呼姆呼姆)   2018-11-27 15:59:00
作者: ksherry (K雪梨)   2018-11-29 11:39:00
fn 裡面 input 沒有 n
作者: AmigoSafin   2018-11-29 21:08:00
我後來做的跟一樓差不多https://reurl.cc/AKn93 像這樣 我沒有用sort我晚點看看一樓大大的跑出來圖是不是一樣 還有我不懂二樓大
作者: andrew43 (討厭有好心推文後刪文者)   2018-11-30 00:00:00
沒定義x啊原PO。
作者: AmigoSafin   2018-11-30 00:11:00
x <-rlaplace(n,scale=theta) 這個是嗎?
作者: andrew43 (討厭有好心推文後刪文者)   2018-11-30 01:46:00
好吧。儘量寫在一起寫完整吧。我和ksherry看到的毛病應該是一樣的,養成習慣把functio自定義的function需用到的參數都寫出來,包括theta,x,n而不是讓R自動到外部去幫你找內部未定義的變數在哪裡這不是錯誤,但養成好習慣以後比較不會出錯,也方便延用
作者: AmigoSafin   2018-11-30 10:30:00
我知道了!!謝謝A大~~我完食這題了 謝謝大家!已修正 抱歉太粗心了

Links booklink

Contact Us: admin [ a t ] ucptt.com