similar to: writing an R code for a given model

Dear all, I am looking for an R package that fits multivariate gaussian or non-gaussian longitudinal outcomes. I am especially interested to non-gaussian outcomes since the outcomes I've got are discrete (some are binomial and some are count data). Many thanks in advance, Abderrahim [[alternative HTML version deleted]]

Dear All, I wonder if there is a way to fit a generalized linear mixed models (for repeated binomial data) via a direct Maximum Likelihood Approach. The "glmm" in the "repeated" package (Lindsey), the "glmmPQL" in the "MASS" package (Ripley) and "glmmGIBBS" (Myle and Calyton) are not using the full maximum likelihood as I understand. The

Dear All, I have the following error message when I fitted lmer to a binary data with the "AGQ" option: Error in family$mu.eta(eta) : NAs are not allowed in subscripted assignments In addition: Warning message: IRLS iterations for PQL did not converge Any help? Thanks in advance, Abderrahim [[alternative HTML version deleted]]

Dear All, I wonder if there is an efficient way to fit the generalized linear mixed model for multivariate outcomes. More specifically, Suppose that for a given subject i and at a given time j we observe a multivariate outcome Yij = (Y_ij1, Y_ij2, ..., Y_ijK). where Y_ijk is a binomial(n_ijk, p_ijk). One way to jointly model the data is to use the following specification: g(p_ijk) =

I am trying to delete directories recursively in 'smbclient' like I typically do with 'rm -r' on a Unix shell, but I am just not able to. Is that in fact possible in 'smbclient'? Youssef Eldakar Bibliotheca Alexandrina

I'm building R-3.1.1 (64 bit) from source on AIX 7.1. It was going well until I hit this: xlc_r -q64 -Wl,-brtl -Wl,-G -Wl,-bexpall -Wl,-bnoentry -lc -L/opt/freeware/lib64 -L/opt/freeware/lib -Wl,-blibpath:/opt/freeware/lib64:/opt/freeware/lib:/usr/lib:/lib -Wl,-bmaxdata:0x80000000 -o tools.so text.o init.o Rmd5.o md5.o signals.o install.o getfmts.o http.o gramLatex.o gramRd.o -lm make[6]:

Thanks Martin to confirm that the issue is due to privilege access. How can I run a domain as non-root and be able to access the cpu information? -----Original Message----- From: Martin Kletzander [mailto:mkletzan@redhat.com] Sent: Tuesday, February 28, 2017 16:18 To: EL FATHI Youssef OBS/OINIS Cc: libvirt-users@redhat.com Subject: Re: [libvirt-users] Redhat 7: cgroup CPUACCT controller is not

Dear All, I am trying to graph a simple scatter plot where the x axis is year and the y axis is a percentage (percentage of infant death). Instead of plotting the raw data, I want to plot summary statistics such as mean and median. Here is the problem: the value range of y is between 0 and 1, but since infant death is a rare event, the mean and median is very low (something like 5%), which shows

Hello, I am trying to fit a nlme (non linear mixed effect). I am using the SelfStart function SSlogis. However the data in my hand contains few observations per subject (4 or less), so the nlsList doesn't work... In this case I should fixe initial parameter estimates. I remark that values of initial estimates have a greater effect on the model fit (i.e. loglikelihood, AIC and also on

Hi everyone, Does the selfstart function SSlogis of the "nlme" library allows the introduction of time varying covariates ? For example how can I interpret the xmid parameter (reperesenting the age at which we reach the half of the asymptote) if I want to explain it by a some time varying covariate? Thanks in adavance, Abderrahim Abderrahim Oulhaj, Phd in Statistics Oxford

I need to retrieve my asterisk to retrieve a values from postgresql, i am looking for some sort of application like *mysql*() app, I found one but it is only available on Suse, is there any way for doing this? Regards, O.Youssef -------------- next part -------------- An HTML attachment was scrubbed... URL:

Dear All, With kind help from several friends on the list, I am getting close. Now here are something interesting I just realized: for random effects, lmer reports standard deviation instead of standard error! Is there a hidden option that tells lmer to report standard error of random effects, like most other multilevel or mixed modeling software, so that we can say something like "randome

Dear All, I am fitting a simple GAM model using the "gam" function in the "mgcv" package. The only independent variable is a continuous variable representing the time of the event (in year and month) coded so that "0" represents January 1960, "1" represents February 1960, etc. Now when I try to plot the results, my x-axis ranges from "-200" to

Hi, With a non-root user account, I am launching virtual machines and would like to get CPU stats for each Core (using python API or not) but face the following problem: - When I issue the command "virsh --readonly cpu-stats MY_DOMAIN" I got the following error: error: Failed to retrieve CPU statistics for domain 'MY_DOMAIN' error: Requested operation is not valid: cgroup

I''ve seen a recent report claiming that it works: http://linux-bsd-sharing.blogspot.pt/2012/10/howto-xen-413-windows-8-hvm-domu-with.html However, in my case, with a slightly different HW configuration the story is totally different. :-（ My config: HW: I7-3770 + AsRock H77M-ITX board Dom0: Debian Wheezy stock xen (v4.1.3) + Debian stock kernel (3.2.x / 3.6.x) 64 bit version. Note that

Dear All, I have been using Sweave (mainly via the Sweave.sh script) and really like it. I am working a paper (using Sweave, of course) which includes several time-consuming computations, and it gets tedious to re-compile the whoel thing every time I made changes. Then I discover the "cacheSweave" package, which seems the right solution to my problem. I only have on problem. Here is

Dear All, I am trying to fiugre out how install packages on the base system. The "FAQ" tells me that more details can be obtained by typing "help(install.packages)", but my R system gives me an error message when I did that. I also tried "help(PACKAGES)", "help(install.packages())", but none worked. Anybody can help? I am running R 1.1.1 on WinNT 4.

Hi, I am looking for any procedure in R to fit "Latent variables models" in particular "Latent trait models (i.e. for binary data)". could anyone help me? Many thanks , Dr Abderrahim Oulhaj Oxford University Department of Pharmacology Mansfield Road Oxford OX1 3QT Tel: +44 1865 224098 Fax: +44 1865 224099 Email: abderrahim.oulhaj@pharmacology.oxford.ac.uk

Dear All, This is an absolutely beginner's question: which compiler should I use to compile R source code under Win NT platform? Thanks a lot. Shige Song Department of Sociology, UCLA -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

I was wondering whether somebody has developed a set of Textpad macro that work with R the same way as WinEdt does. I use Textpad for almost everything (programming, creating LaTex docs, etc.), it will be nice if I can use it as my R program editor. Or maybe someone can tell me how that kind of macro works so that I can develop my own macros. Thank you! Best, Shige Song Department of Sociology,