R语言:我写的一个数据可视化的函数,散点图和线性回归趋势线及公式

在做项目的时候,发现需要多次绘制散点图和图例,以及计算线性回归的拟合公式和R2,为了方便处理,特编写了这个函数,希望对大家有所帮助。
也欢迎大家提出报给的改进建议或者更多的需求组合。具体的可视化函数如下。
其中,data.xts是tidy的时间数据集,调用时只需要给出time和lie这2个变量,lie的格式是data.xts的第几个字段。sandiantu函数直接绘制散点图
qushixian函数将添加趋势线,并显示公式和R2,在Rmd格式中使用效果更好。
sandiantu<-function(time,lie){
        a<-names(data.xts)[lie][1]
        b<-names(data.xts)[lie][2]
        with(as.dataframe(data.xts[time,lie]),
             {
                     plot(get(b) ~ get(a),xlab=a,ylab=b,col='lightgray',lwd=1)
             })
}

qushixian<-function(time,lie){ 
        a<-names(data.xts)[lie][1]
        b<-names(data.xts)[lie][2]
        with(as.dataframe(data.xts[time,lie]),

             {
                      reg<-lm((get(b)~get(a)))
                     abline(reg,col="red",lty=3,lwd=3)
                     gongshi<-paste("y=",substr(summary(reg)$coef[2,1],1,6),"x+",
                                    substr(summary(reg)$coef[1,1],1,6),sep='')
                     rfang<-paste("R2=",substr(summary(reg)$r.squared,1,6),sep='')
                     gongshi2<-paste("y=",summary(reg)$coef[2,1],"x+",
                                     summary(reg)$coef[1,1],sep='')
                     legend("topright", legend = c(gongshi,rfang))
                     print(gongshi2)
                     print(rfang)
             }

        )

}
  • 5
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值