similar to: attach

Sorry, I pasted the wrong file in earlier... this is the correct one: pValCalculator(b, n=20, m=20) { ind <- 1:min(c(n,m)) prob <- (1-pnorm(b,sd=std*sqrt(ind))) prob1 <- sum((n-ind+1)*(m-ind+1)*prob) prob1 } inputData <- scan("/users/lvssso/projects/LAMA/output/pValLamaScore.tmp", list(block1 = "",block2 = "",width1 = 0,width2 = 0,alignment = 0,score

Hi all, I'm currently working on the fitdistrplus package (that basically fit distributions). There is something I do not understand about the generic function summary. In the current version on CRAN, there is no NAMESPACE saying S3method(summary, fitdist) . However if we use summary on an object send by fitdist function it works fine... According to R-lang, we have " The most

Hi all, I tried to a nice legend with small boxes filled with the colors used for the plots. But it does nor work, boxes are always filled with black. An example is here plot(1:4,1:4) lines(1:4,4:1, col="blue") legend("top",leg=c("a","b"),col=c("black","blue"), fill=TRUE) How could I specify the colors? the argument col.box is the

Hi, I have a simple question that I just cannot figure out. I have 2 corresponding columns of data, one column (X-axis) for time (formatted thus: 8:30:01am = 830.1, 12:30:05pm = 1230.5, and one column (Y-axis) for values. When I attempt to plot the data using something like plot(inputdata[,1],inputdata[,2],type="l"); I get breaks in the plot (since the time essentially jumps from

Hello, I have created a .dll file using G77 and MinGW on my PC (Windows 2000). After using dyn.load to bring it into R2.0.1, I then call the .dll through the function ccprox shown below. It returns the correct values. If I run it a second time though it returns different values, so it seems something is being placed oddly in memory. If I unload and reload the .dll it works again the first

Dear list, Don't you think it could be useful to have in R base a function rowCumSums, that compute cumulative sums for each row of a matrix? My implementation of rowCumSums is rowCumSums <- function(x) t(mapply(function(row)cumsum(x[row,]), 1:NROW(x))) I'm sure it can be improved to have other arguments like na.rm or dims. Is there any hope to have this function in R?

Hi all, I use the symbol functions to draw circles. But the argument lty does not work on windows but works correctly on my macbook: x <- -4:4 y <- -4:4 plot(x,y,type="n") symbols(0, 0, add = TRUE, circles = 1, inches = 1, fg="black", lty="solid") symbols(0, 0, add = TRUE, circles = 1, inches = 2, fg="black", lty="dashed") The

Dear useRs, I have just discovered that the R optim function does not return the number of iterations. I still wonder why line 632-634 of optim C, the iter variable is not returned (for the BFGS method for example) ? Is there any trick to compute the iteration number with function call number? Kind regards Christophe -- Christophe Dutang Ph.D. student at ISFA, Lyon, France website:

Hi, I have a problem when checking the package 'probdistr' (on probability distributions). I got this warning * checking index information ... WARNING Vignettes with missing or empty \VignetteIndexEntry: [1] "probdistr-chi" "probdistr-contextra" "probdistr-discrete" [4] "probdistr-discrextra" "probdistr-exp"

Hi everybody, I a, trying to use R as a pipe like this: cat inputData | R --silent RCommandFile >outputData The RCommandFile would contain something like readLines(stdin()). I have tryed various things and none did work cleanly. One possible solution is to use the pipe() function inside R and to pass in the "cat inputData" however this is not very convenient since I would like to

Dear all, I have just found a 'good' tutorial R for datamining. I think it should be on the contributed docs. http://cran.r-project.org/other-docs.html Here is the link http://www.liaad.up.pt/~ltorgo/DataMiningWithR/ What do you think? Kind regards Christophe -- Christophe Dutang Ph. D. student at ISFA, Lyon, France website: http://dutangc.free.fr

Hi useRs, I'm trying to fit a basic AR(2) model with the 'ar' function. And when I try to check the value of the coefficients, I could not find the same value as the 'ar' function. Here is my example: myserie <- c(212, 205, 210, 213, 217, 222, 216, 218, 220, 212, 215, 236) #plot(myserie, type="l") myserieminus0 <- tail(myserie, -2) myserieminus1 <-

Hi useRs, I'm trying to fit a basic AR(2) model with the 'ar' function. And when I try to check the value of the coefficients, I could not find the same value as the 'ar' function. Here is my example: myserie <- c(212, 205, 210, 213, 217, 222, 216, 218, 220, 212, 215, 236) #plot(myserie, type="l") myserieminus0 <- tail(myserie, -2) myserieminus1 <-

Dear useRs, We are happy to announce a new version of the randtoolbox* package (version 1.09) on CRAN. The package is dedicated to Random Number Generation and in a less exhaustive way to RNG tests. The package is source-controlled at r-forge as part of the Rmetrics project http://r-forge.r-project.org/projects/rmetrics/ . We invite you to take a look at the NEWS section to know the

Dear useRs, We are happy to announce a new version of the randtoolbox* package (version 1.09) on CRAN. The package is dedicated to Random Number Generation and in a less exhaustive way to RNG tests. The package is source-controlled at r-forge as part of the Rmetrics project http://r-forge.r-project.org/projects/rmetrics/ . We invite you to take a look at the NEWS section to know the

Hello, Writing to seek help in regard to some unexpected performance anomaly i am observing in using tsoutlers:tso on the mac vs on an AWS cloud server.. I am running the following code with very small dataset of about 208 records. d.dir <- '/Users/darshanpandya/xxxxxx' FNAME <- 'my_data.csv' d.input <- fread(file.path(paste0(d.dir,"/zzz/"),FNAME,fsep =

Hello everyone, Hopefully a simple question. How can I access function argument alignment information? If my code has something like this: void foo(short * __restrict __attribute__ ((aligned(8))) InputData){} ...as I deal with argument (InputData) can I get the alignment info set for it? Thanks! Sergei --- Qualcomm Innovation Center, Inc. is a member of Code Aurora Forum, hosted by

Greetings, My goal is to create a Markov transition matrix (probability of moving from one state to another) with the 'highest traffic' portion of the matrix occupying the top-left section. Consider the following sample: inputData <- c( c(5, 3, 1, 6, 7), c(9, 7, 3, 10, 11), c(1, 2, 3, 4, 5), c(2, 4, 6, 8, 10), c(9, 5, 2, 1, 1) ) MAT <- matrix(inputData,

Hello, Writing to seek help in regard to some unexpected performance anomaly i am observing in using tsoutlers:tso on the mac vs on an AWS cloud server.. I am running the following code with very small dataset of about 208 records. d.dir <- '/Users/darshanpandya/xxxxxx' FNAME <- 'my_data.csv' d.input <- fread(file.path(paste0(d.dir,"/zzz/"),FNAME,fsep =

If you mean this paper by Fahrmeir: http://biomet.oxfordjournals.org/cgi/content/abstract/81/2/317 I would recommend BayesX: http://www.stat.uni-muenchen.de/~bayesx/. BayesX interfaces with R and estimates discrete (and continuous) time survival data with penalized regression methods. If you are looking for a bona fide Bayesian survival analysis method and do not wish to spend a lot of time