Как создать косинорную модель с линейным компонентом в R? - PullRequest
Новая
       8

Как создать косинорную модель с линейным компонентом в R?

0 голосов
/ 07 мая 2018

Привет, сообщество StackOverflow, Я не программист, а пользователь-самоучка для статистики и визуализации данных. Это мой первый вопрос, так как я всегда нахожу сообщения от других участников, ответивших на мои вопросы, так что, спасибо и прошу прощения за любые ошибки этикета в моем вопросе.

Я пытаюсь создать косинор модель с линейным компонентом в R.

Так, что ? = ???? {2? ((? − ∅) / τ)} + ?? + ?. Где A - амплитуда синусоиды, τ - период, установленный на 24 часа, t - время, Ø - акрофаза синусоиды, а C и D - наклон и y-пересечение линейного компонента, соответственно. Я посмотрел документацию для cosinor.lm {cosinor} и cosinor {psych}, но не вижу этой возможности.

Может ли линейный компонент быть добавлен в косинор модель или есть другая модель, которая может выводить значимость и значения для амплитуды, акрофазы, наклона и y-пересечения?

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

library(cosinor)
library(ggplot2)

#data frame
time<- c(2, 6, 10, 18, 22, 26, 30, 34, 38)
resp<- c(2.54, 0.13, -0.38, -0.57,  0.11, -0.20, -0.26, -0.62, -0.73)
df<-data.frame(time,resp)

# data fits both linear and cosinor models (but not very well)
cos.model <- cosinor.lm(resp ~ time(time), period=24, data=df)
summary(cos.model)
l.model <- lm(resp ~ time, data=df)
summary(l.model)

#plot with a loess smoother looks like a combination of the models 
#plot to see fit

ggplot(df, aes(x = time, y = resp))+
  geom_point()+
  geom_smooth(method = "lm", se = FALSE)

linear plot

ggplot.cosinor.lm(cos.model)+
  geom_point(data=df, aes(y=resp, x=time))

cosinor plot 

ggplot(df, aes(x = time, y = resp))+
  geom_point()+
  geom_smooth(method ="loess", se = FALSE)

loess plot

Если бы кто-нибудь мог помочь мне создать модель, включающую косинор и линейные компоненты, я был бы очень признателен!

Заранее спасибо.

...