вычисление y_pred в регрессии наименьших квадратов (R)?

У меня проблемы с вычислением y_pred в регрессии наименьших квадратов. Идея примерно такая:

mydata <- read.csv("G:\\sample.csv",header=T)
x<-rep(mydata$wavelength,each=119)
y<-c(mydata$v1,....mydata$v119)
lm(y~x)

Образцы данных можно загрузить по адресу: https://drive.google.com/file/d/0B86_a8ltyoL3Y3BhU2xFVVo5dnM/view?usp=sharing

В файле переменная «Длина волны» - это x, где для каждого x есть несколько значений y, измеренных в разное время, как указано переменными от V1 до V119.

Я не уверен, что регрессия y (множественный) ~ x (один) ... Может ли кто-нибудь помочь вычислить y_pred в этом случае?

Большое спасибо!


person Vicki1227    schedule 08.10.2014    source источник
comment
Либо используйте predict.lm (но внимательно прочтите страницу справки), либо возьмите значения наклона и перехвата из вашего результата и вручную вычислите новые пары x, y.   -  person Carl Witthoft    schedule 08.10.2014
comment
Начните с чтения страницы справки для lm с помощью команды help('lm').   -  person    schedule 08.10.2014


Ответы (1)


Я думаю, что ты хочешь просто

mydata <- read.csv("G:\\sample.csv",header=T)
lm(Wavelength ~ ., data = mydata)

Это делает регресс Wavelength по отношению ко всем другим столбцам в вашем фрейме данных.

В вашем звонке

x<-rep(mydata$wavelength,each=119)

x заканчивается нулем. Вы должны использовать Wavelength с заглавной буквы.

person Sean Hughes    schedule 08.10.2014
comment
Это не отвечает на вопрос, как рассчитать прогнозируемые значения на основе lm модели. - person ; 08.10.2014
comment
Я использовал прогноз (), чтобы получить прогнозируемое значение. - person Vicki1227; 09.10.2014