similar to: Manually set UCL and LCL

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

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;

Hello again I'm trying to put error bars onto a barplot. I've tried something that Bill Simpson suggested a while ago, ie: x<-c(1,2,3,4,5) y<-c(1.1, 2.3, 3.0, 3.9, 5.1) ucl<-c(1.3, 2.4, 3.5, 4.1, 5.3) lcl<-c(.9, 1.8, 2.7, 3.8, 5.0) plot(x,y, ylim=range(c(lcl,ucl))) arrows(x,ucl,x,lcl,length=.05,angle=90,code=3) #or segments(x,ucl,x,lcl) but I can't get it to work on a

hello, please can anyone help me out. Am a new user of R program. Am having problem with this code below, not getting the expected results. 1. Each m, the cumulative sum should be 1.000 but the 2nd and 3rd m returned 2.000 and 3.000 instead of 1.000. 2. to get the LCL(m) and UCL(m) for each m base on these instructions if out.cum > 0.025 then LCL(m)= y-1 if out.cum >0.975

Folks I have modified an existing function to calculate 'ec/ld/lc' 50 values and their associated Fieller's confidence limits. It is based on EC50.calc (writtien by John Bailer) - but also borrows from the dose.p (MASS) function. My goal was to make the original EC50.calc function flexible with respect to 1) probability at which to calculate the expected dose, and 2) the link

hello, please can anyone help me out. Am a new user of R program. Am having problem with this code below, not getting the expected results. 1. Each m, the cumulative sum should be 1.000 but the 2nd and 3rd m returned 2.000 and 3.000 instead of 1.000. 2. to get the LCL(m) and UCL(m) for each m base on these instructions if out.cum > 0.025 then LCL(m)= y-1 if out.cum

I'm new to Hmisc and trying to get the following to work, but if I un- comment the y-scale list (in order to get a log-scale for the hazard ratio), the error bars become strangely large. The dataframe is simply ODS output from TPHREG in SAS. Can someone point me towards what I'm sure is a naive error? Many thanks! tphreg1_parameterestimates1987$LCL <-

Deadline - 6th February 2015 Five PhD awards to begin September 2015, (one including an additional MRes award in Urban Sustainability and Resilience in collaboration with UCL?s Centre for Urban Sustainability and Resilience), are available for work associated with the Economic and Social Research Council (ESRC) funded Consumer Data Research Centre, based at University College London (UCL) and the

The Child and Adolescent Mental Health Services (CAMHS) Evidence Based Practice Unit (EBPU) is a dynamic, expanding and friendly academic and service development unit, currently consisting of 19 people including: researchers, clinicians and service development specialists (http://www.ucl.ac.uk/clinical-psychology/EBPU/). It is part of University College London (UCL) and the Anna Freud Centre

Hello R-Users, The following questions are not R-technical, but more of general statistical nature. 1. NORMALITY I built a normal linear regression model and now I want to check for the residual normality assumption. If I check the distribution graphically and look at the descriptive characteristics (skewness and kurtosis are below 1), I would confirm that the residuals are normally

I have windows version of R v1.6.2 I have downloaded the qtoolbox.zip from http://www.cmis.csiro.au/S-PLUS/qtoolbox/ My R command line is > install.packages("C:/Program Files/R/qtoolbox.zip", .libPaths()[1], CRAN = NULL) and I get the following error message: updating HTML package descriptions Warning message: error -1 in extracting from zip file > Is there another way to

Hello R users, I would like to create several forestplots with the same X axis, so, if you were to look at the plots lined up all the X axes would be identical (and the different plots could be compared). Here is one version of code I've used: mytk10<-c(0.1, 0.5, 1, 2, 5, 10) pdf(file = "myfile.pdf", pointsize = 7, paper="letter", width=6, height=9)

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

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) *

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

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

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 =