search for: qpcr

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/libra...

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.

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) cat("\n") cat("PRESS =",p...

...r as R has a threshold > of > 2x10^9 even in 64 bit R. > > It runs fine for the fungii variable. > > If you guys want to run the data (attached), the full command is below. > > Thanks. > > --------------------------------------------- > > ##Import data: > > qPCR <- read.delim(file.choose(), > header = TRUE, > dec = ".") > > ##Load package > > library(lme4) > > ##Other steps: > > qPCR$obs <- 1:nrow(qPCR) > qPCR$fID<-as.factor(qPCR$ID) > qPCR$fDiet<-as.factor(qPCR$Die...

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.

...le is not a an integer, and cannot be transformed into an integer as R has a threshold of 2x10^9 even in 64 bit R. It runs fine for the fungii variable. If you guys want to run the data (attached), the full command is below. Thanks. --------------------------------------------- ##Import data: qPCR <- read.delim(file.choose(), header = TRUE, dec = ".") ##Load package library(lme4) ##Other steps: qPCR$obs <- 1:nrow(qPCR) qPCR$fID<-as.factor(qPCR$ID) qPCR$fDiet<-as.factor(qPCR$Diet) ##Run the model: M1 <- glmer (Baci ~ fDiet + Cro...

I am using LOF.test() function from the qpcR package and got the following result: > LOF.test(nlregmod3) $pF [1] 0.97686 $pLR [1] 0.77025 Can I conclude from the LOF.test() results that my nonlinear regression model is significant/statistically significant? Where my nonlinear model was fitted as follows: nlregmod3 <- nlsr(formula=y...

...n answer: http://markmail.org/search/?q=list%3Aorg.r-project.r-help+PRESS#query:list%3Aorg.r-project.r-help%20PRESS+page:1+mid:k2mbz5sov5eo5ejw+state:results I also see that other packages have implemented PRESS at least as reported by others: Subject: [R] I need help computing PRESS statistics (qpcR package) of...: From: Francisco Goes (xico... at hotmail.com) Date: Jun 4, 2014 4:05:03 pm I tried a current search, although I admit that the fact that "press" is an acronym shared by other topics does seem to complicate that process. I counted 9 packages with PRESS functions even after...

.... For instance if the Main_Exp is 6, I would like to group the sub exp's for the main_exp and show that there were two sub_exp's done an ELISA and a FCM and that the ELISA was done once and the FCM was done twice. Similarly for Main_Exp 7 I would like to show that for there were ELISA,FCM,qPCR,Telometry done and that they had a count of 3,3,2,5 respectively. This is what my dataset looks like: Main_Exp_Name Sub_Exp_Name Sub_Exp_Count 1 ELISA 2 6 ELISA 1 6 FCM...

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

...R> source('bugs.R') library(minqa) library(rgl) newuoa(initpar, optimft) => OK <Start R> source('bugs.R') library(rgl) library(minqa) newuoa(initpar, optimft) => Crash: segfault: address 0x18, cause 'memory not mapped' I found the bug using the package qpcR, where rgl is loaded when loading qpcR while minqa is only loaded later, when needed. Running on Debian squeeze 64 bit. R version: R version 2.11.1 (2010-05-31) x86_64-pc-linux-gnu rgl version: 0.91 minqa version: 1.1.9 Rcpp version: 0.8.6 (loaded by minqa) Kind regards, Gaspard Lequeux

...r of rows, which I have to cut back then again. I couldn''t come up with any different way of writing the data without having R either give me an error for different dimensions or overwriting the unmatched NA''s and fill the whole matrix by repeating the row data. I am using the qpcR package and rbind.na. Is there a different way of achieving the same but more directly? I know it''s just one additional line, but it''s still bugging me and it will definitely come in handy at some time. Best regards, Karl ##Loading Packages library(qpcR) data =...

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 =

...thods base other attached packages: [1] lattice_0.17-25 loaded via a namespace (and not attached): [1] grid_2.9.1 tools_2.9.1 ##-------------------------------------------------------------- ## This dataset is too complicated, but it does show the type of plot I want. ## ## Create a fake qPCR dataset: Eight 96-well plates over 4 days (2 per day), ## 2 genes per plate (multiplexed), and 4 "Hi" positive control and ## 4 "Lo" positive controls per plate. ## Create the experimental data; by rights it is all identical, expect for ## experimental errors with in days and be...

Hi I get the following error when I try and get the AICc for a gls regression using qpcR: > AICc(gls1) Loading required package: nlme Error in n/(n - p - 1) : 'n' is missing My gls is like this: > gls1 Generalized least squares fit by REML Model: thercarnmax ~ therherbmax Data: NULL Log-restricted-likelihood: 2.328125 Coefficients: (Intercept) therherbmax 1.6...

.... For instance if the Main_Exp is 6, I would like to group the sub exp's for the main_exp and show that there were two sub_exp's done an ELISA and a FCM and that the ELISA was done once and the FCM was done twice. Similarly for Main_Exp 7 I would like to show that for there were ELISA,FCM,qPCR,Telometry done and that they had a count of 3,3,2,5 respectively. This is what my dataset looks like: Task Datatype Count 1 ELISA 2 6 ELISA 1 6 FCM...

I am currently trying to write a program that minimises the amount of work required for “auditable” qPCR data. At the moment I am using an Excel (.csv) spreadsheet as source data that has been transposed to the column format required for R to read. Unfortunately, this means I have* *to manually confirm the whole data set prior to doing any analysis, which is taking a considerable amount of time! My id...

Hello Bernardo, --------- If I understood your problem this script solve your problem: q<-0.15 + c(-.1,0,.1) h<-10 + c(-.1,0,.1) 5*q*h [1] 2.475 7.500 12.625 --------- OK, this solves the simple example. But what if the example is not that simple. E.g. P = 5 * q/h Here, to get the maximum tolerances for P, we need to divide the maximum value for q by the minimum value for h, and

Hi, I have several matrix in a list, for example: e [[1]] [,1] [,2] [1,] 1 3 [2,] 2 4 [[2]] [,1] [,2] [1,] 1 4 [2,] 2 5 [3,] 3 6 [[3]] [,1] [,2] [1,] 2 1 I would like to join them by column i.e. [,1] [,2] [,3] [,4][,5] [,6] [1,] 1 3 1 4 2 1 [2,] 2 4 2 5 NA NA [3,] NA NA 3 6 NA NA I have tried