search for: emax

Dear list, I am running some linear regressions through lapply, >lapply(c('EMAX','EC50','KOUT','GAMMA'),function(x)confint(lm(get(x)~RR0,dataset2))) I got results like [[1]] 2.5 % 97.5 % (Intercept) 0.6595789212 0.8821691261 RR0 -0.0001801771 0.0001489083 [[2]] 2.5 % 97.5 % (Intercept) -63...

Colleagues, I am fitting an Emax model using nls. The code is: START <- list(EMAX=INITEMAX, EFFECT=INITEFFECT, C50=INITC50) CONTROL <- list(maxiter=1000, warnOnly=T) #FORMULA <- as.formula(YVAR ~ EMAX - EFFECT * XVAR^GAMMA / (XVAR^GAMMA + C50^GAMMA)) ## alternate version of formula FORMULA <- as.formula(YVAR...

...G 1 = PK, FLAG 2=PD, ;CMT 1=GUT, 2=CENTRAL ;DV =PK ORDERED BY ID, TIME $DATA ../projectRSD1.csv IGNORE=C $SUBROUTINE ADVAN2 TRANS2 INFN=../MItoolsRunlogNM6.for $ABB COMRES=8 $PK FLA2=FLAG TVKA=THETA(1) TVCL=THETA(2) TVV=THETA(3) CL=TVCL*EXP(ETA(1)) V=TVV*EXP(ETA(2)) KA=TVKA*EXP(ETA(3)) S2=V $ERROR EMAX=THETA(4)*EXP(ETA(4)) EC50=THETA(5)*EXP(ETA(5)) N = THETA(6)*EXP(ETA(6)) EFF= (EMAX*F**N)/(EC50**N+F**N) + ERR(2) CONC= F + F*ERR(1) Y = CONC*(2-FLAG) + EFF*(FLAG-1) IPRED=Y " LAST " COM(1)=G(1,1) " COM(2)=G(2,1) " COM(3)=G(3,1) " COM(4)=G(4,1) &...

...rs quite stable. In the process of moving to new servers I noticed beta 4 and compiled it for FreeBSD 10 on a dell poweredge. I have been unable to get the administrative interface working properly. You can view screenshots showing the problems and a copy of the error log at http://vm3.emax.ca/icecast-bug/ Regards, -- Barry Friedman Emax Computer Systems Inc., 480 Tweedsmuir Ave., Ottawa, Ont. Canada K1Z 5N9 bfriedman at emax.ca 613-725-3198 Mobile: 613-261-0767

Dear R Users: I have a list function as: Flat: y = 0 Linear: y = -(1.65/8)d Logistic: y = 0.015 - 1.73/{1+exp[1.2(4-d)]} Umbrella: y= -(1.65/3)d + (1.65/36)d^2 Emax: y = -1.81d/(0.79+d) Sigmoid Emax: y = -1.70d^5/(4^5+d^5) And want draw the figure as attached (those material are extracted from a paper). Could anyone give me a sample code to do this? Thanks

...nlme object "mod.nlme". I have two questions 1.How do I extract the degrees of freedom? 2.How is this degrees of freedom calculated in an nlme model? Thanks. Jun Shen Some sample code and data ================================================================= mod.nlme<-nlme(RESP~E0+(Emax-E0)*CP**gamma/(EC50**gamma+CP**gamma),data=Data, fixed=E0+Emax+gamma+EC50~1, random=list(pdDiag(EC50+E0+gamma~1)), groups=~ID, start=list(fixed=c(E0=1,Emax=100,gamma=1,b=50)) ) The "Data" object structure(list(ID = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L,...

...timism index.corrected n Dxy 0.6425 0.7017 0.6131 0.0887 0.5539 1000 R2 0.3270 0.3716 0.3335 0.0382 0.2888 1000 Intercept 0.0000 0.0000 0.0821 -0.0821 0.0821 1000 Slope 1.0000 1.0000 1.0548 -0.0548 1.0548 1000 Emax 0.0000 0.0000 0.0263 0.0263 0.0263 1000 > validate(full.approx.lrm, bw=F, B=1000) index.orig training test optimism index.corrected n Dxy 0.6446 0.6891 0.6265 0.0626 0.5820 1000 R2 0.3245 0.3592 0.3428 0.0164...

Hi there, I am attempting to analyse data from a crossover study using a mixed effect Emax model in R using nlme. I am having problems ascertaining how to fit the model. Subjects are randomly assigned to 1 of 4 treatments in a complete block design. The response variable of interest is recorded once within each period. So, subjects have a response for each period. Specifically, I wan...

I have tried multiple different methods to figure out how to get a date axis of my preference (start date of each month). Any assistance would be appreciated. The below section is not producing a date axis: plot(totaldays$totaldays,totaldays$y,type="n",ylim=c(0,Emax),xaxt="n") *xlabels<-(strptime(totaldays$totaldays,"%Y-%m-%d",tz="")) xlabels<-xlabels[xlabels$mday==1] axis.Date(1,at=xlabels,format="%b-%Y")* The above section is used in the following plot: par(mfrow=c(4,1)) par(mar=c(0.8,0,0,0)) par(oma=c(2,4.5,3,...

Hello I am fitting a Cox PH model using the function coxph(). Does anyone know how to obtain the estimate of the covariance matrix under the null hypothesis. The function coxph.detail() does not seem to be useful for this purpose. Thanks, KD. [[alternative HTML version deleted]]

Hello, I have this data: Time AMP 0 0.2000000 10 0.1958350 20 0.2914560 40 0.6763628 60 0.8494534 90 0.9874526 120 1.0477692 where AMP is the concentration of this metabolite with time. If you plot the data, you can see that it could be fitted using a logistic regression. For this purpose, I used this code: AMP.nls <- nls(AMP~SSlogis(Time,Asym, xmid, scal), data

...test optimism index.corrected n Dxy 0.9932 0.9940 0.9905 0.0035 0.9897 363 R2 0.9291 0.9364 0.9163 0.0202 0.9089 363 Intercept 0.0000 0.0000 0.0233 -0.0233 0.0233 363 Slope 1.0000 1.0000 0.7836 0.2164 0.7836 363 Emax 0.0000 0.0000 0.0582 0.0582 0.0582 363 D 0.9118 0.9190 0.8915 0.0275 0.8844 363 U -0.0110 -0.0110 0.0124 -0.0234 0.0124 363 Q 0.9228 0.9299 0.8791 0.0508 0.8720 363 B 0.0205 0.0172 0.02...

...measure" for each one. That way I can pick the "best" model from all my experiments. As mentioned above, I'm mainly interested in the accuracy of my probability predictions. Does anyone have an opinion about which measure I should look at?? (I see Dxy, C, R2, D, U, Briar, Emax, Eavg, etc.) Thanks!! -N

This is to announce the release of Icecast 2.4 beta5 (aka 2.3.99.5). We're looking forward to testing and feedback. As this is a beta release please discuss it on the icecast-dev mailing list. (Please do NOT cross-post to both lists, do NOT reply to the icecast mailing list, PLEASE reply to icecast-dev at xiph.org only). * Updated web interface to be fully XHTML compliant and mobile

This is to announce the release of Icecast 2.4 beta5 (aka 2.3.99.5). We're looking forward to testing and feedback. As this is a beta release please discuss it on the icecast-dev mailing list. (Please do NOT cross-post to both lists, do NOT reply to the icecast mailing list, PLEASE reply to icecast-dev at xiph.org only). * Updated web interface to be fully XHTML compliant and mobile

...ndex.corrected n Dxy 0.3861 0.4081 0.3699 0.0382 0.3479 200 R2 0.1537 0.1716 0.1378 0.0339 0.1198 200 Intercept 0.0000 0.0000 -0.0585 0.0585 -0.0585 200 Slope 1.0000 1.0000 0.8835 0.1165 0.8835 200 Emax 0.0000 0.0000 0.0375 0.0375 0.0375 200 D 0.1160 0.1315 0.1030 0.0285 0.0875 200 U -0.0063 -0.0063 0.0021 -0.0084 0.0021 200 Q 0.1223 0.1378 0.1010 0.0369 0.0855 200 B 0.2073 0.2035...

[This email is either empty or too large to be displayed at this time]