similar to: Still avoiding loops

I cannot operate lattice plots, 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 -----------------------------------------------------

Dear ALL, I fear my question has already been answered many times before, but I haven't fund that in archives... I am working on spatial datasets and, in most arrays I'm handling, there are two columns dedicated to (x,y)-coords. For different reasons - notably to draw image() plots, I need to have these two columns sorted in increasing order. But sort() and order() seem to apply to

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

Dear R-users, My question is how to get right F tests for random effects in random effect models (I hope this question has not been answered too many times yet - I didn't find an answer in rhelp archives). My data are in mca2 (enc.) : names(mca2) [1] "Lignee" "Pollinisateur" "Rendement" dim(mca2) [1] 100 3 replications(Rendement ~ Lignee *

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 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

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

Dear R-users, Is there an easy way to avoid points one upon another when ploting rows and columns 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

JGR has an fairly intelligent editor, and it works on Mac. Andy > From: Sean Davis > > Consider using ESS and xemacs or emacs. You get syntax highlighting, > auto-indent, command 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

Hi, I got the following error message when running a function of mine doing intensive computations: Erreur dans .Call("R_lazyLoadDBfetch", key, file, compressed, hook, PACKAGE = "base") : r?f?rence d'argument par d?faut r?cursive I haven't found neither where the problem lies nor what could be recursive. Besides, I wonder whether it might or not be a problem

Dear R-users, Sorry for reposting. I put it in another way : I want to test random effects in this random effect model : Rendement ~ Pollinisateur (random) + Lignee (random) + Pollinisateur:Lignee (random) Of course : summary(aov(Rendement ~ Pollinisateur * Lignee, data = mca2)) gives wrong tests for random effects. But : summary(aov1 <- aov(Rendement ~ Error(Pollinisateur * Lignee), data =

Dear R-users, Someone, who uses R under Mac, wants to insert a couple of small plots (each with several lines) in an article, but he has to reduce plots' size significantly. He did it (in pdf or enc. ps) but, unfortunately, everything is reduced but lines' width. Besides, 'lwd' argument in par() or plot() can't be less than one. I am not a Mac user and I don't know

Dear all, I tried to use layout argument in xyplot to get one panel per page. I have a dataframe named 'data' with the following variables: x, y = coords, sub, bloc = 2-level factors, etat = 5-level factor, I did : > lset(theme = col.whitebg()) > xyplot(y ~ x | bloc*sub , data=data, groups=etat, + layout=c(0,1,4), + main="Etat des plantes dans chaque

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 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

Dear friends, suppose i have two datasets: A and B A: id<-1:6 x<-c(1,2,3,4,5,6) y<-c(2,4,6,8,3,2) xy<-data.frame(id,x,y) B m<-c(1,1,3,3,5,5) n<-c(2,2,6,6,3,3) mn<-data.frame(m,n) Now, i want to perfomr two tasks: 1. get a subset of B,no duplicate values,: C: m n 1 2 3 6 5 3 2.Extract the values in A on the conditions that x=m and y=n the results should be: id x y 1 1 2 3

Hello, >I am hoping for some advice as to how I might create a subset of a >matrix. The matrix is 176 x 3530. The rows are individual records >and the columns words. I want to create a new matrix that only >consists of words which occur in at least 5 records. For example, >if column 7 is "charges" and this only appears in 4 records/rows >this variable would not

Here is an example of the data.frame that I have, df<-data.frame("case"=rep(1:5,each=9),"id"=rep(1:9,times=5),"x"=round(runif(length(rep(1:5,each=9))))) "case" represents the cases, "id" the persons, and "x" is the binary state. I would like to know in how many cases any two persons a. both have "1", b. the first has

Dear all, I'm not getting what I'm doing wrong. The line below from my read.fsa.bin function throws an error when just loading my AFLP package and disappears when I load the zoo package as well. #the line that throws the error Index <- which(Peak == rollmax(Peak, k = 1 + 2 * floor((min(diff(SizeStandard)) * Fs - 1) / 2), fill = -Inf)) #the error Error in UseMethod("rollmax")

Hello All, I have had considerable bad luck with attempting the following with for loops. Here is the problem: # Suppose we have a data.frame with the following data, which can be considered a type of edgelist (for those with networks backgrounds): # # V1 V2 # 1 A # 1 A # 1 B # 2 A # 3 C # 3 A # 3 C # 3 B # # I want the output of the function to produce a matrix, such that #each factor of