iterGlm {mseq}R Documentation

iterative generalized linear regression

Description

Fit the Poisson linear model by iteratively fitting glm and gene expression levels.

Usage

iterGlm(data, thrd = 0.01, max_iter = 10)

Arguments

data the data frame generated by expData.R
thrd the threshould of improvement of R squared
max_iter maximum number of iterations

Value

an glm structure. Summarize it will list the coefficients and their p-values.

Examples

 data(g1_part) # for real data, please use read.csv, like g1 <- read.csv("g1.csv")
 data <- expData(g1_part, 2, 3)

 data.glm <- iterGlm(data)
 summary(data.glm)

[Package mseq version 1.1 Index]