Ограниченная взвешенная линейная регрессия в R - PullRequest
1 голос
/ 28 мая 2019

Я пытаюсь настроить ограниченную взвешенную линейную регрессию.То есть у меня есть набор данных i наблюдений и три разных значения x.Каждое наблюдение имеет вес.Я хочу выполнить взвешенную множественную линейную регрессию, используя ограничения, что weighted mean каждого значения x должно быть равно нулю, а weighted standard deviation должно быть равно единице.

Поскольку я новичок и не имею репутации, я не могу публиковать изображения с латексными формулами.Поэтому я должен записать их таким образом.

Первое ограничение $\sum_{i} w_{i} X_{i,k} = 0$ for k = 1,2,3.

Второй: $\sum_{i} w_{i} X_{i,k}^2 = 1$ for k = 1,2,3.

Это примерный набор данных:

y <- rnorm(10)
w <- rep(0.1, 10)
x1 <- rnorm(10)
x2 <- rnorm(10)
x3 <- rnorm(10)
data <- cbind(y, x1, x2, x3, w)
lm(y ~ x1 + x2 + x3, data = data, weigths = data$w)

Веса не обязательно должны быть равными для каждого наблюдения, но должны быть добавлены к одному.Я хотел бы включить эти ограничения в регрессию.Есть ли способ сделать это?

1 Ответ

0 голосов
/ 01 июня 2019

Возможно, вы могли бы использовать обобщенную линейную модель:

glm(y ~ x1 + x2 + x3, weights = w, data=data)

Данные должны быть data.frame (...).

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...