PTT
Submit
Submit
選擇語言
正體中文
简体中文
PTT
R_Language
[問題] 使用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
是這樣嗎
https://pastebin.com/eScbMQck
作者:
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大~~我完食這題了 謝謝大家!已修正 抱歉太粗心了
繼續閱讀
[問題] 有重複觀測值的回歸該如何寫
leo790124
[問題] 能否ㄧ次匯出多個檔案呢? (write.table)
ronald736
[問題] 讀取多年維度為矩陣的變數值(nc檔)
AndrewShi
[問題] stata 指令轉成R
warpdrive
[問題] 將資料時間按月份計算(迴圈)
AndrewShi
[問題] 選取資料特定區域作圖
AndrewShi
[問題] 原始檔案為數字,R讀入後變成時間格式,求轉換法
paranoia5201
Re: [問題] 輸出的結果,如何完整貼到文件呢?(求救)
TreeMan
[問題] 限制CPU使用上限?
Edster
[問題] 輸出的結果,如何完整貼到文件呢?(求救)
ronald736
Links
booklink
Contact Us: admin [ a t ] ucptt.com