similar to: Pharmacokinetics using R

library(nlme) data(Phenobarb) na.include <- function(x)x phe1 <- nlme(conc~phenoModel(Subject, time, dose, lCl, lV), data = Phenobarb, fixed = lCl+lV~1, random= pdDiag(lCl+lV~1), start = c(-5,0), na.action = na.include, naPattern = ~!is.na(conc)) phe.ranef <- ranef(phe1,augFrame=TRUE) plot(phe.ranef, form=lCl~Wt+ApgarInd) [Error in max(length(x0),

Hi, In <- as.numeric(strsplit("11111011111101111111111111111111", "")[[1]]) Sets <- rep(letters[1:4], each = 8) ## your sets Sp.In <- split(In, Sets) logical <- sapply(Sp.In, function(x) any(x == 0)) c(as.matrix(data.frame(Sp.In))* rep(!logical, each = 8)) ## for equal sets This could give you a hint Thomas -----Urspr?ngliche Nachricht----- Von: Pravin

This concerns the "Clinical Study of Quinidine" example on page 380 of the book "Mixed-Effects Models in S and S-PLUS" by Pinheiro and Bates (2000). I have tried to reproduce the example, but get an error: > library(nlme) > fm1Quin.nlme <- nlme(conc ~ quinModel(Subject, time, conc, dose, interval, lV, lKa, lCl), + data=Quinidine, +

Hello, I am a beginner in R programming and recently heard about this mailing list. Currently, I am trapped into a simple problem for which I just can't find a solution. I have a huge dataset (~81,000 observations) that has been analyzed and the final result is in the form of 0 and 1(one column). I need to write a code to process this column in a little complicated way. These 81,000

Hello all, I am conducting a randomization test on a given dataset. One of the covariates, gender, is randomly assigned 1000 times to create a randomization dataset (rand.data). To these 1000 datasets, I fit a full model (see below) and the aim is to generate a distribution of LRT statistic under the null. Here are some of the questions I have: 1. I am using for() loop to fit the dataset

I cannot add users to my server and after going through many google searches and trying several things here is an issue I have found at minimum a LDAP issue when comparing my primary and secondary domain controllers. I do not know how to repair it and that is why I am contacting the list. Below is the output of the comparison. Master: 4.2.0pre1-GIT-58865d9 Slave: 4.2.0pre1-GIT-815bde2 Thank

Hello all. I'm having a problem debugging an IVR I'm building. I can't see any reason this shouldn't be working. Firstly the asterisk version is: Asterisk SVN-trunk-r7230 built by root @ localhost.localdomain on a i686 running Linux on 2006-02-17 22:44:48 UTC Basically the problem is this. While the playbacks are happening you can push any one of the options and to happily

I'm working on an IVR that gives the users the option (number 5 in the main menu) to dial by extension: exten => 5,1,Set(TIMEOUT(digit)=5) ; Dial Extension exten => 5,2,Set(TIMEOUT(response)=10) exten => 5,3,Background(LCL/prompt-60) exten => 5,4,WaitExten(15) When going option 5 you can dial some extensions such as 2802, it goes to the extension (all extens start with 28 on the

A couple of weeks ago there was a question regarding apparent convergence in nls when using the SSfol selfStart model for fitting a first-order pharmacokinetic model. I can't manage to find the original message either in my archive or in the list archives but the data were time conc dose 0.50 5.40 1 0.75 11.10 1 1.00 8.40 1 1.25 13.80 1 1.50 15.50 1

Dear listers, I am trying to fit a model using nlsList() using alternately a SSfol() selfstart function or its developped equivalent formulae. This preliminary trial works well mydata<-groupedData(Conc~Tps|Organ,data=mydata) mymod1<-nls(Conc~SSfol(Dose,Tps,lKe,lKa,lCl),data=mydata) as well as a developped form: mymod2<-nls(Conc~Dose * exp(lKe+lKa-lCl) *

Hi, I wonder if anybody could help me in converting this easy SAS program into R. (I'm still trying to do that!) PROC IMPORT OUT= WORK.CHLA_italian DATAFILE= "C:\Documents and Settings\carleal\My Documents\REBECCA\stat\sas\All&nutrients.xls" DBMS=EXCEL2000 REPLACE; GETNAMES=YES; RUN; data chla_italian; set chla_italian;

Dear friends of nlme, Running quinModel (Pinheiro/Bates page 380) on R (current release, windows) gives: Nonlinear mixed-effects model fit by maximum likelihood Model: conc ~ quinModel(Subject, time, conc, dose, interval, lV, lKa, lCl) Data: Quinidine Log-likelihood: -497 Fixed: lV + lKa + lCl ~ 1 lV lKa lCl 5.382 -0.273 2.470 Random effects: Formula: list(lV ~ 1, lCl ~ 1)

Dear all, I've been struggling for days now with a piece of code that I have posted here before, that has a really obscure bug. I think I may have isolated it, but I have no idea what it is.... It might also be a bug in R I guess, as it seems that one or several of list elements are not passed when a function is called, but quite rarely. I have been hacking rather wildly on the histogram

Using samba 3.05 I'm trying to upload printer drivers for point and print using an xp client and I'm getting an error. Unable to install driver ...... Operation could not be completed. I've followed the Samba How To [global] workgroup = local server string = CH File/Print Server passdb backend = ldapsam:ldap://ldap.mycompany.lcl log file =

I have moved to 4.1.16 on both controllers. It did not solve my issue (though it reduced the errors as seen in the new attached ldap comparison). I am still getting the error telling me to see the administrator when trying to add a user via the Windows Active Directory Users and Computers utility. Gilbert [root at pdc ~]# /usr/local/samba/bin/samba-tool ldapcmp ldap://pdc ldap://sdc domain

On 22/01/15 00:11, Gilbert Gutierrez wrote: > I have moved to 4.1.16 on both controllers. It did not solve my issue > (though it reduced the errors as seen in the new attached ldap > comparison). I am still getting the error telling me to see the > administrator when trying to add a user via the Windows Active > Directory Users and Computers utility. > > Gilbert > >

Pravin, You don't necessarily need a login script, but it is nice for mounting shared drives and printers, or just copying/installing files. Use an empty login script to start, then use the Windows "command line" (the "net use /?" command is helpful here) to mount and attach to services you want. Put those commands into the login script and use environment

<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> <head> <title></title> </head> <body> <a class="moz-txt-link-abbreviated" href="mailto:speex-dev@xiph.org">speex-dev@xiph.org</a> wrote: <blockquote type="cite" cite="mid1050398492.11049.TMDA@motherfish-II.xiph.org"> <pre

Hi Gents, I need to get the control limits from qcc function. As follows: qcc(MDI, type = "xbar.one") Call: qcc(data = MDI, type = "xbar.one") xbar.one chart for MDI Summary of group statistics: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.3266 0.4249 0.4371 0.4333 0.4451 0.4858 Group sample size: 1 Number of groups: 383 Center

On Wed, 1 Mar 2017 16:24:36 +0100 Tony Peña <emperor.cu at gmail.com> wrote: > Hi thanks by answer quickly > > yes is a ad pdc, i refers to PDC i think will be the same, now i saw > isn't. anyway, this is the smb.conf and krb5.conf I suggest you change your smb.conf to: [global] workgroup = sambadc realm = SAMBADC.LCL netbios name = samba-dc server string