[R] help me
kafi dano
kafi_dano at yahoo.com
Sun Apr 27 13:30:01 CEST 2014
please help me and correct why this program give me error
bcos I want to apply the condition follow the attached below
# program (R) used estimate beta of robust ridge (weighted stacklose)
rm(list=ls())
library(MASS)
library(mvoutlier)
library(robustbase)
library(car)
library(quantreg)
n<-21
p<-3
data<-read.csv('C:/Users/Kafi/Documents/stacklose.csv',header=TRUE)
cell<-data
y<-as.matrix(cell[,1])
x1<-cell[,2]
x2<-cell[,3]
x3<-cell[,4]
stack<- as.matrix(data.frame(y,x1,x2,x3))
stack
d <- lm(y ~ x1+x2+x3)$residuals
d1<-lm(y ~ x1+x2+x3)$fit
d2<-abs(d)
k<-4.685
w<-NULL
for(i in 1:n){
if(d2<=k){w[i]=d{1-(d/k)^2}^2} else{w[i]=0}}
w
Kafi Dano Pati
Ph.D candidate ( mathematics/statistics)
Department of mathematical Science/ faculty of Science
University Technology Malaysia
81310 UTM, Johor Bahru, Johor, Malaysia
IC. NO. 201202F10234
Matric No. PS113113
HP. No. 00601117517559
E-mail: kafi_dano at yahoo.com
supervisor- Assoc. Prof. Robiah Binti Adnan
More information about the R-help
mailing list