fit.li.wong(data.matrix, remove.outliers=TRUE, normal.array.quantile=0.5, normal.resid.quantile=0.9, large.threshold=3, large.variation=0.8, outlier.fraction=0.14, delta=1e-06, maxit=50, outer.maxit=50,verbose=FALSE, ...)
li.wong(data.matrix,remove.outliers=TRUE, normal.array.quantile=0.5, normal.resid.quantile=0.9, large.threshold=3, large.variation=0.8, outlier.fraction=0.14, delta=1e-06, maxit=50, outer.maxit=50,verbose=FALSE)
normal.array.quantile
of all SDs x
large.threshold
.normal.resid.quantile
quantile of all residuals x
large.threshold
is considered an outlier.TRUE
information is given of
the status of the algorithm.express
which is no longer part of
the package.fit.li.wong
returns much more. Namely, a list containing the
fitted parameters and relevant information.
TRUE
).TRUE
)FALSE
the algorithm did
not converge when fitting the phis and thetas.FALSE
the algorithm did
not converge in deciding what are outliers.Notice that this iterative algorithm will not always converge. If you run the algorithm on thousands of probes expect some non-convergence warnings. These are more likely when few arrays are used. We recommend using this method only if you have 10 or more arrays.
Please refer to references for more details.
Li, C. and Wong, W.H. (2001) Proc. Natl. Acad. Sci USA 98, 31--36.
li.wong
, expresso
x <- sweep(matrix(2^rnorm(600),30,20),1,seq(1,2,len=30),FUN="+")
fit1 <- fit.li.wong(x)
plot(x[1,])
lines(fit1$theta)
Run the code above in your browser using DataCamp Workspace