trimProjReg2d: trimProjReg2d

Description Usage Arguments Author(s) Examples

View source: R/trimprojreg2d.R

Description

Computes projection trimmed regression in 2 dimensions.

Usage

1
trimProjReg2d(x, y, alpha = 0.1)

Arguments

x

Independent variable

y

Dependent variable

alpha

Percentage of trimmed observations

Author(s)

Zygmunt Zawadzki from Cracow University of Economics.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
#EXAMPLE 1
 data(pension)
 plot(pension)
 abline(lm(Reserves~Income,data = pension), lty = 3, lwd = 2) #lm
 abline(trimProjReg2d(pension[,1],pension[,2]), lwd = 2) #trimprojreg2d
 legend("bottomright", c("OLS","TrimLS"), lty = 1:2)

 #EXAMPLE 2
 data(under5.mort)
 data(inf.mort)
 data(maesles.imm)

 data2011=na.omit(cbind(under5.mort[,22],inf.mort[,22],maesles.imm[,22]))
 x<-data2011[,3]
 y<-data2011[,2]
 plot(x,y,cex=1.2, ylab="infant mortality rate per 1000 live birth",
 xlab="against masles immunized #'
 percentage", main='Projection Depth Trimmed vs. LS regressions')
 abline(lm(x~y,data = pension), lwd = 2, col='black') #lm
 abline(trimProjReg2d(x,y), lwd = 2,col='red') #trimmed reg
 legend("bottomleft",c("LS","TrimReg"),fill=c("black","red"),cex=1.4,bty="n")

DepthProc documentation built on May 31, 2017, 4:49 a.m.