小编Jav*_*ier的帖子

我正在尝试存储以下代码的结果，但是我只能想出一种解决方案，以最小的残差平方和保存模型的结果。这是有用的，直到结果在c和γ的范围内为止，因此我需要评估其他点的特征。为此，我需要存储每次迭代的结果。有谁知道在这种情况下该怎么做？

提前致谢！

dlpib1 <- info$dlpib1
scale <- sqrt(var(dlpib1))
RSS.m <- 10

for (c in seq(-0.03,0.05,0.001)){
  for (gamma in seq(1,100,0.2))
    {
    trans <- (1+exp(-(gamma/scale)*(dlpib1-c)))^-1
    grid.regre <-lm(dlpib ~ dlpib1 + dlpib8 + trans + trans*dlpib1 + 
                  + I(trans*dlpib4) ,data=info) 
coef <- grid.regre$coefficients
RSS <- sum(grid.regre$residuals^2)

if (RSS < RSS.m){
  RSS.m <- RSS
  gamma.m <- gamma
  c.m <- c
  coef.m <- coef
  }
 }
}
grid <- c(RSS=RSS.m,gamma=gamma.m,c=c.m,coef.m)
grid`

我有下图.正如你所看到的,它由两条实线组成geom_line,并用灰色线条构成geom_point.图表正是我需要的图形,但是,图例并未显示灰线也有点 - 需要包括它,因为图像将以黑白打印.

这是我正在使用的代码

library(ggplot2)
library(reshape2)

ggplot(base_fin, aes(x = fecha, y = valor)) +
     geom_line(aes(colour = Serie, linetype= Serie, size = Serie)) + 
     facet_wrap(~ indicador_f, scales = "free") + theme_bw() + 
     theme(legend.title = element_blank(), 
           legend.position="bottom", 
           text=element_text(family="Calibri"),
           plot.caption=element_text(size =8)) +
     labs( x =" ", 
           y = "Balance %",
           caption = "Fuente: Fedesarrollo - Encuesta de Opinión Empresarial ") + 
     geom_point(data = pick(~ Serie == "Desestacionalizada"),
                size = 1,
                shape=21, 
                fill = "#696969") +
  scale_color_manual(values=c("#696969","#A52A2A")) + 
  scale_linetype_manual(values=c("solid", …

我想知道如何在R中做这个图？

我考虑过使用ggplot包中的geom_bar.但是我不知道如何绘制变化的厚度和同时反映多年增长的颜色.

我很欣赏任何想法.

谢谢.