Methekar, Pushpa (GE Transportation, Non-GE)
2014-Dec-30 08:43 UTC
[R] to create a fuction for removing outlier
Hi all,
Thanks for your support. sorry for inconvenience cause due to lack of
information. here is my all r- program .
#### read in data
raw=read.csv(file.choose(),header = T, stringsAsFactors=FALSE)
str(raw)
####conversion from factor to numeric of xs
raw$Engine.Speed <- as.numeric(raw$Engine.Speed)
raw$Engine.Power<- as.numeric(raw$Engine.Power)
raw$Diesel.LHV <- as.numeric(raw$Diesel.LHV)
raw$NG.LHV <- as.numeric(raw$NG.LHV)
raw$Compressor.Inlet.Temperature <-
as.numeric(raw$Compressor.Inlet.Temperature)
raw$Int.Mfld.Temp<- as.numeric(raw$Int.Mfld.Temp)
raw$Cmd.Advance.Angle<- as.numeric(raw$Cmd.Advance.Angle)
raw$Cmd.Fuel.Rail.Press.Manual<- as.numeric(raw$Cmd.Fuel.Rail.Press.Manual)
raw$X..NG.by.Energy<-
as.numeric(sub("%","",raw$X..NG.by.Energy))/100
####conversion from factor to numeric of Ys
raw$Pre.Turb.Temp.L <- as.numeric(raw$Pre.Turb.Temp.L)
raw$Avg.PCP <- as.numeric(raw$Avg.PCP)
raw$BSHC <- as.numeric(raw$BSHC)
raw$BSCO <- as.numeric(raw$BSCO)
raw$GETS.cBSNOx <- as.numeric(raw$GETS.cBSNOx)
raw$Avg.COV.GIMEP <- as.numeric(raw$Avg.COV.GIMEP)
raw$Intake.Manifold.Pressure <- as.numeric(raw$Intake.Manifold.Pressure)
raw$Emiss..1..EPA.MAF..Dry. <-
as.numeric(sub(",","",raw$Emiss..1..EPA.MAF..Dry.))
raw$Turbo.Speed <- as.numeric(raw$Turbo.Speed)
#### create a matrix plot of scatter plot
xs=raw[c("Engine.Speed","X..NG.by.Energy","Int.Mfld.Temp","Cmd.Advance.Angle","Engine.Power","Cmd.Fuel.Rail.Press.Manual","Compressor.Inlet.Temperature","Diesel.LHV","NG.LHV")]
ys=raw[c("Pre.Turb.Temp.L","Avg.PCP","BSHC","BSCO","GETS.cBSNOx","Avg.COV.GIMEP","Turbo.Speed","Intake.Manifold.Pressure","Emiss..1..EPA.MAF..Dry.")]
pairs(xs,main="x's")
pairs(ys,main="y's")
### calculate correlation coefficient
xs[is.na(xs)] <- 0 #replace all NA with 0
cor(xs)
ys[is.na(ys)] <- 0 #replace all NA with 0
cor(ys)
cor(xs,ys,method = "spearman")
pairs(xs,ys)
save.image(file="plotting.RData")
## make data frame of x's and y's
xsys=data.frame(xs,ys)
library(MASS)
#### first order multiple linear regression
fitted.modely1.temp.l <- rlm(
Pre.Turb.Temp.L~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=xsys)
fitted.modely2.avg.pcp <- rlm(
Avg.PCP~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=xsys)
fitted.modely3.bshc <- rlm(
BSHC~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=xsys)
fitted.modely4.bsco <-
rlm(BSCO~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=xsys)
fitted.modely5.gets <-
rlm(GETS.cBSNOx~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=xsys)
fitted.modely6.gimep <-
rlm(Avg.COV.GIMEP~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=xsys)
fitted.modely7.ts <- rlm( Turbo.Speed
~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=xsys)
fitted.modely8.imp <- rlm(Intake.Manifold.Pressure
~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=xsys)
fitted.modely9.maf <- rlm( Emiss..1..EPA.MAF..Dry.
~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=xsys)
#### comlete second order multiple linear regression
secfitted.modely1.temp.l<-lm(Pre.Turb.Temp.L
~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Manual,2)
+poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
+Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Speed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV+NG.LHV)
+Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed+NG.LHV)
+NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engine.Speed),data=xsys)
secfitted.modely2.avg.pcp<-lm(Avg.PCP
~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Manual,2)
+poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
+Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Speed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV+NG.LHV)
+Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed+NG.LHV)
+NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engine.Speed),data=xsys)
secfitted.modely3.bshc<-lm( BSHC
~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Manual,2)
+poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
+Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Speed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV+NG.LHV)
+Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed+NG.LHV)
+NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engine.Speed),data=xsys)
secfitted.modely4.bsco<-lm( BSCO
~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Manual,2)
+poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
+Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Speed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV+NG.LHV)
+Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed+NG.LHV)
+NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engine.Speed),data=xsys)
secfitted.modely5.gets<-lm( GETS.cBSNOx
~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Manual,2)
+poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
+Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Speed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV+NG.LHV)
+Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed+NG.LHV)
+NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engine.Speed),data=xsys)
secfitted.modely6.gimep<-lm( Avg.COV.GIMEP
~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Manual,2)
+poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
+Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Speed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV+NG.LHV)
+Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed+NG.LHV)
+NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engine.Speed),data=xsys)
secfitted.modely7.ts<-lm( Turbo.Speed
~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Manual,2)
+poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
+Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Speed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV+NG.LHV)
+Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed+NG.LHV)
+NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engine.Speed),data=xsys)
secfitted.modely8.imp<-lm( Intake.Manifold.Pressure
~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Manual,2)
+poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
+Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Speed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV+NG.LHV)
+Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed+NG.LHV)
+NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engine.Speed),data=xsys)
secfitted.modely9.maf<-lm(Emiss..1..EPA.MAF..Dry.~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Manual,2)
+poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
+Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Speed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV)
+Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV+NG.LHV)
+Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed+NG.LHV)
+NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engine.Speed),data=xsys)
####stepwise forward regression for first order
step(fitted.modely1.temp.l,direction="forward")
step(fitted.modely2.avg.pcp,direction="forward")
step(fitted.modely3.bshc,direction="forward")
step(fitted.modely4.bsco,direction="forward")
step(fitted.modely5.gets,direction="forward")
step(fitted.modely6.gimep,direction="forward")
step(fitted.modely7.ts,direction="forward")
step(fitted.modely8.imp,direction="forward")
step(fitted.modely9.maf,direction="forward")
step(fitted.modely1.temp.l,direction="backward")
step(fitted.modely2.avg.pcp,direction="backward")
step(fitted.modely3.bshc,direction="backward")
step(fitted.modely4.bsco,direction="backward")
step(fitted.modely5.gets,direction="backward")
step(fitted.modely6.gimep,direction="backward")
step(fitted.modely7.ts,direction="backward")
step(fitted.modely8.imp,direction="backward")
step(fitted.modely9.maf,direction="backward")
now after this step I want to find out outliers in each fitted model and remove
them from xsys data set and update it .
Ii want a function which will find one by one outlier and remove it .
The condition is that outlier would be consider as outlier if its r student I
greater than +3 and -3.
I want all in loop and repeat loop till no outlier is remained.
Hi
Without going to deeply to your code; your data seems to be numeric, but for
some reason you read them as character. You shall consult read.csv help and your
original data. They probably have some feature which prevent them to be read
correctly.
To your question, what you want can be achieved by while loop.
something like this
while (abs(studres(object))>3) {
ms <-max(abs(studres(object)))
del <- which(ms== abs(studres(object)))
#do fitting with del removed from data
#save the fit
}
you can also use rstudent(object)
Cheers
Petr
> -----Original Message-----
> From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of
> Methekar, Pushpa (GE Transportation, Non-GE)
> Sent: Tuesday, December 30, 2014 9:44 AM
> To: r-help at r-project.org
> Subject: [R] to create a fuction for removing outlier
>
> Hi all,
> Thanks for your support. sorry for inconvenience cause due to lack of
> information. here is my all r- program .
>
>
> #### read in data
> raw=read.csv(file.choose(),header = T, stringsAsFactors=FALSE)
>
> str(raw)
>
> ####conversion from factor to numeric of xs raw$Engine.Speed <-
> as.numeric(raw$Engine.Speed)
> raw$Engine.Power<- as.numeric(raw$Engine.Power) raw$Diesel.LHV <-
> as.numeric(raw$Diesel.LHV) raw$NG.LHV <- as.numeric(raw$NG.LHV)
> raw$Compressor.Inlet.Temperature <-
> as.numeric(raw$Compressor.Inlet.Temperature)
> raw$Int.Mfld.Temp<- as.numeric(raw$Int.Mfld.Temp)
> raw$Cmd.Advance.Angle<- as.numeric(raw$Cmd.Advance.Angle)
> raw$Cmd.Fuel.Rail.Press.Manual<-
> as.numeric(raw$Cmd.Fuel.Rail.Press.Manual)
> raw$X..NG.by.Energy<-
as.numeric(sub("%","",raw$X..NG.by.Energy))/100
>
>
> ####conversion from factor to numeric of Ys raw$Pre.Turb.Temp.L <-
> as.numeric(raw$Pre.Turb.Temp.L) raw$Avg.PCP <- as.numeric(raw$Avg.PCP)
> raw$BSHC <- as.numeric(raw$BSHC) raw$BSCO <- as.numeric(raw$BSCO)
> raw$GETS.cBSNOx <- as.numeric(raw$GETS.cBSNOx) raw$Avg.COV.GIMEP <-
> as.numeric(raw$Avg.COV.GIMEP) raw$Intake.Manifold.Pressure <-
> as.numeric(raw$Intake.Manifold.Pressure)
> raw$Emiss..1..EPA.MAF..Dry. <-
> as.numeric(sub(",","",raw$Emiss..1..EPA.MAF..Dry.))
> raw$Turbo.Speed <- as.numeric(raw$Turbo.Speed)
>
>
> #### create a matrix plot of scatter plot
>
xs=raw[c("Engine.Speed","X..NG.by.Energy","Int.Mfld.Temp","Cmd.Advance.
>
Angle","Engine.Power","Cmd.Fuel.Rail.Press.Manual","Compressor.Inlet.Te
> mperature","Diesel.LHV","NG.LHV")]
>
>
ys=raw[c("Pre.Turb.Temp.L","Avg.PCP","BSHC","BSCO","GETS.cBSNOx","Avg.C
>
OV.GIMEP","Turbo.Speed","Intake.Manifold.Pressure","Emiss..1..EPA.MAF..
> Dry.")]
> pairs(xs,main="x's")
> pairs(ys,main="y's")
>
> ### calculate correlation coefficient
> xs[is.na(xs)] <- 0 #replace all NA with 0
> cor(xs)
> ys[is.na(ys)] <- 0 #replace all NA with 0
> cor(ys)
> cor(xs,ys,method = "spearman")
> pairs(xs,ys)
> save.image(file="plotting.RData")
>
>
> ## make data frame of x's and y's
> xsys=data.frame(xs,ys)
>
> library(MASS)
> #### first order multiple linear regression fitted.modely1.temp.l <-
> rlm(
> Pre.Turb.Temp.L~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.
> Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LH
> V+NG.LHV,data=xsys)
>
> fitted.modely2.avg.pcp <- rlm(
> Avg.PCP~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cm
> d.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV
> ,data=xsys)
>
> fitted.modely3.bshc <- rlm(
> BSHC~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.F
> uel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,da
> ta=xsys)
>
> fitted.modely4.bsco <-
> rlm(BSCO~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LH
> V,data=xsys)
>
> fitted.modely5.gets <-
> rlm(GETS.cBSNOx~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.
> Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LH
> V+NG.LHV,data=xsys)
>
> fitted.modely6.gimep <-
> rlm(Avg.COV.GIMEP~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engin
> e.Power+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.
> LHV+NG.LHV,data=xsys)
>
> fitted.modely7.ts <- rlm( Turbo.Speed
> ~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.
> Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=x
> sys)
>
> fitted.modely8.imp <- rlm(Intake.Manifold.Pressure
> ~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.
> Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=x
> sys)
>
> fitted.modely9.maf <- rlm( Emiss..1..EPA.MAF..Dry.
> ~X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+Cmd.Fuel.
> Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+NG.LHV,data=x
> sys)
>
>
>
> #### comlete second order multiple linear regression
> secfitted.modely1.temp.l<-lm(Pre.Turb.Temp.L
> ~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+pol
> y(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Ma
> nual,2)
>
> +poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
>
> +Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.S
> peed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.
> Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advanc
> e.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV
> +NG.LHV)
>
> +Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Pow
> er+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed
> +NG.LHV)
>
> +NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engin
> e.Speed),data=xsys)
>
> secfitted.modely2.avg.pcp<-lm(Avg.PCP
> ~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+pol
> y(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Ma
> nual,2)
>
> +poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
>
> +Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.S
> peed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.
> Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advanc
> e.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV
> +NG.LHV)
>
> +Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Pow
> er+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed
> +NG.LHV)
>
> +NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engin
> e.Speed),data=xsys)
>
> secfitted.modely3.bshc<-lm( BSHC
> ~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+pol
> y(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Ma
> nual,2)
>
> +poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
>
> +Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.S
> peed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.
> Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advanc
> e.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV
> +NG.LHV)
>
> +Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Pow
> er+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed
> +NG.LHV)
>
> +NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engin
> e.Speed),data=xsys)
>
>
> secfitted.modely4.bsco<-lm( BSCO
> ~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+pol
> y(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Ma
> nual,2)
>
> +poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
>
> +Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.S
> peed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.
> Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advanc
> e.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV
> +NG.LHV)
>
> +Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Pow
> er+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed
> +NG.LHV)
>
> +NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engin
> e.Speed),data=xsys)
>
>
>
> secfitted.modely5.gets<-lm( GETS.cBSNOx
> ~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+pol
> y(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Ma
> nual,2)
>
> +poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
>
> +Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.S
> peed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.
> Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advanc
> e.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV
> +NG.LHV)
>
> +Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Pow
> er+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed
> +NG.LHV)
>
> +NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engin
> e.Speed),data=xsys)
>
>
>
> secfitted.modely6.gimep<-lm( Avg.COV.GIMEP
> ~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+pol
> y(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Ma
> nual,2)
>
> +poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
>
> +Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.S
> peed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.
> Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advanc
> e.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV
> +NG.LHV)
>
> +Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Pow
> er+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed
> +NG.LHV)
>
> +NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engin
> e.Speed),data=xsys)
>
> secfitted.modely7.ts<-lm( Turbo.Speed
> ~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+pol
> y(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Ma
> nual,2)
>
> +poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
>
> +Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.S
> peed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.
> Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advanc
> e.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV
> +NG.LHV)
>
> +Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Pow
> er+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed
> +NG.LHV)
>
> +NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engin
> e.Speed),data=xsys)
>
>
>
> secfitted.modely8.imp<-lm( Intake.Manifold.Pressure
> ~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)+poly(Int.Mfld.Temp,2)+pol
> y(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+poly(Cmd.Fuel.Rail.Press.Ma
> nual,2)
>
> +poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
>
> +Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.S
> peed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.
> Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advanc
> e.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV
> +NG.LHV)
>
> +Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Pow
> er+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed
> +NG.LHV)
>
> +NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engin
> e.Speed),data=xsys)
>
>
> secfitted.modely9.maf<-
> lm(Emiss..1..EPA.MAF..Dry.~poly(Engine.Speed,2)+poly(X..NG.by.Energy,2)
> +poly(Int.Mfld.Temp,2)+poly(Cmd.Advance.Angle,2)+poly(Engine.Power,2)+p
> oly(Cmd.Fuel.Rail.Press.Manual,2)
>
> +poly(Compressor.Inlet.Temperature,2)+poly(Diesel.LHV,2)+poly(NG.LHV,2)
>
> +Engine.Speed*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +X..NG.by.Energy*(Engine.Speed+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Int.Mfld.Temp*(X..NG.by.Energy+Engine.Speed+Cmd.Advance.Angle+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Advance.Angle*(X..NG.by.Energy+Int.Mfld.Temp+Engine.Speed+Engine.P
> ower+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Engine.Power*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.S
> peed+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Cmd.Fuel.Rail.Press.Manual*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.
> Angle+Engine.Power+Engine.Speed+Compressor.Inlet.Temperature+Diesel.LHV
> +NG.LHV)
>
> +Compressor.Inlet.Temperature*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advanc
> e.Angle+Engine.Power+Cmd.Fuel.Rail.Press.Manual+Engine.Speed+Diesel.LHV
> +NG.LHV)
>
> +Diesel.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Pow
> er+Cmd.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Engine.Speed
> +NG.LHV)
>
> +NG.LHV*(X..NG.by.Energy+Int.Mfld.Temp+Cmd.Advance.Angle+Engine.Power+C
> md.Fuel.Rail.Press.Manual+Compressor.Inlet.Temperature+Diesel.LHV+Engin
> e.Speed),data=xsys)
>
>
> ####stepwise forward regression for first order
>
> step(fitted.modely1.temp.l,direction="forward")
> step(fitted.modely2.avg.pcp,direction="forward")
> step(fitted.modely3.bshc,direction="forward")
> step(fitted.modely4.bsco,direction="forward")
> step(fitted.modely5.gets,direction="forward")
> step(fitted.modely6.gimep,direction="forward")
> step(fitted.modely7.ts,direction="forward")
> step(fitted.modely8.imp,direction="forward")
> step(fitted.modely9.maf,direction="forward")
>
>
> step(fitted.modely1.temp.l,direction="backward")
> step(fitted.modely2.avg.pcp,direction="backward")
> step(fitted.modely3.bshc,direction="backward")
> step(fitted.modely4.bsco,direction="backward")
> step(fitted.modely5.gets,direction="backward")
> step(fitted.modely6.gimep,direction="backward")
> step(fitted.modely7.ts,direction="backward")
> step(fitted.modely8.imp,direction="backward")
> step(fitted.modely9.maf,direction="backward")
>
>
> now after this step I want to find out outliers in each fitted model
> and remove them from xsys data set and update it .
> Ii want a function which will find one by one outlier and remove it .
> The condition is that outlier would be consider as outlier if its r
> student I greater than +3 and -3.
> I want all in loop and repeat loop till no outlier is remained.
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.
________________________________
Tento e-mail a jak?koliv k n?mu p?ipojen? dokumenty jsou d?v?rn? a jsou ur?eny
pouze jeho adres?t?m.
Jestli?e jste obdr?el(a) tento e-mail omylem, informujte laskav? neprodlen? jeho
odes?latele. Obsah tohoto emailu i s p??lohami a jeho kopie vyma?te ze sv?ho
syst?mu.
Nejste-li zam??len?m adres?tem tohoto emailu, nejste opr?vn?ni tento email
jakkoliv u??vat, roz?i?ovat, kop?rovat ?i zve?ej?ovat.
Odes?latel e-mailu neodpov?d? za eventu?ln? ?kodu zp?sobenou modifikacemi ?i
zpo?d?n?m p?enosu e-mailu.
V p??pad?, ?e je tento e-mail sou??st? obchodn?ho jedn?n?:
- vyhrazuje si odes?latel pr?vo ukon?it kdykoliv jedn?n? o uzav?en? smlouvy, a
to z jak?hokoliv d?vodu i bez uveden? d?vodu.
- a obsahuje-li nab?dku, je adres?t opr?vn?n nab?dku bezodkladn? p?ijmout;
Odes?latel tohoto e-mailu (nab?dky) vylu?uje p?ijet? nab?dky ze strany p??jemce
s dodatkem ?i odchylkou.
- trv? odes?latel na tom, ?e p??slu?n? smlouva je uzav?ena teprve v?slovn?m
dosa?en?m shody na v?ech jej?ch n?le?itostech.
- odes?latel tohoto emailu informuje, ?e nen? opr?vn?n uzav?rat za spole?nost
??dn? smlouvy s v?jimkou p??pad?, kdy k tomu byl p?semn? zmocn?n nebo p?semn?
pov??en a takov? pov??en? nebo pln? moc byly adres?tovi tohoto emailu p??padn?
osob?, kterou adres?t zastupuje, p?edlo?eny nebo jejich existence je adres?tovi
?i osob? j?m zastoupen? zn?m?.
This e-mail and any documents attached to it may be confidential and are
intended only for its intended recipients.
If you received this e-mail by mistake, please immediately inform its sender.
Delete the contents of this e-mail with all attachments and its copies from your
system.
If you are not the intended recipient of this e-mail, you are not authorized to
use, disseminate, copy or disclose this e-mail in any manner.
The sender of this e-mail shall not be liable for any possible damage caused by
modifications of the e-mail or by delay with transfer of the email.
In case that this e-mail forms part of business dealings:
- the sender reserves the right to end negotiations about entering into a
contract in any time, for any reason, and without stating any reasoning.
- if the e-mail contains an offer, the recipient is entitled to immediately
accept such offer; The sender of this e-mail (offer) excludes any acceptance of
the offer on the part of the recipient containing any amendment or variation.
- the sender insists on that the respective contract is concluded only upon an
express mutual agreement on all its aspects.
- the sender of this e-mail informs that he/she is not authorized to enter into
any contracts on behalf of the company except for cases in which he/she is
expressly authorized to do so in writing, and such authorization or power of
attorney is submitted to the recipient or the person represented by the
recipient, or the existence of such authorization is known to the recipient of
the person represented by the recipient.