similar to: Formula 1.0-0: Model formulas with multiple parts and responses

This update adds an intercept option (by popular request) - now one can fit a model without an intercept Glmnet is a package that fits the regularization path for a number of generalized linear models, with with "elastic net" regularization (tunable mixture of L1 and L2 penalties). Glmnet uses pathwise coordinate descent, and is very fast. The current list of models covered are:

This update adds an intercept option (by popular request) - now one can fit a model without an intercept Glmnet is a package that fits the regularization path for a number of generalized linear models, with with "elastic net" regularization (tunable mixture of L1 and L2 penalties). Glmnet uses pathwise coordinate descent, and is very fast. The current list of models covered are:

Is there a way to specify a Z matrix using the lmer function, where the model is written as y = X*Beta + Z*u + e? I am trying to reproduce smoothing methods illustrated in the paper "Smoothing with Mixed Model Software" my Long Ngo and M.P. Wand. published in the /Journal of Statistical Software/ in 2004 using the lme4 and Matrix packages. The code and data sets used can be found at

# plyr plyr is a set of tools for a common set of problems: you need to __split__ up a big data structure into homogeneous pieces, __apply__ a function to each piece and then __combine__ all the results back together. For example, you might want to: * fit the same model each patient subsets of a data frame * quickly calculate summary statistics for each group * perform group-wise

# plyr plyr is a set of tools for a common set of problems: you need to __split__ up a big data structure into homogeneous pieces, __apply__ a function to each piece and then __combine__ all the results back together. For example, you might want to: * fit the same model each patient subsets of a data frame * quickly calculate summary statistics for each group * perform group-wise

Dear All, I revise my question about the problem I have. Take a look at the article http://www.jstatsoft.org/v09/i01 and download the attached code. try to run one of the codes for example section 2.1 in R here is the code fossil <- read.table("fossil.dat",header=T) x <- fossil$age y <- 100000*fossil$strontium.ratio knots <- seq(94,121,length=25) n <- length(x) X <-

If you mean this paper by Fahrmeir: http://biomet.oxfordjournals.org/cgi/content/abstract/81/2/317 I would recommend BayesX: http://www.stat.uni-muenchen.de/~bayesx/. BayesX interfaces with R and estimates discrete (and continuous) time survival data with penalized regression methods. If you are looking for a bona fide Bayesian survival analysis method and do not wish to spend a lot of time

# plyr plyr is a set of tools for a common set of problems: you need to __split__ up a big data structure into homogeneous pieces, __apply__ a function to each piece and then __combine__ all the results back together. For example, you might want to: * fit the same model each patient subsets of a data frame * quickly calculate summary statistics for each group * perform group-wise

# plyr plyr is a set of tools for a common set of problems: you need to __split__ up a big data structure into homogeneous pieces, __apply__ a function to each piece and then __combine__ all the results back together. For example, you might want to: * fit the same model each patient subsets of a data frame * quickly calculate summary statistics for each group * perform group-wise

# plyr plyr is a set of tools for a common set of problems: you need to __split__ up a big data structure into homogeneous pieces, __apply__ a function to each piece and then __combine__ all the results back together. For example, you might want to: * fit the same model each patient subsets of a data frame * quickly calculate summary statistics for each group * perform group-wise

# plyr plyr is a set of tools for a common set of problems: you need to __split__ up a big data structure into homogeneous pieces, __apply__ a function to each piece and then __combine__ all the results back together. For example, you might want to: * fit the same model each patient subsets of a data frame * quickly calculate summary statistics for each group * perform group-wise

Dear useRs, version 2.2-2 of the "betareg" package has just been released on CRAN http://CRAN.R-project.org/package=betareg accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i02/ The package provides beta regression for data in the unit interval (0, 1) such as rates and proportions. The manuscript replicates several practical

Dear useRs, version 2.2-2 of the "betareg" package has just been released on CRAN http://CRAN.R-project.org/package=betareg accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i02/ The package provides beta regression for data in the unit interval (0, 1) such as rates and proportions. The manuscript replicates several practical

Hi all, I'm running samba-3.0.25b w/Linux 2.6.16.27 (SUSE 10.1), newly compiled. I've joined the AD domain (which is a subdomain) whose PDC is Win2K3 server and which is running in native mode but allowing NT4 members. From the HOWTO, I followed the instructions - manually created the computer account, then did the "net rpc join ...." and successfully joined the AD

=== actuar: An R Package for Actuarial Science === We are pleased to announce the immediate availability of version 1.0-0 of actuar. This release follows publication of our papers in JSS (*) and R News (**). From the NEWS file: Version 1.0-0 ============= NEW FEATURES o Improved support for regression credibility models. There is now an option to make the computations with the

=== actuar: An R Package for Actuarial Science === We are pleased to announce the immediate availability of version 1.0-0 of actuar. This release follows publication of our papers in JSS (*) and R News (**). From the NEWS file: Version 1.0-0 ============= NEW FEATURES o Improved support for regression credibility models. There is now an option to make the computations with the

Dear Achim, >>>>> Achim Zeileis <Achim.Zeileis at r-project.org> >>>>> on Fri, 10 Mar 2017 15:02:38 +0100 writes: > Hi, we came across the following unexpected (for us) > behavior in terms.formula: When determining whether a term > is duplicated, only the order of the arguments in function > calls seems to be checked but not

Hi, Using arima on data that are in a data frame, especially when adding xreg, would be much easier if the arima function contained 1) a "data=" argument 2) the possibility to include the covariate(s) in a formula style. Ideally the call could be something like > arima(symptome, order=c(1,0,0), xreg=~trait01*mesure0, data=anxiete) ( or arima(symptome~trait01*mesure0,

We are using res_fax_digium with a Sangoma PRI card on asterisk 1.8.13 The docs at http://docs.digium.com/FAX/fax_for_asterisk_admin_manual.pdf indicate v34 is supported, but when I enable it I get the message "res_fax_digium.c:1624 dgm_fax_new: V.34 not supported, will be ignored." Is v34 only supported with SpanDSP? Also, the res_fax.conf.sample does not indicate v34 as a valid

Hi, we came across the following unexpected (for us) behavior in terms.formula: When determining whether a term is duplicated, only the order of the arguments in function calls seems to be checked but not their names. Thus the terms f(x, a = z) and f(x, b = z) are deemed to be duplicated and one of the terms is thus dropped. R> attr(terms(y ~ f(x, a = z) + f(x, b = z)),