- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 18999 个
- 通用积分
- 1029.7754
- 学术水平
- 146 点
- 热心指数
- 166 点
- 信用等级
- 135 点
- 经验
- 36357 点
- 帖子
- 541
- 精华
- 0
- 在线时间
- 887 小时
- 注册时间
- 2015-9-25
- 最后登录
- 2025-12-4
|
|
- f=read.table('ggdata.txt',header=TRUE)[,-1]
- library(tidyverse)
- library(ggplot2)
- library(splines)
- library(ggthemes)
- #ggplot(f,aes(x=TI,y=RAIN,colour=MODEL,group=MODEL))+geom_line()
- #ggbs <- lm(RAIN ~ bs(TI, df = 3), data = f)
- data1=spread(f,MODEL,RAIN)
- ggbs1 <- lm(obs ~ bs(TI, df = 10), data = data1)
- ggbs2 <- lm( rr ~ bs(TI, df = 10), data = data1)
- ggbs3 <- lm( xaj ~ bs(TI, df = 10), data = data1)
- TI1=seq(1,length(data1$TI),length.out=100)
- RAIN1=c(predict(ggbs1,data.frame(TI=TI1)),predict(ggbs2,data.frame(TI=TI1)),
- predict(ggbs3,data.frame(TI=TI1)))
- MODEL1=rep(names(data1)[-1],rep(100,3))
- data2=data.frame(TI=rep(TI1,3),RAIN=RAIN1,MODEL=MODEL1)
- ggplot(data2,aes(x=TI,y=RAIN,colour=MODEL,group=MODEL))+geom_line()+
- theme_wsj()
复制代码
|
-
总评分: 论坛币 + 10
查看全部评分
|