multiscan (version 1.32.0)

plot.multiscan: Plot fitted model and standardised residuals for a multiscan object

Description

This function provides plots for fitted model and standardised residual for multiscan fitted objects.

Usage

"plot"(x, residual=FALSE,...)

Arguments

x

a multiscan fitted object.

residual

Logical. Should residuals be plotted instead of the fitted model?. Defaults to FALSE.

...

Further arguments passed to the plot function.

Value

Returns either one plot of the fitted model (residual=FALSE) on the input data or m plots of the residuals (residual=TRUE) corresponding to each scan of data.

Details

Fitted model on the input data, after rescaling by the corresponding scanning effects, is plotted against the estimated gene expressions. Standardised residuals for each scan of data are plotted against the rank of estimated gene expressions.

References

Khondoker, M. R., Glasbey, C. A. and Worton, B. J. (2006). Statistical estimation of gene expression using multiple laser scans of microarrays. Bioinformatics 22, 215--219.

See Also

multiscan

Examples

data(murine)

fit<-multiscan(murine)

## plot the fitted model

plot(fit)

## plot the residuals

op<-par(mfrow=c(2,2))
plot(fit, residual=TRUE)
par(op)