Description Usage Arguments Details Value Note Author(s) References See Also Examples

The parameters of the random coefficients model are estimated by means of the maximum likelihood method. The implemented maximum likelihood procedure has been optimized with respect to computational efficiency and memory usage.

1 2 3 | ```
RCMestimation(Y, X, R, hypothesis = "H2", shrinkType = "none",
estType = "normal", corType = "unif", maxNoIt = 100,
minSuccDist = 0.005, verbose = FALSE)
``` |

`Y` |
The |

`X` |
The design matrix (number of rows equal to number of samples, number of columns equal to number of covariates). |

`R` |
The linear constraint matrix (number of columns equal to the number of covariates). |

`hypothesis` |
The hypothesis under which the model is fitted: |

`shrinkType` |
The type of shrinkage to be applied to the error variances: |

`estType` |
Type of estimation, either |

`corType` |
Correlation structure to be used, either |

`maxNoIt` |
Maximum number of iterations in the ML procedure. |

`minSuccDist` |
Minimum distance between estimates of two successive iterations to be achieved. |

`verbose` |
Logical indicator: should intermediate output be printed on the screen? |

Details on the type of random coefficients model that is actually fitted are specified in the reference below.

Object of class `rcmFit`

.

In case a covariate for the intercept is included in the design matrix `X`

we strongly recommend the center, per feature, the data around zero.

Wessel N. van Wieringen: [email protected]

Van Wieringen, W.N., Berkhof, J., Van de Wiel, M.A. (2010), "A random coefficients model for regional co-expression associated with DNA copy number", *Statistical Applications in Genetics and Molecular Biology*, Volume 9, Issue1, Article 25, 1-28.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | ```
# load data
data(pollackCN16)
data(pollackGE16)
# select features belonging to a region
ids <- getSegFeatures(20, pollackCN16)
# extract segmented log2 ratios of the region
X <- t(segmented(pollackCN16)[ids[1], , drop=FALSE])
# extract segmented log2 ratios of the region
Y <- exprs(pollackGE16)[ids,]
# center the expression data (row-wise)
Y <- t(Y - apply(Y, 1, mean))
# specify the linear constraint matrix
R <- matrix(1, nrow=1)
# fit the random coefficients model to the random data
RCMresults <- RCMestimation(Y, X, R)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.