similar to: NLME

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, 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, 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) =

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]]

Hi All: This may turn out to be very simply, but I can't seem to add the name of the school to a chart. The loop I created is below that subsets a dataframe and creates a chart for each school based on certain variables. As it stands now, they title includes the school's ID number. Instead, I want to replace this with the school's actual name, which is stored in a variable called

Hi All, I finally signed up for the list after years of using Samba successfully - a testament to the quality of Samba. Yet now I have a problem with the point-and-print functionality. I am able to authenticate against my server (Solaris 8, Samba 3.0.7, OpenLDAP 2.1.25) as user 'chuck' in my LDAP directory and browse the shares, but when I right-click on the printer and select

Ok.. Just for the record, I STAND CORRECTED! The problem does appear to be NT.. The thing that irritated me was that nobody was interested in giving any backup data! Everyone just wanted to blindly admit that it must be NT. Ok.. Enuf bitching.. On to the better stuff... Anyway, David Mansfield and I are trying to narrow down the machine configurations where this happens.. Currently, we've

Dear R help list, I am desperately looking for any reference explaining by examples how to write R codes in order to fit the parameters of a given model using maximum likelihood or any other criteria function. I know the general structure: First write a code for the maximum likelihood function and afterwards write a code to maximize it using optim and then invert the Hessian to get the

Ernie, Thanks for the great info.. Based on what you were seeing, would you suspect that part of the problem is with the combination of NT and the 3com card (3c595)? Currently, we are ONLY seeing this under NT Workstation (not server) and almost all of our machines have this same 3com 3c595 card (only if using 100BaseT).. -- Rick -----Original Message----- From: Ernie Oporto

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]]

Hi all, I am trying to interparete the result of the following output from lm; fit1 =lm(Feed _Intake ~ weight + season + weight*season) Season has three classes(x,y,z) Reults are Estimate (Intercept) 21.51559 weight 2.13051 factor(season)y 10.59739 factor(season)z 1.30421 weight:factor(season)y

I have a system with many processes writing to a common data file using pwrite. These processes are each writing to existing blocks in the file, not changing the file size, and the file has no holes. When the processes get going, they seem to bottleneck at log_wait_common (according to ps alnx). That is, one process is uninterruptible in log_wait_common, the rest are uninterruptible in down.

Dear R-experts, I am running a script that has the following structure: windows(height=5.5,width=8) dat<-read.csv("myfile.csv") names(dat)<-c('a','b','c','d') dat<-dat[,1:4] xyplot(dat$a~dat$b) Then I usually save the plot as a PDF (from the menu in the R console). I can save the PDF twice in a row; but at the third time, right when I

Dear R helpers, I am working with R 2.4.1 GUI 1.18 (4038) for MacOSX. I have a matrix of 10 000 genes and try to run the following commands: > model.mix<-makeModel (data=data, formula=~Dye+Array+Sample+Time, random=~Array+Sample) > anova.mix<-fitmaanova (data, model.mix) > test.mix<-matest (data, model=model.mix, term="Time", n.perm=100, test.method=c(1,0,1,1))

Hi, I was having question about mysql driver support ( not odbc). Do we still need the asterisk-add-on to be installed for mysql support.? If yes, Which version should be used and from where I should get it? Thanks in adavance. ---- Thanks & Regards, PrashantAbhang -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hello, What is the recommended way to offer high availability to Guests? For example, let host1, host2 and host3 be libvirt and KVM/Qemu enabled. If I start a guest in host1, how can I guarantee that it will stay online if host1 goes down? GlusterFS, for example, would take care of storage; but what about CPU and RAM? How can this be accomplished? Thank you for any feedback or comment in

An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20040112/e9dceeb4/attachment.htm -------------- next part -------------- Hello One of the option in VoicemailMain() is "Adavance Options". Could anyone explain what are these ?. Because whenever I select Advance Options, it repeats the same process of asking "Change Folders,Advance

I have been regularly crashing my NT boxen trying to back them up via an smbtar (smbclient). The crash is not complete - the network layer on NT goes away. You can still operate the console, but no net traffic. This is NT workstation and server 4.0 with sp3. The UNIX samba server is Solaris 2.5 with the latest samba version (1.9.18p8?) also with other samba versions. I tried updating the NT

Hello everybody, I would like to know whether it is possible to run Dial and MeetMe commands simultaneoously on the same channel. I am using a C AGI as below but it seems to me that only the first command that is called in the agi is executed. ........... // Pr?paration de la commande pour l'appel du client fprintf(stderr,"%s%s",numtocall," is the number to