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@Epsi2(b) B <- hubPsi@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: arasheedbe...@yahoo.com supervisor- Assoc. Prof. Robiah Binti Adnan [[alternative HTML version deleted]] ______________________________________________ R-help@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.