лул, R - rextester

лул

dat<-read.table(file("stdin"))
dat
t<-as.numeric(dat[,1])
x<-as.double(dat[,2])
t
x
f<-as.data.frame(dat)
f
n<-length(x)
n
averageLevel<-mean(x)
averageLevel
for( i in 1:n) if ( x[i] < averageLevel ) f[i,3] = 0 else f[i,3] = 1
f
lowbound<-(n+1-3*sqrt(n-1))/2
upperbound<-(n+1+3*sqrt(n-1))/2
lowbound
upperbound
f[1,4]=(5*f[1,2]+2*f[2,2]-f[3,2])/6
f[n,4]=(5*f[n,2]+2*f[n-1,2]-f[n-2,2])/6
for ( i in 3:n-1 ) f[i, 4]=(x[i-1]+x[i]+x[i+1])/3
f
f[1,5] <- 
for ( j in 2:n ) f[j,5]=f[j-1,5]+1
for ( j in 1:n ) f[j,6]=f[j,5]*f[j,5]
for ( j in 1:n ) f[j,7]=f[j,2]*f[j,5]
f
a0<-sum(f[,2]) / n
a1<-sum(f[,7])/sum(f[,6])
a0
a1
for ( j in 1:n ) f[j,8]=f[j,5]*a1 + a0
f
prognoz <- a1*9+ a0
prognoz
x<-as.numeric(f[,1])
y<-as.numeric(f[,2])
z<-as.numeric(f[,4])
u<-as.numeric(f[,8])
# ËÒıÓ‰Ì˚Â ‰‡ÌÌ˚Â
plot(x, y, ylim=range(c(y,z, u)), type="l", col="red", lwd = 2)
par(new = TRUE)
# ÚÂıÛÓ‚ÌÂ‚‡ˇ ÒÍÓÎ¸Áˇ˘‡ˇ ÒÂ‰Ìˇˇ
plot(x, z, ylim=range(c(y,z,u)), axes = FALSE, type="l",  col="blue", xlab = "", ylab = "", lwd = 2)
par(new = TRUE)
# ÏÂÚÓ‰ ‡Ì‡ÎËÚË˜ÂÒÍÓ„Ó ‚˚‡‚ÌË‚‡ÌËˇ
plot(x, u, ylim=range(c(y,z,u)), axes = FALSE, type="l",  col="green", xlab = "", ylab = "", lwd = 2)


run \| edit \| history \| help	0

Please log in to post a comment.