國立成功大學

#請益 #R語言

2020年6月26日 20:48
小妹已請教好久google大神了,但真的還沒找到跟我類似情況的人😭想請問各位懂R的版友,關於我跑程式以下的小疑問🙏 上圖為卜瓦松分配,套入9個自變數,一個應變數(共十個變數的資料)時所跑出來的,一切正常 Number of iterations in BFGS optimization: 31 Log-likelihood: -1610 on 20 Df Warning message: In value[[3L]](cond) : 系統計算上是獨特的: 互反條件數=3.31014e-54FALSE > summary(n1<- glm(J~A+B+C+D+E+F+G+H+I,data = `V二級發布區自變數+應變數for.R`)) 同一份資料,零膨脹模型的跑出來怎麼長這樣🥺?想請問底下Warning message的部分是什麼意思呢? 感恩!
7
回應 12
文章資訊
共 12 則留言
國立成功大學
換個package看看,我記得有專們在寫這個的書;之前也搞死我 不過個人的直覺是自變數太多導致,或是你研究設計或資料不符合零膨脹的數學假設
原 PO - 國立成功大學
B1 自變數太多有可能耶!因為五筆就跑得出來!然後我自變數F(也就是第六個自變數)的部分在輸入的時候他自己變成藍色的,輸入到E就去跑,資料就都跑得出來,不知道什麼意思😂😂 謝謝學長的指教!書的部分我來找找😆
國立成功大學
B2 書名:Mixed Effects Models and Extensions in Ecology with R 相關章節:Chapter 11 Zero-Truncated and Zero-Inflated Models for Count Data
國立成功大學
B2 這部分我就不清楚了耶⋯拍謝 考慮把f拿掉看看嗎XD 不過我敢肯定的是⋯用常用套件跑出一堆na,不是假設錯誤或是資料不乾淨等
國立嘉義大學
如果你是要用特定的package跑 可以多嘗試幾個組合看看哪裡會有問題 如果你確定自己應該對 可以寫信去問package的作者是不是有bug... 會這樣說是因為本身碰過類似問題 某package在餵入特定變數數量時會運行發生問題 這個問題後來也被作者證實有 到這個方向的話除非你有能力去他github之類的抓原始碼下來改或等他改 不然另外找package可能比較快
原 PO - 國立成功大學
B3 哇賽!!謝謝學長嗚嗚嗚 人真的太好了,章節都出來了🥺❤️連假完我回學校去借這本書來讀!٩(˃̶͈̀௰˂̶͈́)و 不過跟你說的好像一樣,好像是資料不乾淨的問題,零不夠多或是資料裡有很多0但非零的資料數值很大! 我試過把F拿掉,也試過其他各種組合,唯獨A+B+C+D+E 可以全部跑出來,其他怎麼配都沒辦法跑出完整的東西,這部分我可能要再研究一下,謝謝你!
原 PO - 國立成功大學
B5 組合的部分蠻有趣的,我可能要研究是不是需要換個模型跑,或是看怎麼搭配那個組合了😭,寫信去給作者的部分感覺好厲害!!不過我的東西程度應該還不到那邊,應該是我的問題不是作者的問題XD(資料的問題),弱弱的請問R要怎麼新增package呢?另外找package是什麼意思呢?
國立嘉義大學
B7 新增不就是用install的語法就可以了?我不知道你的問題點 另外找的部分我以計算生物多樣性指數為例好了 線上計算生物多樣性指數的package可能有好多個 當你用的出問題換一個就有可能解決問題 只是要注意function背後參考的文獻公式是否相同或是有變形 模型我比較不熟,不過要換的話就大概要確認運作原理都一樣之類的吧 或是你能理解的
原 PO - 國立成功大學
B4 B5 謝謝你們!結果我發現他沒辦法跑資料超過四位數的,超過四位數我把他改單位,像是改成千個,就可以了,超奇怪的😂
國立成功大學
B9 硬體問題吧XD
原 PO - 國立成功大學
B10 蛤?哈哈哈 原來這個跟硬體有關係呀?是不是我該換電腦了 MacBook Air 用了五年呀 XD
國立成功大學
B11 那真的該換一台了!