similar to: add labelled contour lines to filled.contour plot

Hi Alejo, According to my knowledge the two plots are different because in the first one a point belongs to a group depending on its group in the data whereas in the second plot a point belongs to the group predicted by the linear discriminant analysis. I hope somebody will correct me if I am wrong. Alain Alejo C.S. wrote: > Hi Alain, this is the code: > > > library(MASS) >

I'm trying to assign NAs to values that satisfy certain conditions (more complex than shown below) and it gives the right result, but breaks the loop having done the first one viz: new<-c(rep(5,4),6) for (i in 1:6) {new[new[i]>5.5][i]<-NA} gives the correct result, though an error message appears which causes a break if it's in a loop. If I can get rid of the error message and

Hello, I am using the recode() function in Rcmdr and the result is not what I expect so I am almost sure I did something wrong but what... > test <- data.frame(x=1:10) > library(car) > recode(test$x,'1:5=0 ; else=1', as.factor.result=TRUE) [1] 0 0 0 0 0 1 1 1 1 1 Levels: 0 1 BUT > library(Rcmdr) # recode from the car package is now masked Now I recode test$x

Dear R-users, I am using R(a package igraph) to calculate certain topological features of networks. When I try to draw a plot between these features I get an error. Following is the code I am using : *> plot(met_eco_deg,met_eco_bet) > lmout<-lm(met_eco_bet ~ met_eco_deg) > abline(lmout) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new

Dear all, i'm trying to label specific individuals (supplementary ones) after a PCA with the FactoMiner package. There is not much details (possibilities?) in the R-help of the plot.pca function. There is indeed a "label" parameter but i could only manage to label the supplementary individuals with there "row.names" (i.e. label="indiv.sup") and not with the

Hello, I have another stupid question. I hope you can give me a hint how to solve this: I have a list and one element is again a list containing matrices, all of the same dimensions. Now, I'd like to set the dimnames for all matrices: example code: m1 <- matrix(1:25, nrow=5) m2 <- matrix(26:50, nrow=5) # ... there can be much more than two matrices l <- list() l[[1]] <-

Dear Forum, I have a class variable 1 (populations A-E), and two other class variables, variable 2 and variable 3. What I want is to see if the combination of var 2 and var 3, will give me a pattern that allows to distinguish populations. I found several packages like ade4, with pcaiv function and factoMineR. but there are not working. Using the ade4 package, when I try to build the pca: pca1

I'm confused on how is the right way to plot a discriminant analysis made by lda function (MASS package). (I had attached my data fro reproduction). When I plot a lda object : X <- read.table("data", header=T) lda_analysis <- lda(formula(X), data=X) plot(lda_analysis) #the above plot is completely different to: plot(predict(lda_analysis)$x,

Hi all, I hope it's not too trivial for the list - I'm trying to concatenate two factor arrays, and obtain the following: > f1<-factor(c("a","a","b")) > f1 [1] a a b Levels: a b > f2<-factor(c("b","b","a")) > f2 [1] b b a Levels: a b > c(f1,f2) [1] 1 1 2 2 2 1 Instead of getting: [1] a a b b b a Levels: a

I am using the order function and the result seems to be incorrect: > a<-c(20,30,15,40) > order(a) [1] 3 1 2 4 Any suggestions? Thanks, Laura -- View this message in context: http://www.nabble.com/Incorrect-order-tp20224993p20224993.html Sent from the R help mailing list archive at Nabble.com.

Dear fellow users Are there any Vietnamese language resources for beginners of R? If so, I would be interested in hearing from people who have had experience with them and which are better (if there is more than one). I am involved with an aid project in Vietnam, and would like to move the scientists involved from using Excel for 'analysis' to R. Thanks .... Peter Alspach The

Dear list, I have a list of three matrices : i = list(matrix(1:4,2,2), matrix(3:6,2,2), matrix(9:12,2,2)) I would like to sum the matrices, as follows : [,1] [,2] [1,] 13 19 [2,] 16 22 I used this code : k <- i[[1]] for (j in (2:length(i))) { k <- k + i[[j]]} But, is it possible to sum without a loop ? Thanks in advance, Carlos [[alternative HTML version deleted]]

Hi everyone, I 'm building a function, in the middle it controls the sign of a variable x. If x < 0 the function write a warning (Error: negative value!). At this point I want the function stops without execute the remaining code. How can I do to terminate the function before your ending? Thanks in advance. Paolo

Hi, I am trying to create a function for a goodness-of-fit test for the Pareto Distribution for some loss data that I have. So far I have the following: function(X=OTOL) { n <- length(X)-1 #calculated the number of values (extra as 0 included) i <- 2:640 #values of i j <- 1:639 #values of i-1 Y <- (n-j+1)*((X[i])-(X[j])) #First part of GoF model Y } Where OTOL is the ordered loss

Dear all, let suppose I have following vector:   > dat1 <- c(rep("asd", 5), rep("xyz", 12), rep("erd", 17)) > dat1 <- dat1[sample(1:length(dat1), length(dat1), replace=F)] > dat1  [1] "erd" "xyz" "erd" "asd" "asd" "erd" "xyz" "asd" "erd" "erd"

Dear R users, When I use aggregate with table as FUN, I get what I would call a strange behaviour if it involves numerical vectors and one "level" of it is not present for every "levels" of the "by" variable: --------------------------- > df <- data.frame(A=c(1,1,1,1,0,0,0,0),B=c(1,0,1,0,0,0,1,0),C=c(1,0,1,0,0,1,1,1)) >

Hi, Does tapply change the order when applied on a factor? Below is the code I tried. > mylevels<-c("IN0020020155","IN0019800021","IN0020020064") >

Hello everyone. When I execute the following in R > (18-46)/(45-93) [1] 0.5833333 I get small precision for what I am trying to deal with . Is it possible to increase the precision for this and for other operations? For example openoffice calc for this operation returns 0.58333333333333300000 I I would like to thank you for your help [[alternative HTML version

Hey listmembers, I am desperately trying to write a data frame to a file. Not in CSV, but as they "appear" on the screen (nice, easy to read tables). I've read that the sink function is the way to go. I have tried the following code inside a function. sink("ABC.txt") MyFrameA MyFrameB sink() It gives the result I need when I use it outside a function, but when I use it

Hi, I created a principal component plot using the first two principal components. I used the function princomp() to calculate the scores. now, I would like to superimpose an ellipse representing the center and the 95% confidence interval of a series of points in my plot (as to illustrate the grouping of my samples). I looked at the ellipse() function in the ellipse package but can't get it