Параметричне, напівпараметричне та непараметричне завантаження для змішаних моделей

Наступні трансплантати взяті з цієї статті . Я новачок у завантажувальній програмі та намагаюся реалізувати параметричне, напівпараметричне та непараметричне завантажувальне завантаження для лінійної змішаної моделі з R bootпакетом.

R код

Ось мій Rкод:

library(SASmixed)
library(lme4)
library(boot)

fm1Cult <- lmer(drywt ~ Inoc + Cult + (1|Block) + (1|Cult), data=Cultivation)
fixef(fm1Cult)


boot.fn <- function(data, indices){
 data <- data[indices, ]
 mod <- lmer(drywt ~ Inoc + Cult + (1|Block) + (1|Cult), data=data)
 fixef(mod)
 }

set.seed(12345)
Out <- boot(data=Cultivation, statistic=boot.fn, R=99)
Out

Запитання

1. Як зробити параметричну, напівпараметричну та непараметричну завантажувальну роботу для змішаних моделей з bootпакетом?
2. Я думаю, я роблю непараметричні завантажувальні роботи для змішаної моделі в своєму коді.

Я знайшов ці слайди , але не може отримати пакет R merBoot. Будь-яка ідея, де я можу взяти цей пакет. Будь-яка допомога буде високо оцінена. Заздалегідь дякую за допомогу та час.

Завантаження в змішаних лінійних моделях дуже схоже на завантаження в регресії, за винятком того, що у вас є складність у тому, що ефекти поділяються на фіксовані та випадкові. При регресії, щоб зробити параметричний завантажувальний пристрій, ви підходите параметричну модель до даних, обчислюєте залишки моделі, завантажуєте залишки, беруть залишки завантажувальної програми та додаєте їх до встановленої моделі, щоб отримати зразок завантажувальної інформації для даних, а потім підходити до моделі до даних завантаження для отримання оцінок зразкових параметрів завантажувальної програми. Ви повторюєте процедуру, повторно завантажуючи початкові залишки, а потім повторюючи інші кроки в процедурі, щоб отримати ще одну оцінку параметрів вибірки завантажувальної програми. Для непараметричної завантажувальної стрічки ви створюєте вектор значень відповіді та коваріату та завантажуєте вибір векторів для вибірки завантажувальної програми. Із зразка завантажувальної машини ви підходите до моделі, щоб отримати параметри, і ви повторюєте процес. Єдина відмінність параметричної та непараметричної завантажувальної програми полягає в тому, що ви завантажуєте залишки для параметричного завантажувального пристрою, тоді як непараметрична завантажувальна програма завантажує вектори. У випадку змішаної моделі ви також можете мати напівпараметричну завантажувальну систему, обробляючи одні ефекти параметрично, а інші непараметрично. Якщо ваш код - це вектори завантаження, ви робите непараметричну завантажувальну систему. У мене немає конкретного рішення для забезпечення цього в R, але якщо ви подивитеся на книгу Ефрона та Тібшірані або мою книгу з Робертом Лабуддом, ви побачите код R для аналогічних типів моделей лінійної змішаної моделі.

Ви можете ознайомитись з bootMerфункцією у розробній версії lme4,

install_github("lme4",user="lme4")
library(lme4)

що може робити модельне (напів) параметричне завантаження змішаних моделей ... Просто перевірте ?bootMer