ggplot2 - zoomando su geom_smooth automaticamente (usando coord_cartesian)

Question

geom_smooth è grande, in gran parte perché fa una grande variazione. Tuttavia, a causa di ciò, è difficile vedere come varia sull'asse x quando viene rimpicciolito. Sto producendo circa 1000 grafici in cui ho bisogno di avere lo zoom ggplot2 via coord_cartesian. Tuttavia, ogni grafico avrebbe diversi limiti di zoom. C'è un modo che posso chiedere a ggplot2 di ingrandire per adattarsi al liscio? Sono interessato a soluzioni sia per lo zoom avanti sulla linea geom_smooth sia per la linea geom_smooth più l'area ombreggiata SE.

Per esempio, sarei interessato a sapere come avrei potuto trasformare questo:

ggplot(data=mtcars, aes(y=qsec,x=wt)) + geom_point() + geom_smooth() 

in qualcosa di simile a questo:

ggplot(data=mtcars, aes(y=qsec,x=wt)) + geom_point() + geom_smooth() + coord_cartesian(ylim = c(15,20)) 

senza specificare esplicitamente i limiti.

Answer 1

L'adattamento manuale dei modelli di livellamento offre molta più flessibilità per raggiungere questo e altri tipi di personalizzazione. Per la maggior parte dei progetti, inizio usando l'interfaccia in-line, ma di solito finisco per dover passare al calcolo manuale quando ho bisogno di altre modifiche.

Vedere anche § 9.1.1 nel libro di Hadley.

require(ggplot2) 

# Fit smooth manually 
fit = loess(qsec ~ wt, data=mtcars) 
newx = data.frame(wt=with(mtcars, seq(min(wt), max(wt), len=100))) 
pred = predict(fit, newdata=newx, se=T) 
pred = cbind(wt=newx, qsec=pred$fit, se=pred$se.fit) 

# Calculate limits based on extent of smooth geom 
ylims = with(pred, c(floor(min(qsec-se)), ceiling(max(qsec+se)))) 

# Plot 
dev.new(width=5, height=4) 
ggplot(data=mtcars, aes(y=qsec, x=wt)) + 
    geom_point() + 
    geom_smooth(aes(ymax=qsec+se, ymin=qsec-se), data=pred, stat='identity') + 
    coord_cartesian(ylim = ylims) 

Tuttavia, questo ancora non funziona per sfaccettature perché è possibile specificare solo, per esempio, scales='free', e non gli attuali limiti direttamente.