mixed-models

    0熱度

    1回答

    我想識別lmer模型中的異常值(lme4包)。我是而不是有興趣去除它們(什麼LMERConvenienceFunctions包) - 我只是想看到列出的異常值。 一個模式,我用的例子: model1<-lmer(Value~ Moisture + Planting + (day|plot), data=plants1)

    0熱度

    1回答

    調整意味着我運行模式是這樣的: model<-lme(Juice~1+Fruit ,random=list(Place=pdIdent(~1)), method="REML", control=lmeControl(niterEM=150 ,msMaxIter=200), na.action=na.omit, data=customData, keep.

    0熱度

    1回答

    我打算使用像varIndent(from = ~1 |Sex)這樣的結構權重和一個在R中的lme函數中稱爲w的定製矢量權重。所以我想知道是否有人知道如何組合函數中的權重。這是一個例子(只爲這個問題作出) n=25 data=data.frame( y = rnorm(n), Sex = as.factor(sample(c('Male','Female','Female'),s

    3熱度

    1回答

    在Julia中可能存在LMM中隨機效應的相互作用嗎? 這給出了一個錯誤 model = fit!(lmm(@formula(response ~ 1 + A*B + (1+A*B|sub)), data) ERROR: MethodError: no method matching getindex(::DataFrames.DataFrame, ::Expr) 拆開方面也沒有幫助。 mod

    1熱度

    1回答

    我想循環混合模型和Tukey測試。 我想要做的就是重複擬合和3列(每個包含一個響應)和4個亞組(總數= 12)的比較。 類似的數據框可在此處獲得:https://drive.google.com/open?id=0Bwrsa11LAnrgTXMzWk1fYXR1MHM。這3個部分是「RESP_1」,「RESP_2」和「RESP_3」列,這些子組是「層」列的變量。 我獲得由模型和一個單個響應的adj

    0熱度

    1回答

    一個混合效應式I具有式: my_formula <- a ~ b + (1|d) + e 欲編程方式從該公式中刪除e變量。 我爲固定效應僅公式進行的是: drop_es <- function(f) { e_idxs <- grep(pattern = 'e', x = all.vars(f[[3]])) stats::formula(stats::dr

    1熱度

    2回答

    我想進行重複測量分析/縱向數據到以下問題: 「有在的區域被分析的16個樹和在B地區有16個。在每個地區,冬季分析8棵樹,夏季分析8棵樹,但它們不是同一棵樹。考慮到在對每棵樹直徑五個不同深度的starch's知覺。」 tree Region season depth starch 1 A W 1 0.07 1 A W 2 0.10 1 A W 3 0.13 1 A

    2熱度

    1回答

    我有60棵樹的數據框。根據DOY,我需要爲EWMZ找到一個斜坡。我使用了混合模型,發現了一個負斜率,但是在繪製曲線時,它給了我一個積極的迴歸曲線。我不知道我用ggplot進行繪圖的方式是否正確。 DOY EWMZ TREE 247 13,01 1 262 11,01 1 274 23,07 1 288 23,09 1 310 20,77 1 247 28,47 2 262 22,55

    0熱度

    1回答

    因此,當我認爲存在與特定特定組之間的依賴關係時,將使用混合效果迴歸模型。我附加了維基鏈接,因爲它的解釋比我好。 (https://en.wikipedia.org/wiki/Mixed_model) 儘管我認爲在很多場合我們需要考慮混合效應,但沒有太多模塊支持這一點。 R有lme4和Python似乎有一個類似的模塊,但它們都是統計驅動的;他們不使用梯度增強等成本函數算法。 在機器學習設置中,您將如

    0熱度

    2回答

    我有一個鳥類數量的數據框。我有參加者的身份證號碼,他們計算的鳥的數量,他們計算的年份,他們的緯度和長度座標以及他們的努力。我製作了這個模型: model = lmer(count~year+lat+long+effort+(1|participant), data = df) 我現在希望模型繪製來自同一數據集的預測值。所以,這些數據是1997 - 2017年的,我希望模型能夠給出我每年的預測值