task1

2
setwd("C:/Users/fng134/Dropbox/KTH/Management/") ##install.packages("stats") library(stats) periodic.x=read.csv("X.csv") periodic.y=read.csv("Y.csv") y1=periodic.y$DispFrames y2=periodic.y$NoAudioPlayed y3=periodic.y$NoRTPPkts xx=periodic.x[c("all_. .usr","X..memused","X..s wpused","proc.s","cswch.s ","file.n r","sum_intr.s","rtps","ldavg.1","tcpsck","bread.s","pgfree.s")] x=model.matrix(~.,xx)[,-1] #y=subset(periodic.y, select = c("DispFrames") ) set.seed(0) train=sample(50000,35000) test=(-train) ##### Y1 prediction Exectime <- proc.time() lm.fit = lm(y1~ x,subset = train) xf=data.frame(x[test,]) xtest=x[test,] ypredic=predict(lm.fit,xf) proc.time() - Exectime plot(ypredic[test],col="blue") par=(new=T) lines(y1[test],col="red") resY1=(sum(abs(y1[test]-ypredic[test]))/15000)/mean(y1[test]) cat(resY1) ### For Y2 lm.fit = lm(y2~ x,subset = train) y2predic=predict(lm.fit,xf) plot(y2predic[test],col="blue") par=(new=T) lines(y2[test],col="red") res=(sum(abs(y2[test]-y2predic[test]))/15000)/mean(y2[test]) cat(res) ### For Y3

Upload: constantine

Post on 08-Jan-2016

212 views

Category:

Documents


0 download

DESCRIPTION

fgdfgdfg

TRANSCRIPT

Page 1: Task1

7/17/2019 Task1

http://slidepdf.com/reader/full/task1-568ef52055fdc 1/2

setwd("C:/Users/fng134/Dropbox/KTH/Management/")##install.packages("stats")library(stats)periodic.x=read.csv("X.csv")periodic.y=read.csv("Y.csv")

y1=periodic.y$DispFramesy2=periodic.y$NoAudioPlayedy3=periodic.y$NoRTPPkts

xx=periodic.x[c("all_..usr","X..memused","X..swpused","proc.s","cswch.s","file.nr","sum_intr.s","rtps","ldavg.1","tcpsck","bread.s","pgfree.s")]

x=model.matrix(~.,xx)[,-1]#y=subset(periodic.y, select = c("DispFrames") )

set.seed(0)train=sample(50000,35000)

test=(-train)

##### Y1 prediction

Exectime <- proc.time()

lm.fit = lm(y1~ x,subset = train)xf=data.frame(x[test,])xtest=x[test,]

ypredic=predict(lm.fit,xf)

proc.time() - Exectime

plot(ypredic[test],col="blue")

par=(new=T)lines(y1[test],col="red")

resY1=(sum(abs(y1[test]-ypredic[test]))/15000)/mean(y1[test])cat(resY1)

### For Y2

lm.fit = lm(y2~ x,subset = train)

y2predic=predict(lm.fit,xf)

plot(y2predic[test],col="blue")par=(new=T)lines(y2[test],col="red")

res=(sum(abs(y2[test]-y2predic[test]))/15000)/mean(y2[test])cat(res)

### For Y3

Page 2: Task1

7/17/2019 Task1

http://slidepdf.com/reader/full/task1-568ef52055fdc 2/2

lm.fit = lm(y3~ x,subset = train)

y3predic=predict(lm.fit,xf)

plot(y3predic[test],col="blue")par=(new=T)lines(y3[test],col="red")

resY3=(sum(abs(y3[test]-y3predic[test]))/15000)/mean(y3[test])cat(resY3)