我的互動術語並未涵蓋所有群體
[英]My interaction term doesn't cover all groups
我正在嘗試以 3 個因素為組來執行線性混合模型。 但是運行語法，我只得到 2 個術語：time:grup_int2 或 time:grup_int3，如圖所示。 time:grup_int1應該有另一個術語。 我的參考因素組是1級，但是我不知道如何在循環中重新調整，我打算用package提供的數據 ...
glmmLasso 可以與 Tweedie 發行版一起使用嗎？
[英]Can glmmLasso be used with the Tweedie distribution?
我有一個線性混合效果 model，我正在嘗試進行變量選擇。 model正在測試1000個采樣點的森林退化程度。 大多數點沒有退化，因此因變量高度偏斜且有許多零。 因此，我使用 Tweedie 分布來擬合 model。我的主要問題是：Tweedie 分布實際上可以用在glmmLasso funct ...
如何報告混合效應中的隨機效應 model
[英]How to report randomeffect in the mixed effects model
我使用 R 中的iris數據集作為下面的示例。 隨機截距和斜率如下所示。 這意味着 setosa 的截距為 -0.49549054（固定 + 隨機截距），斜率為 0.78331501（固定 + 隨機斜率），對嗎？ 因此，存在三對截距和斜率。 在一般線性 model 中，我們可以說 y = 截距 + ...
可視化混合效果 model
[英]Visualising a mixed effects model
因此，正如標題所暗示的那樣，我想形象化我擁有的這兩個模型，以便更好地理解它們之間的關系。 第一個是固定斜率 model，第二個是隨機斜率 model。我想 plot 它們分別，但我不知道該怎么做。 任何幫助將不勝感激。 數據示例： 我試過給他們 plot 但似乎無法讓他們工作 ...
如何根據日期在 R 中創建天數變量？
[英]How can I create a day number variable in R based on dates?
我想用參與者參加調查的日期（第一天、第二天、第三天等）創建一個變量。問題是有些參與者是在午夜后參加調查的。 例如，這是它的樣子： ID 日期 1個 08/03/2020 08:17 1個 08/03/2020 12:01 1個 08/04/2020 15:08 1個 08/04/2020 ...
如何正確解釋具有大 z 統計信息/沖突錯誤消息的 glmmTMB 模型？
[英]How to correctly interpret glmmTMB models with large z statistics/conflicting error messages?
我正在使用glmmTMB運行一個零膨脹的雙分量障礙 model 以確定某些協變量可能如何影響（1）魚的胃中是否有食物以及（2）如果胃中有食物，協變量會影響數量在它的胃里發現的獵物。 我的數據包括一條魚被捕獲的年份、捕獲的季節、性別、狀況、原產地、總海齡（1SW = 海上一年，MSW = 海上多年） ...
R 中的 GLMM 與 SPSS（收斂性和奇異性問題消失）
[英]GLMM in R versus SPSS (convergence and singularity problems vanish)
不幸的是，我在 R 中計算我的 GLMM 分析模型時遇到了收斂（和奇點）問題。當我在 SPSS 中嘗試時，我沒有收到這樣的警告消息，結果只是略有不同。 這是否意味着我可以毫無顧慮地解讀 SPSS 的結果？ 或者我是否必須測試奇點/收斂問題才能確定？ ...
如何在 R 中運行具有混合效應的多項邏輯回歸？
[英]How do I run a multinomial logistic regression with mixed effects in R?
我正在嘗試運行具有混合效果的多項式邏輯回歸。 假設我有以下變量： 參與者（十個參與者，每個參與者有 10 個觀察）單詞（十個不同的單詞，每個參與者看到每個單詞一次）IV（一些兩級分組變量）DV（可以是 1、2 或 3） 我將如何運行以 ppt 和 word 作為隨機變量的多項式邏輯回歸？ 這是一些 ...
如何將隨機效應納入非線性混合效應（nlme）model？
[英]How to incorporate a random effect into a nonlinear mixed effect (nlme) model?
我想建立一個非線性混合效應 model 來描述兩個變量“x”和“y”之間的關系，這兩個變量隨第三個變量“r”隨機變化，使用指數上升到最大值，如等式所述： y = theta(1-exp(-beta*x))。 我已經能夠使用 nls() 為 x 和 y 創建非線性 model，但是我沒有成功地將隨機 ...
如何在 lmer 中包含 3 個與單個預測變量的交互
[英]How to include 3 interactions with a SINGLE predictor in lmer
我知道我可以使用lm(a ~ (b + c + d)^2)來獲得模型中所有可能的雙向交互，但我只需要與單個預測變量的交互。 假設我想要B + C + D與預測變量E的可能交互。 我試過了： 我知道我可以用*或:對每個交互進行硬編碼，但我想有一種簡單的方法可以同時完成所有操作，不是嗎？ 提前致謝。 ...
從循環函數內的 lmerTest 模型中提取 p 值
[英]Extracting the p-value from an lmerTest model inside a loop function
我有許多線性混合模型，我已經安裝了lmerTest庫，因此函數的summary()會為我提供固定效果的 p 值。 我寫了一個循環函數，提取性別的固定效應：時間和性別：時間：感興趣的解釋變量。 現在還嘗試提取性別的 p 值：時間固定效應（第 1 步）以及性別：時間：解釋變量（第 2 步）。 通常我可以 ...
lme4 方差分析結果 + 效應量的導出表
[英]Export table of lme4 anova results + effect sizes
我正在嘗試導出一個表，其中包含我使用lme4運行的混合方差分析的結果。 這是模型。 我想導出方差分析統計數據加上各自的 eta 平方效應大小。 格式可以是 html、pdf 甚至 png。 原來是一場噩夢。 嘗試了幾個庫 xtable、stargazer、apa、knitr、apatable。 沒有 ...
如何運行混合指數衰減 model？
[英]How to run an exponential decay mixed model?
我不熟悉非線性回歸，希望對在 R 中運行指數衰減 model 提供一些幫助。請查看圖表了解數據的外觀。 我的直覺是指數 model 可能是一個不錯的選擇。 我有一種固定效應和一種隨機效應。 y ~ x + (1|隨機因子)。 如何獲取 R 中指數 model 的起始值（請假設我對非線性回歸一無所 ...
R_lme4_mixed-effects_modelling_SPSSexample_replication_cross-platform
[英]R_lme4_mixed-effects_modelling_SPSSexample_replication_cross-platform
你好 stackoverflow 社區！ 我是混合效果建模 (MEM) 或mixed-models的新手。 為了更好地理解 MEM，我決定復制Howard J. Seltman 博士的教科書《實驗設計與分析》中r （ lme4包）中的兩個例子。 教材中作者用SPSS解決了兩個例子，並收錄了相關的o ...
nlme::lme() 中隨機術語語法的解釋
[英]Explanation of random term syntax in nlme::lme()
我正在嘗試確定治療是否對兩個timepoints （兩個時間點 x 三個treatment ）之間的微生物組diversity有顯着影響。 有人可以解釋一下如何使用 R 中的nlme庫使用線性混合模型R嗎？ 特別是如何隨着時間的推移處理同一主題的重復采樣。 我已經看到使用了以下三種語法，但並不真正理 ...
glmmTMB model 的似然比？
[英]Likelihood ratio for glmmTMB model?
我有一個混合的 model，我試圖在其中找到隨機效應的重要性。 model 是我使用 R package glmmTMB 構建的具有零膨脹 beta 分布的混合 model，具有以下 function： 我試圖找出變量“Xnumber”的意義。 我讀過我需要做的是似然比檢驗，但不知道如何使用 g ...
從 R 調用 julia GeneralizedMixedModel
[英]Calling julia GeneralizedMixedModel from R
我正在嘗試使用 R 中的 julia 來擬合龐大的數據集。我最初嘗試使用 jglmm ( https://github.com/mikabr/jglmm )，但生成的對象很難使用。 我最初使用 JellyMe4 ( https://github.com/palday/JellyMe4.jl ) 時 ...
難以理解“固定效應”和“混合效應模型中的相關性”
[英]Trouble understanding "Fixed effects" and "Correlation in mixed effect models"
我在混合效應模型中有一個任務，但幾乎沒有關於“lmer”output 的信息。任務中使用的數據涉及植物高度的兩種處理類型（處理 OTC 和控制 CTL）。 這是使用的代碼： `mod_vivipara <- lmer(值 ~ 處理 + (1|站點), 數據 = bistorta_vivip ...
在 SPSS 中使用 GLM 重復測量構建自定義 Model
[英]Building Custom Model using GLM Repeated Measures in SPSS
我正在尋找使用 SPSS 版本 28 的支持，特別是關於如何使用 GLM 重復測量構建自定義 model。 我正在嘗試使用一個主體內變量和兩個主體間變量對協方差進行混合分析。 我在 IBM 網站上找到了關於我的主題的 SPSS 文檔： https://www.ibm.com/docs/en/sps ...
准泊松混合效應 model 對來自 R 中多個估算數據集的過度分散計數數據
[英]Quasi-Poisson mixed-effect model on overdispersed count data from multiple imputed datasets in R
我正在處理我可以單獨解決的三個部分的問題，但現在我需要一起解決它們： 極度偏斜、過度分散的相關計數變量（做某事時發生的事件數量）， 必須包括隨機效應， 許多缺失值 -> 多重插補 -> 10 個插補數據集。 為了解決前兩部分，我選擇了一個准泊松混合效應 model。由於stats::gl ...