[英]how to add an intercept for linear regression when using matrix as an input in GLM Julia
我正在嘗試在 Julia 的 GLM 中使用線性回歸,使用矩陣作為輸入而不是 DataFrame。 輸入是: 但是當我嘗試使用lm function 來擬合它時,我發現攔截不是默認的: 我查看了 GLM 的官方文檔,但他們只解釋了 DataFrames 作為輸入的用法。 在不改變輸入的情況下使用矩 ...
[英]how to add an intercept for linear regression when using matrix as an input in GLM Julia
我正在嘗試在 Julia 的 GLM 中使用線性回歸,使用矩陣作為輸入而不是 DataFrame。 輸入是: 但是當我嘗試使用lm function 來擬合它時,我發現攔截不是默認的: 我查看了 GLM 的官方文檔,但他們只解釋了 DataFrames 作為輸入的用法。 在不改變輸入的情況下使用矩 ...
[英]Apply R glm predict function on dataframe by group
目前,我正嘗試在 dataframe 上應用 GLM 預測。dataframe 非常大,因此我想按塊應用預測。 我找到了一個解決方案,但它很不方便。 我首先創建一個空的 dataframe,然后使用 rbind。 有沒有更有效的方法來做到這一點? ...
[英]Why will glm model converge but caret train version of the same model will not?
據我所知,我已經使用基本 glm function 和插入符序列 function 以相同的方式指定了這個簡單的 GLM。但是,插入符版本不會收斂。 我指定火車 model 的方式是否遺漏了什么? ...
[英]Fitted values and confidence intervals on grouped bar chart
我有一個數據集,其中包含 5 個處理級別和 4 個時間級別,以及這 20 種組合的發芽值。 我正在嘗試制作一個分組條形圖,它將顯示我制作的 glm 中的擬合值和置信區間。 然而,這只是給我一個沒有置信區間的分組條形圖,並且會出現關於沒有正確尺寸的置信區間的錯誤。 置信區間似乎被制成向量而不是矩陣。 ...
[英]Using an exponential distribution in a predictive glm model in R
我正在使用 glm model 來預測我的因變量。 為此,我需要選擇一個變量分布族。 不幸的是,指數分布不是參數“家族”的可用對象的一部分。 現在我不知道如何繼續我的研究。 這是我的 model。有人知道我能做什么嗎? ...
[英]Issues related to Gamma Regression
我想通過迭代重新加權方法(手動)導出伽馬回歸系數。 當我在沒有for{}循環的情況下運行此代碼時,它可以正常工作,但在循環中它會產生NaN 。 我的代碼是: 請幫忙。 我的第二個問題與glm()有關。 有什么方法可以描述glm()使用了多少次迭代? 問候。 更新借助於此,我更新了這段代碼,但它不起作 ...
[英]R GLM Predict Error - factor has new levels
我正在使用 glm() 進行基本的邏輯回歸我將數據分成訓練和測試,使用 glm 構建 model,然后嘗試使用測試數據運行 predict()。 這是代碼 我遇到 我知道這個錯誤消息意味着 A 列中有一個值未在 model 中計算。但我檢查了訓練和測試數據中 A 列的唯一值,兩者具有完全相同的值 兩 ...
[英]Using proportions as an offset in a binomial (glm) model
我正在嘗試測試一種處理方法對候鳥種群中幼鳥比例的影響。 每天對這些鳥進行計數並識別為幼鳥或成鳥,但僅每隔一天進行一次處理。 沒有治療的天數被用作對照。 問題在於,預計幼魚在種群中的比例不僅會受到處理的影響,還會受到遷徙物候的影響。 例如,在某一天可能有更多的青少年遷移到研究區域,因此,不僅是治療, ...
[英]How to create a cglm Mat4 from a float[16] array?
我知道 C++ 有 glm::make_mat4() 和 #include <glm/gtc/type_ptr.hpp>,但 CGLM 中似乎沒有可用的 make_mat4() 函數。 有人可以告訴我如何從 float[16] 數組創建 CGLM Mat4 嗎? 我嘗試了以下方法,但 ...
[英]Error in MEEM(object, conLin, control$niterEM) in lme function
我正在嘗試將 lme 函數應用於我的數據,但模型給出了以下消息: dput 數據,從https://pastebin.com/tv3NvChR復制(太大而無法包含在這里) 有什么問題? 我看到一個應用lme4::lmer函數的解決方案,但是還有另一個選項可以繼續使用lme函數嗎? ...
[英]get an angle between 2 points and rotate a point about an other point with this angle - C++
我基本上是在嘗試根據貝塞爾曲線的點創建二維線。 貝塞爾曲線的所有點都放置得很好,一切看起來都井井有條。 從這些點開始,我在 z 軸上創建了另外 2 個點,它們將成為直線的邊界: 然后我通過圍繞主要點旋轉它們來更改這些點的位置: pm 是圍繞固定點 pf 旋轉的移動點 angle 是貝塞爾曲線上 ...
[英]How can I apply interaction term with linear predictors in mgcv package
我有一個關於在帶有 2 個線性預測變量的mgcv package 中使用交互項的問題。 我編寫了代碼以適應x1和x2之間的交互, mgcv::gam(y ~ x1 + x2 + ti(x1, x2, k = 3), method = "REML", data = DATA) 但我真的不確定這是在mg ...
[英]Predicting from glm function gives an "Error in log(time)"
我正在嘗試使用泊松回歸來預測故障率。 但是,我不知何故得到了一個我不明白的錯誤。 這是我的代碼: 我得到的錯誤說明如下: 這是什么意思,我該如何解決? ...
[英]Text Argument in R Function With GLM: Object Not Found
當我嘗試向 function 提供權重時,包含在 package 中的 function 會拋出錯誤。需要像這樣指定 package 調用的部分: 在function內部,下面兩行代碼用來指定一個名為weight的數據框object: 但是,function 的分析部分嘗試使用 glm 使用提供的權 ...
[英]plot(glm.out) using the wrong type of residuals for scale-location plot?
我正在嘗試為 R 中的負二項式 model 生成比例位置 plot。理想情況下,此圖將 plot 是 y 軸上標准化Pearson殘差的平方根。 但是,就我而言,y 軸被標記為“標准化偏差殘差”。 我似乎無法以任何方式強制執行此操作。 這是令人擔憂的,因為我的朋友運行完全相同的代碼獲得了不同的 p ...
[英]Getting Warning: " 'newdata' had 5 rows but variables found have 750 rows" even though column names are the same
我有一個名為PimaDiabetes的數據集。 並且可以從這里提取數據集。 我從中得出邏輯 model: 但是,每當我嘗試針對名為 ToPredict 的新數據集運行它時: 懷孕葡萄糖血壓皮膚厚度胰島素體重指數糖尿病家系年齡 4個 136 70 0 0 20 31.2 22 1個 12 ...
[英]How Do I Use bestglm For Subset Selection In Logistic Regression?
我有一個 dataframe 如下: 我想在針對h運行時對變量a - g執行詳盡的最佳子集選擇。 我試過以下方法: 結果如下: 我收到了所有這些我不確定其含義的警告。 此外,我不確定 model 是否自動讀取h變量作為目標變量。 我已經檢查了bestglm CRAN 手冊,但我仍然不確定錯誤的原因 ...
[英]How can I get a lme4::glmer to run
我有這個 csv 文件名稱 ID 年游客反饋阿努比斯 1個 2014 真的 0:01:17 雅典娜 2個 2014 錯誤的 0:01:53 我正在嘗試使用 lme4 運行一般線性混合 model 但我正在努力讓它運行。 我的代碼是 我不斷得到 我在沒有家庭的情況下嘗試過,也嘗試過二項式。 如果我更改 ...
[英]How to change ci_method to "wald" in a logistic regression in R
我需要對包含大約 1,000,000 個數據點的數據集運行邏輯回歸。 我跑了一個 model 我想使用plot_model()可視化結果 function 這需要永遠運行(我已經等了大約一個小時,但還沒有完成)。 我收到一條消息Profiled confidence intervals may t ...
[英]Is there a way to plot hurdle model results pscl package or a way to plot the count - zero truncated negbin part of hurdle model in r?
有沒有辦法讓 plot 跨過 model 結果為 R? 我能夠 plot 障礙 model(下方)的零部分(帶 logit 鏈接的二項式),但我不知道如何 plot model 的計數部分(帶對數鏈接的截斷負二項式)。 我正在使用 pscl package 作為障礙 model。 示例數據(df ...