search for: leony

Dear List I'm a newbie R user. I'm utilizing the list function in order to make a var like this: clusters<-list(a=var1,b=var2....) My problem is that the total numer of variables that I need to include in my list is up to 200. I've the text string with the complete list of my variables but is too long to cut and paste in my bash shell. So is there a way too import the list from a

Dear list I've a generic question about how to tune an SVM I'm trying to classify with caret package some population data from a case-control study . In each column of my matrix there are the SNP genotypes , in each row there are the individuals. I correctly splitted my total dataset in training(132 individuals) and test (50 individuals) (respecting the total observed genotypic

Dear list I'm testing a predictor and I produced nice performance plots with ROCR package utilizing the 3 standard command pred <- prediction(predictions, labels) perf <- performance(pred, measure = "tpr", x.measure = "fpr") plot(perf, col=rainbow(10)) The pred object and the perfo object are S4 with the following slots An object of class "performance"

Dear All, I am planning to run a Monte-Carlo experiment which involves to do roughly 100.000 times the following 1- Generating a sample of, say, 50.000 numbers from an ARMA or GARCH 2- Doing some regressions on the series 3- On each regression storing one special value from the results into a sequence that will be analyzed later. The experiment is calculation-intensive, and I suspect some

Hi - I am doing some monte carlo simulations comparing bayesian (using Plummer's jags) and maximum likelihood (using lmer from package lme4 by Bates et al). I would like to know if there is a way I can flag nonconvergence and exceptions. Currently the simulations just stop and the output reads things like: Error in optim(.Call("lmer_coef", x, 2, PACKAGE = "Matrix"), fn,

Dear list I wish to extract from a population genotypized for 10 SNP a subsample of the same population of size n with similar allele frequencies. Essentially i have a matrix of 200 rows (df) like this Name,Condition,rs1385699_X,rs6625163_X,rs962458_X,Rs4658627_1, sample01,Case,1,1,1,-1 sample02,Control,1,1,1,1 sample06,Control,1,-1,1,0 sample10,Case,1,1,1,0 sample11,Control,1,1,1,1

Hola a todos, En la empresa donde trabajo estamos interesados en contratar un curso ad-hoc de R en nuestras oficinas de Madrid. El nivel del curso requerido es de iniciación a intermedio, en modalidad intensiva (2 a 3 días completos) y de ser posible utilizando la herramienta R Studio. En esta ocasión, las búsquedas en Google no están siendo exitosas y solamente hemos encontrado empresas /

Hi - I am doing a monte carlo experiment that requires to do a linear regression of a matrix of vectors of dependent variables on a fixed set of covariates (one regression per vector). I am wondering if anyone has any idea of how to speed up the computations in R. The code follows: #regression function #Linear regression code qreg <- function(y,x) { X=cbind(1,x) m<-lm.fit(y=y,x=X)

Hi - I need to export data to openoffice base, where one of the elements is a long character vector (>255 characters.) write.dbf exports it as varchar, truncating the data. Any idea how to do this? thanks, -eduardo

Hi all, How could I set a timer in R, so that at fixed interval, the R program will invoke some other functions to run some tasks? Thank you very much!

Hello, I have a bit of a tricky puzzle with trying to implement a logit model as described in a paper. The particular paper is on horseracing and they explain a model that is a logit trained "per race", yet somehow the coefficients are combined across all the training races to come up with a final set of coefficients. My understanding is that they maximize log likelihood across the

Dear user is the first time that I use wine. I'm trying to install with wine 1.31 a technical software for analyzing some biological data. Starting from a shell I retrieve always the same error Please could someone give to me some suggestion about what this error mean? thank you very much Guido fixme:sfc:SFC_3 0 fixme:advapi:SaferCloseLevel ((nil)) stub err:ole:create_server class

Dear list Is there a way to extract a random sample without duplicated row from a dataframe ?. a=c(1,2,3,1,1,1,2,1) b=c(1,2,3,1,2,1,2,1) c=c(1,1,1,1,1,1,1,1) d=c(1,2,3,1,1,1,2,1) prov<-data.frame(a,b,c,d) prov2<-prov[sample(1:nrow(prov),5,replace=T),] prov2 a b c d 3 3 3 1 3 6 1 1 1 1 3.1 3 3 1 3 5 1 2 1 1 8 1 1 1 1 I tryed the above code but as you can see sample function

Dear list I'd like to have your opinion about my case study. I'm analizing a dataset of 9 experiments and 15 variables with the aim to highlight the variables that can majorly explain the variance between the experiments. This is an example with only 3 rows and 5 variables var1 var2 var3 var4 var5 sample5 0,067 0,005 0,008 0,100 0,005 sample6 0,069 0,001

Same behavior was noted on R 2.9.0 OS X and WIndows. > t(sapply(seq(0.01, .99, length=20), function(x) c(x, length(agrep("Staatssekretar im Bundeskanzleramt","Bundeskanzler",max.distance=x))))) [1,] 0.01000000 0 [2,] 0.06157895 0 [3,] 0.11315789 0 [4,] 0.16473684 0 [5,] 0.21631579 0 [6,] 0.26789474 0 [7,] 0.31947368 0 [8,] 0.37105263 0 [9,]

Dear List I'm trying to install a package not present in cran named A2R ( http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=79) After running the demo script I retrieve the following error: cannot change value of locked binding for '._a2r_counter' Please could someone give to me a tip about my error? Thank you very much Here is my sessionInfo() R version 2.15.0

I write about R every weekday at the Revolutions blog: http://blog.revolution-computing.com . In case you missed them, here are some articles from last month of particular interest to R users. http://bit.ly/un680 demonstrated reader Paul Bleicher's code for visualizing a time series as a heat-map calendar. http://bit.ly/5fgis0 and http://bit.ly/mepBH showed (with thanks to Drew Conway) how

Dear all, I try to show a subset of coefficients in my presentation. It seems that a "standard" table is not a good way to go. I found figure 9 (page 9) in this file ( http://www.destatis.de/jetspeed/portal/cms/Sites/destatis/Internet/DE/Content/Wissenschaftsforum/Kolloquien/VisualisierungModellierung__Beitrag,property=file.pdf ) looks pretty good. I wonder if there is any function for

Hi, I'm trying to get data from web page and modify it in R. I have a problem with encoding. I'm not able to get encoding right in htmlTreeParse command. See below > library(RCurl) > library(XML) > > site <- getURL("http://www.aarresaari.net/jobboard/jobs.html") > txt <- readLines(tc <- textConnection(site)); close(tc) > txt <- htmlTreeParse(txt,