比赛中用到的一些R语言命令

最新推荐文章于 2024-01-20 12:42:56 发布

Jerk_Lee

最新推荐文章于 2024-01-20 12:42:56 发布

阅读量1k

点赞数

文章标签： r语言

本文链接：https://blog.csdn.net/Jerk_Lee/article/details/46654505

版权

test<-scan("F:/1.txt")
demo<-ts(test,start=c(1))
demodiff1<-diff(demo,difference=1)
plot(demodiff1)
acf(demodiff1,lag.max=50,plot=FALSE)
pacf(demodiff1,lag.max=50,plot=FALSE)
r<-arima(demo,order=c(28,1,19))
s<-forecast.Arima(r,h=30,level=c(99.5))
r=arima(demo,order=c(28,1,34),list(order=c(0,1,2),period=7))
r=arima(demo,order=c(13,1,28),list(order=c(7,1,12)))

s<-forecast.Arima(r,h=30,level=c(99.5))
write.csv(s, file = "F:/555.csv", row.names = F, quote = F)

dataordered<-data[order(data[,2],data[,1]),]

user<- read.csv("F:/新数据/user1.csv",header=F)
data<-read.csv("F:/新数据/class.csv",header=F)
datap <- do.call("paste",data[1])
userp <- do.call("paste",user)
exist <- ! datap %in% userp
userN <- data[exist,]
userY <- data[!exist,]
write.csv(userN, file = "F:/非土豪赎回.csv", row.names = F, quote = F)
write.csv(userY, file = "F:/土豪赎回.csv", row.names = F, quote = F)

test<-scan("F:/1.txt")
demo<-ts(test,start=c(1))
lg<-log10(demo)
lgd1<-diff(lg,difference=1)
plot(lgd1)

r=arima(demo,order=c(26,1,26),list(order=c(9,1,14)))

acf(lgd1,lag.max=50,plot=FALSE)