similar to: Difficult with round() function

Hi, I am newbie of R. I a currently using multdrc object to generate fitting curve and IC50. My 384 well format raw data contains multi dose response curves. My script goes through set of data then produce curve and ic50. Here is my sudo code: For (plateid in platelist) { Input data (plateid) as matrix Curve fitting model4logistic <- multdrc(rdata ~ ld, logDose=10) }

Buenos dias R-help-es, Estoy interesado en estimar una curva dosis-respuesta para un conjunto de datos y para ello, estoy utilizando la libreria "drm". Hasta ahi todo bien. Me gustaria automatizar algunas cosas y el primer paso para ello es la estimacion del modelo. Si la estimacion funciona, todo lo demas funciona; de lo contrario, todo fallara. Tengo algunas lineas que mitigan un

I found this by accident yesterday. The program that crashes is the first two lines of the example from the help page for nlmer. That example hasn't changed in a long time, so I assumed that it is an R-devel issue. It could also be a long latent nlmer bug. The second run with valgrind is puzzling. Terry T. > library(lmer) > sessionInfo() R Under development (unstable)

Thanks, posted to https://github.com/lme4/lme4/issues/412 for further discussion ... On Wed, Feb 22, 2017 at 10:03 AM, Therneau, Terry M., Ph.D. <therneau at mayo.edu> wrote: > I found this by accident yesterday. The program that crashes is the first > two lines of the example from the help page for nlmer. That example hasn't > changed in a long time, so I assumed that it is

Hi, I'm trying to extract a histogram over the results from a bootstrap. However I keep receiving the error message "Error in hist.default(boot.lrtest$ll, breaks = "scott") : 'x' must be numeric". The bootstrap I'm running looks like: > boot.test <- function(data, indeces, maxit=20) { + y1 <- fit1+e1[indeces] + mod1 <- glm(y1 ~ X1-1, maxit=maxit) +

Hi all! I've read http://reactor-core.org/ogg-tagging.html and http://www.vorbis.com/ot/20030303.html#id2726753 about the goals and non-goals of vorbis comments, but I'm still unsure: 1) Is it by purpose to forget about CDA's index feature because it's so rarely used or did nobody here ever thought about it? -> Add an INDEXNUMER field? 2) If "Vorbis comments are the

Hello, I need to regress data like the example below. The data points represent friction factors derived from observed trip length data. The function used to describe this data is a gamma function of the form, f(t) = a * t^b * e^(c*t) and I need to regress the data to obtain the a,b, and c coefficients. The gamma function can also be expressed in the log-linear form, ln[f(t)] = ln[a] + b

Dear all, I am currently working on a function in which I would like to avoid using the command sample(). Therefore, I am now trying to make a for loop that does the same thing as the in built function sample does: Rearranging the items of a object randomly. So, the output I want to you get is the same as sample() would give me: e.g.: data <- c(5,4,6,7,8) sample(data) > data <-

Hi, I would like to know what is the easiest way to compile two or more set of vectors or data frame, according to their index. They are interrelated to one another by their assigned index. for example: #data set 1 abc #output: X403 X408 X410 X415 X418 X419 X420 X423 X424 X425 X426 X427 549.58 541.91 544.18 549.37 555.54 540.83 543.26 544.26 546.85 548.98 553.10 556.49

Full_Name: Isto Aho Version: 2.0.0 OS: Windows XP Submission from: (NULL) (80.186.9.143) Hi, The bug can be reproduced as follows. Install R 2.0.0 and RExcel. Take an excel sheet having 65000 rows and 11 columns, integer values. Give the 'context' command 'put R val' for the large enough selection of rows. (I didn't check, how many are needed, but it seems that 10000 rows

Hi, Jean-Marc Valin wrote: > How about this: the codebook search gives you the N best (in order) > combinations of entries for the current sub-frame. You can simply > compute a one-bit checksum on all the codebook entries and choose > which of the N-best to use based on that and the "message" you have. So you're suggesting, I can either write into the stream: -

Dear all, Are there any functions to convert decimals to fractions in R? I have the result: > summary(as.factor(complete.ID)) 0 0.0133333333333333 0.04 2256 488 230 0.0666666666666667 0.0933333333333333 0.106666666666667 2342 310 726 0.133333333333333

Hi every one, I am using the package 'drc' to model root elongation using dose response data. I don't know which model I should use. Though I don't know which model I should use, I tried the following codes given below. But it produced the error messages.Can any one tell me the code in 'drc' package to find out the EC (Effective Concentration) values and Confidence

Dear Contributors, I have an easy question for you which is puzzling me instead. I am running loops similar to the following: for (i in c(100,1000,10000)){ print((mean(i))) #var<-var(rnorm(i,0,1)) } This is what I obtain: [1] 100 [1] 1000 [1] 10000 In this case I ask the software to print out the result, but I would like to store it in an object. I have tried a second loop, because if I

Hi all, >From Monty's letter I know that ResVQTrain runs on vqd files that contain numbers separated by commas. I know that the .vqd files of this format generated by Vorbis compiled with RES_TRAIN contains partwords, i.e. codebook indeces assigned to partitions. I think that ResVQTrain should run on real data, not codebook indeces. Did I misunderstand something? Thanks Ervin

Hi I'm trying to develop some C code to find the fixpoint of a contraction mapping, the code compiles and gives the right results when executed in R. However R-gui session is frequently terminated. I'm suspecting some access violation error due to the exception code 0xc0000005 In the error report windows 10 gives me. It is the first time I'm writing any C-code so I'm guessing I

Hello R Gurus, This is a simple enough question, but I am curious as to whether there's an answer... Can R generate a random variable uniformly distributed on -Inf to Inf? Philosophically this doesn't seem possible, and if not, as I imagine so, is there some sort of generally accepted factor I should be multiplying by a Unif(-1,1) rv to sample from the real line? Thanks, Dan Polhamus

Hi all, I'm trying to check a package using R CMD check --as-cran <package>. When I do so, the check fails upon not finding the dependencies which I have installed in a non-standard location. Per the manual, I created ~/.R/check.Renviron : $ cat ~/.R/check.Renviron R_LIBS=/usr/lib64/R/library:/usr/lib64/R/var R_LIBS_SITE=/usr/lib64/R/library:/usr/lib64/R/var No luck. I tried

Hi, I'm having trouble understanding how to construct a random number generator for a parametric bootstrap. My aim is to bootstrap a Likelihood Ratio statistic (under the null) for a linear model. The function at this point is given by boot.test.n01 <- function(data, indeces, maxit=20) { y1 <- fit1+se(e2)*rnorm(314) mod1 <- glm(y1 ~ X1-1, maxit=maxit) y2 <-

Hello, I gave a list of 2 x 2 matrices called matlist. I have about 5000 2 x 2 matrices. I would like to count how many of each 2 x 2 unique matrix I have. So I am thinking that I need a list of the unique 2 x 2 matrices and their counts. Can anyone help. -- Thanks, Jim. [[alternative HTML version deleted]]