2013-11-28 3 views
1

Я использую qplot построить и установить наклон к моим данным:получение коэффициентов регрессии от qplot

qplot(log(X),log(Y),geom=c("point","smooth"),method="gam",formula=y~ns(x,2)) 

Он работает отлично. Но как мне получить коэффициенты наклонной кривой регрессии? Я знаю, что могу получить наклон явно используя, например, nls. Тем не менее, я хотел бы знать, какой уклон qplot подходит для моих данных. Я благодарен за любой совет.

+1

см http://stackoverflow.com/a/9790803/1317221 –

+3

[** Это **] (http://stackoverflow.com/ вопросы/8845279/предсказанные значения для логистики-регрессии-from-glm-and-stat-smooth-in-ggplot2-are/8846361 # 8846361) также могут быть релевантными, то есть «geom_smooth не должен заменять фактическое моделирование " – Henrik

ответ

1

Как об этом:

g<-qplot(log(X),log(Y),geom=c("point","smooth")) 
lookinside<-ggplot_build(g)$data[[2]] 
smooth<-data.frame(cbind(x=lookinside$x,y=lookinside$ymax-lookinside$ymin)) 
lm(smooth$x~smooth$y) 

    Call: 
    lm(formula = smooth$x ~ smooth$y) 

    Coefficients: 
    (Intercept)  smooth$y 
    1.25482  0.06275 


qplot(smooth$x,smooth$y)+geom_abline(intercept=1.25482,slope=0.06275,color="red") 

enter image description here

Смежные вопросы