search for: veslot

...r avoding loops; for that, my first idea was to add this matrix: X1=X[,rep(1:46,46:1)] to this one: res=NULL for (i in (2:47)) res=c(res,i:47) X2=X[,res] (Is it a nice alternative way ?) Is there a way to create the second matrix X2 without a loop, such as for X1 ? Thanks in advance, Jacques VESLOT

...w image() plots, I need to have these two columns sorted in increasing order. But sort() and order() seem to apply to vectors - or one single column - only. Could somedy please let me know how to sort two columns of a dataframe, with priority to one of them, just like in access ? Thanks, Jacques VESLOT CIRAD Another related question : How to put legend on image() plots - or image-based plots like image.kriging() in geoR ? Despite help pages, I haven't succeeded. ----------------------------------------------------- CIRAD 3P Reunion - MailScanner -------------------------------------------...

...e), I can "see" standard deviations of random effects (but don't know how to find them) with : library(lme4) summary(lmer(Rendement ~ (1 |Pollinisateur) + (1 | Lignee) + (1 | Pollinisateur:Lignee), data=mca2)) but I can't get F tests. Thanks in advance. Best regards, Jacques VESLOT

...mand auto-complete, transcripts for your session, > integrated help, and the tools of one of the most powerful > text editors > on the planet. The learning curve is a bit steep, but if you use R > much, it is worth it. > > Sean > > On Jan 21, 2005, at 5:35 AM, Jacques VESLOT wrote: > > > Dear all, > > > > Could someone please make me know if there is a nice script editor > > available > > under Mac, similar to Crimson, that offers R syntax > highlighting (and > > pairs > > of parentheses underlining) ? > > >...

Hi, I would like to be able to use R with a connection to a PostgreSQL database. I am using R 2.0.1 on windows XP. I have tried installing both the DBI and rDBI packages (which is better?) but in either case I run into the same problem - when i try to use either the dbDriver or dbConnect functions i'm told the function couldn't be found. Seems like the packages aren't installing

Hi all, I have a very large binary vector, I wish to calculate the number of 1's over sliding windows. this is my very slow function slide<-function(seq,window){ n<-length(seq)-window tot<-c() tot[1]<-sum(seq[1:window]) for (i in 2:n) { tot[i]<- tot[i-1]-seq[i-1]+seq[i] } return(tot) } this works well for for reasonably sized vectors. Does

...ts, notably those in package gstat, within for() loops. (I can neither read such graphics included in a .R (or .txt) file with the source() function). However same commands executed one by one don't fail. Could you please help me find a means to execute lattice plots inside a loop ? Jacques VESLOT CIRAD ----------------------------------------------------- CIRAD 3P Reunion - MailScanner ----------------------------------------------------- Ce message a ??t?? v??rifi?? par MailScanner et analyse par 4 antivirus : ClamAv - BitDefender - F-Prot - Panda. Aucun virus,ni rien de suspect n'...

Hi, Use Newdata <- subset(d2004, select=-c(concentration,stade)) See ?subset for details Best, Matthias > -----Urspr??ngliche Nachricht----- > Von: Jacques VESLOT [mailto:jacques.veslot at cirad.fr] > Gesendet: Donnerstag, 09. September 2004 12:38 > An: r-help at stat.math.ethz.ch > Betreff: [R] Indexing dataframe > > > I am sorry to ask such question, but I can't find a solution... > > I have a dataframe 'd2004' and I...

My friend Veslot proposed me the Raven ! http://www.teteamodeler.com/allopass/images/corbeau.jpg As a matter of fact the intelligence of this bird is comparable to one of a monkey, although its brain is close to a reptile brain. R is quite the same : small, compact, but so clever. This is a canadian website ex...

...plot() can't be less than one. I am not a Mac user and I don't know whether it is a Mac-related problem or not. Could somebody please give me a way to solve this problem, either to reduce line width under R, or to get reduced lines' width when reducing plot size ? Thanks, Jacques VESLOT Cirad

...olumns of 'dudi' objects ? Maybe there is a function in ade4 or in an other package, or maybe someone has his or her own function to do this (for example to automatically modify a little the coordinates of these points to get a readable plot ?). Thanks in advance. Best regards, Jacques VESLOT

...the problem lies nor what could be recursive. Besides, I wonder whether it might or not be a problem of memory size. Could you please give me any suggestion on how to interpret this message? Thanks in advance, jacques ------------------------------------------------------------------- Jacques VESLOT CNRS UMR 8090 I.B.L (2?me ?tage) 1 rue du Professeur Calmette B.P. 245 59019 Lille Cedex Tel : 33 (0)3.20.87.10.44 Fax : 33 (0)3.20.87.10.31 http://www-good.ibl.fr

Dear all, I have the following data frame: > raw.count Var1 Freq 1 AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 707 2 AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAC 14 3 AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAT 3 But why when printint it using 'cat', it doesn't print the desired string "AAA" ? > cat(raw.count$Var1,

Dear All, I have a data, suppose it is an N*M matrix data. All I want is to classify it into, let see, 3 classes. Which method(s) do you think is(are) appropriate for this purpose? Any reference will be welcome! Thanks! Best, Baoqiang Cao

...51), c(-1, 15))) and received this error message : Error in if (!any(cond.max.level - cond.current.level < 0) && (row - 1) * : missing value where TRUE/FALSE needed I tried some changes in arguments - notably layout=c(0,1), but anything works. Thanks for helping... Jacques VESLOT CIRAD R??union

Hi R guru coders I wrote a bit of code to add a new column onto a "topTable" dataframe. That is a list of genes processed using the limma package. I used a for loop but I kept feeling there was a better way using a more vector oriented approach. I looked at several commands such as "apply", "by" etc but could not find a good way to do it. I have this feeling there

...ot;see" standard deviations of random effects (I don't know how to find them) but I can't find F tests for random effects. I only want to know if there is an easy way (a function ?) to do F tests for random effects in random effect models. Thanks in advance. Best regards, Jacques VESLOT Data and output are as follows : > head(mca2) Lignee Pollinisateur Rendement 1 L1 P1 13.4 2 L1 P1 13.3 3 L2 P1 12.4 4 L2 P1 12.6 5 L3 P1 12.7 6 L3 P1 13.0 > names(...

Dear all, I have a dataset. I want to make barplot from this data. Zero1 <- " V1 V2 V3 V4 V5 V6 V7 V8 V9 1 1 0 0 0 1 0 0 0 Positive 2 0 0 1 0 1 0 1 1 Negative 3 0 0 1 0 0 0 1 1 Positive 4 0 1 0 1 1 1 0 1 Negative 5 0 0 1 0 1 1 0 0 Positive 6 0 1 0 0 1 1 1 1 Negative 7 1 0 1 1 1 1 1 1 Negative 8 0 0 0 0 1 0 0 1

...'d like to use column names (to avoid finding column numbers each time). I cannot find an easy way to operate... I wonder why that works: > d2004[, "concentration"] and this don't: > d2004 <- d2004[, -c("concentration", "stade")] Thanks... Jacques VESLOT / CIRAD

...Test data: pheno by groups 1 < 2 < 3 chi-squared = 0.2268, df = 1, p-value = 0.6339 > tabletrend(with(dat2, table(pheno, geno))) [1] 0.6338308 -------- Message original -------- Sujet: Trend test and test for homogeneity of odd-ratios Date: Wed, 16 Aug 2006 17:39:33 +0200 De: Jacques VESLOT <jacques.veslot at good.ibl.fr> Pour: R-Help <r-help at stat.math.ethz.ch> Dear r-users, I am looking for some R functions to do Cochran-Armitage trend test for 2*3 tables (binary phenotype vs. genotypes) and for testing the homogeneity of odds ratios within 2*3*k tables (binary pheno...