Я пытаюсь выполнить bootstrap в R для неотрицательных коэффициентов наименьших квадратов.
Коэффициенты ядра соответствуют оценочным значениям Y для хроматографического чарда, которые даны синим цветом в на этом графике .
Теперь мне нужно вычислить доверительные интервалы для 9 самых больших коэффициентов, используя непараметрическую c самозагрузку, и построить их (таким образом, 9 графиков) , которые должны выглядеть следующим образом .
У меня есть этот код:
B = 1000 # we use 1000 bootstrap replicates
boot_coefs = matrix(0, nrow=B, ncol=ncol(data$X)) # empty matrix to store nnls coefficients calculated on bootstrapped datasets
for (i in 1:B) {...}
Я пытался использовать sample(data, replace=TRUE) для точек, но когда я использую оценочные коэффициенты в качестве данных, R возвращает матрицу, заполненную нулями. Так что моя начальная загрузка просто не сработает. Кажется, я не могу понять, что мне нужно рассчитывать доверительные интервалы по коэффициентам (а не по таким параметрам, как среднее).
Может ли кто-нибудь объяснить мне, что означают эти графики и доверительные интервалы и, возможно, даже как мне их рассчитать с помощью начальной загрузки?
Заранее спасибо!