Fit Intercept to Vector of Gamma Distributed Variates
Fit Intercept to Vector of Gamma Distributed Variates
fitGammaIntercept(y,offset=0,maxit=1000)
y |
numeric vector of positive response values. |
offset |
numeric vector giving known part of the expected value of |
maxit |
maximum number of Newton iterations to be done. |
The values y are assumed to follow a gamma distribution with common shape parameter and with expected values given by x+offset.
The function implements a globally convergent Newton iteration to estimate x.
Numeric value giving intercept.
Gordon Smyth and Belinda Phipson
Phipson, B. (2013). Empirical Bayes modelling of expression profiles and their associations. PhD Thesis. University of Melbourne, Australia.
This function is called by genas.
offset <- runif(10) x <- 9 mu <- x+offset y <- rgamma(10,shape=20,scale=mu/20) fitGammaIntercept(y,offset=offset)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.