[R] R-Help

[Bello A rasheed]

Dear R users 

Kindly help me on how run this code because there is error when ever i run it.
It give me an error at Error: lent >= 1 is not TRUE 
Error: object 'hubPsi' not found
Error in cbind(b = b, A = A, B = B, V = A/B^2, e = B^2/A, gamma. = b/B,  :   dims [product 3] do not match the length of object [21] Thank you in advance********************************************************
R-CODE
Hf0 <- function(x, c=1.35) pmax2(-c,pmin2(x,c))

Hf <- new("functionX",Hf0)

stopifnot(validObject(Hf)) # ok !


 

 
### psiFunc() examples## classical {trivial, notinteresting}:

F1 <- function(x) rep.int(1,length(x))

cPsi <- psiFunc(rho = function(x)x^2 / 2, psi = function(x) x,

                wgt = F1, Dpsi = F1,

                Erho = function(x) rep.int(1/2,length(x)),

                Epsi2 = F1, EDpsi = F1)


 
### MASS  -- ?rlm --- has

##

##-     psi.huber   (u, k = 1.345, deriv =0)

##-     psi.hampel  (u, a = 2, b = 4, c =8, deriv = 0)

##-     psi.bisquare(u, c = 4.685, deriv = 0)

##  where deriv = 0 :  psi(x)/xi.e.  'wgt'


 
## MM has more in psi-funs.R (seeabove)

## Reproduce Table 1, p.138 of Hampel,Rousseeuw, Ronchetti, Stahel (1986):
b <- c(seq(0, 3, by = 0.1), 4, 5,Inf)

A <- hubPsi at Epsi2(b)

B <- hubPsi at EDpsi(b)


 
options(width = 100, digits = 7)


 
cbind(b = b, A = A, B = B, V = A/B^2, e= B^2/A,

     gamma. = b/B, k. = 1 + b^2/A, lambda. = 1/B)


 
(hubPsi2 <- chgDefaults(hubPsi, k =2)) # works too!

 Bello Abdulkadir RasheedPh.D candidate ( mathematics/statistics)
Department of mathematical Science/ faculty of Science
University Technology Malaysia
81310 UTM, Johor Bahru, Johor, Malaysia
IC. NO. 201204M10026
Matric No. PS113124
HP. No.  0060149144837E-mail: arasheedbello at yahoo.com
supervisor- Assoc. Prof. Robiah Binti Adnan
	[[alternative HTML version deleted]]