similar to: Calculating with tolerances (error propagation)

Dear list, I am from an engineering background, accustomed to work with tolerances. For example, I have measured Q = 0.15 +- 0.01 m^3/s H = 10 +- 0.1 m and now I want to calculate P = 5 * Q * H and get a value with a tolerance +- What is the elegant way of doing this in R? Thank you, Jan

I am using the MASS library function fitdistr(x, dpois, list(lambda=2)) but I get Error in optim(start, mylogfn, x = x, hessian = TRUE, ...) : Function cannot be evaluated at initial parameters In addition: There were 50 or more warnings (use warnings() to see the first 50) and all the first 50 warnings say 1: non-integer x = 1.452222 etc Can anyone tell me what I am doing

I want to use function cbind.na at library(qpcR) I install package qpcR and I can use functions such m1 <- pcrfit(reps, 1, 2, l5) > AICc(m1) [1] -102.5843 but when i try cbind.na(1, 1:7) i take message Error: could not find function "cbind.na" Thanks -- View this message in context: http://r.789695.n4.nabble.com/library-qpcR-cbind-na-tp4023339p4023339.html Sent from the

I am using R on a Windows XP professional platform. The following code is part of a bigger one CODE press=function(y,x){ library(qpcR) models.press=numeric(0) cat("\n") dep=y print(dep) indep=log(x) print(indep) yfit=dep-PRESS(lm(dep~indep))[[2]] cat("\n yfit\n") print(yfit) yfit.orig=yfit presid=y-yfit.orig press=sum(presid^2)

Ok. I`m trying to run a Poisson glmm with an observation-level random intercept. But I`m getting the following error for the 'Baci' variable: 'Error: (maxstephalfit) PIRLS step-halvings failed to reduce deviance in pwrssUpdate'. I guess this message is because the baci variable is not a an integer, and cannot be transformed into an integer as R has a threshold of 2x10^9 even in

Dear all, Is there anyway too generate MA plot for 2 qPCR assays (an array of 2x 400). -- View this message in context: http://r.789695.n4.nabble.com/Bioconductor-MA-plot-for-qPCR-array-tp4182805p4182805.html Sent from the R help mailing list archive at Nabble.com.

Inline. -- Bert On Tue, Apr 26, 2016 at 2:25 PM, Andr? Luis Neves <andrluis at ualberta.ca> wrote: > Ok. I`m trying to run a Poisson glmm with an observation-level random > intercept. But I`m getting the following error for the 'Baci' variable: > > 'Error: (maxstephalfit) PIRLS step-halvings failed to reduce deviance in > pwrssUpdate'. I guess this message

> On Apr 10, 2016, at 9:38 AM, David Winsemius <dwinsemius at comcast.net> wrote: > >> >> On Apr 10, 2016, at 3:11 AM, Murray Efford <murray.efford at otago.ac.nz> wrote: >> >> Martin - >> Thanks, but although hatvalues() is useful for calculating PRESS, I can't find anything directly relevant to my question in the influence help pages. After

Hi! Is There a way to manually create an affybatch object from qPCR array data? -- View this message in context: http://r.789695.n4.nabble.com/Creating-affybatch-objects-from-matrix-data-from-qPCR-array-tp3921559p3921559.html Sent from the R help mailing list archive at Nabble.com.

Can you explain why you need them as 'integer', A floating point representation can hold a value upto ~4.5e15 as an "integer" keeping the precision that you might need. Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it. On Tue, Apr 26, 2016 at 1:11 PM, Andr? Luis Neves <andrluis at

> On Apr 10, 2016, at 3:11 AM, Murray Efford <murray.efford at otago.ac.nz> wrote: > > Martin - > Thanks, but although hatvalues() is useful for calculating PRESS, I can't find anything directly relevant to my question in the influence help pages. After some burrowing in the literature I'm doubting there is an answer out there (PRESS R^2 is always presented in a fairly

Hi, Many apologies for sending this twice. I accidentally hit the send button before I finished writing my mail. I am new to R and I hope someone can help me with my problem. I am trying to draw a side by side barplot. There is a main experiment and there are many sub experiments within the main experiment. I would like to draw a bar plot showing the number and type of sub_experiments done for

Dear Helpers   I have a sample frame and i have sampled from it using three methods and now i want to calculate the statistics but i only get the population parameters.   H <- matrix(rnorm(100, mean=50000, sd=5000)) sampleframe=data.frame(type=c(rep("H",100)),value=c(H)) sampleframe   str=strata(sampleframe,c("type"),size=c(20,), method="srswor")

Hej, Calling newuoa (from the minqa package) makes R crash when the package rgl is loaded first. This however only on certain selected data. The data used for testing (saved to 'bugs.R'): xvals = c(1,2,4,5,7,8,9,10,11,12,14,15,16,18,19,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36) yvals =

I have the clinical study data. Year 0 Year 3 Retinol (nmol/L) N Mean +-sd Mean +-sd Vitamin A group 73 1.89+-0.36 2.06+-0.53 Trace group 57 1.83+-0.31 1.78+-0.30 where N is the number of male for the clinical study. I want to test if the mean serum retinol has increased over 3 years among subjects in the vitamin A group. > 1.89+0.36

Dear all: I converted the columns (Baci, Meti, Fungii, Protozoai) into integers (using excel) and then imported the data (.txt) into R. Interestingly, the other three variables were loaded as INT, but the 'Baci' one continued as Num. I imported the data using the following command line: X <- read.delim(file.choose(), header = TRUE, dec =

Hi folks, I am planning to write some more time-consuming matrix manipulations in c++. What is the experience with the existing c++ matrix libraries? Do you have some recommendations? Are some libraries more compatible with R than the others? All suggestions welcome! Best, Ott

Hi, The paper @ARTICLE{AlmironSilvaMM:2009, author = {Almiron, M. and Almeida, E. S. and Miranda, M.}, title = {The Reliability of Statistical Functions in Four Software Packages Freely used in Numerical Computation}, journal = {Brazilian Journal of Probability and Statistics}, year = {in press}, volume = {Special Issue on Statistical Image and Signal Processing}, url =

I have tightened the tolerances in glm.control in R-devel (aka 1.8.0 Under Development) from epsilon = 1e-4 to 1e-8, and increases maxit from 10 to 25. Normally the effect is to do one more iteration and get more accurate results. However, in cases of partial separation several more iterations will be done and it will be clearer from the results which parameters are converging and which are

Hi, I am looking to build even quintiles for a set of data. I managed to get it done, but I would like to know if there is a more direct way to write the data from my loop output x in the bottom of the code into the "empty" matrix p1, which I filled with NA''s. The way I am doing it at the moment is, more or less adding the matrix x after p1 and then deleting in a second step