similar to: how to label the som notes by the majority vote

Dear Whom it may concern, I need help to figure the "macro" function in R: I need to plot the different data sets by a plotxyf function, I want the title to be different for different data set. # get the data set final.xyf<- xyf(data=as.matrix(my.final), Y=classvec2classmat(final$outcome), xweight = 0.5, grid=somgrid(5, 4, "hexagonal")) #plot function

I would really like some or any advice on how I can improve (or fix??) the following analysis. I hope I have provided a completely runnable code - it doesn't produce any errors for me. The resulting plot at the end shows a pretty poor correlation (just speaking visually here) to the test set. How can I improve the performance of the mapping and prediction? Here are some of the data

Hi all, I am reading a huge data set(12M rows) that contains family information, Offspring, Parent1 and Parent2 Parent1 and parent2 should be in the first column as an offspring before their offspring information. Their parent information (parent1 and parent2) should be set to zero, if unknown. Also the first column should be unique. Here is my sample data set and desired output. fam

Does anyone know of a method that I can get the intersection where the red and blue curves meet i.e. the value on the x-axis? x <- 1:10 y <- 10:1 plot(x,y) abline(lm(y~x),col="blue") abline(h=2.5,col="red") Muhammad

HI, Dear R community, My original file has 1932 lines, but when I read into R, it changed to 1068 lines, how comes? cdu@nuuk:~/operon$ wc -l id_name_gh5.txt 1932 id_name_gh5.txt > gene_name<-read.table("/home/cdu/operon/id_name_gh5.txt", sep="\t", skip=0, header=F, fill=T) > dim(gene_name) [1] 1068 3 -- Sincerely, Changbin -- Changbin Du DOE Joint Genome

Dear community, I have one file named ca_boost_feature.txt, Feature selection (Boosting:0.0025,5)! H.2.C C.1.D C.3.R E.0.N C.2.S C.0.G H.3.G log file: ep If I want to use the second line of this file, how to read it into R? varr<-read.table("/home/cdu/operon/carbonic/ca_boost_feature.txt", sep=" ", skip=1, header=F, strip.white=TRUE, nrows=1) Warning message: In

HI, Dear community, I am building the following loop, ww<-function(file) { lossw<-vector() for (x in seq(0.1, 0.9, by=0.1)) { cat('xweight ', x, '\n') lossw[i] <- cross.validation(file, x)$avg } return(lossw) } MY question is how to index the lossw[i]? for (i in 1:9) for (x in seq(0.1, 0.9, by=0.1)) Thanks so much!

HI, Dear R community, I have a large data set names dd.txt, the columns are: there are 2402 variables. a1, b1, ..z1, a11, b11, ...z11, a111, b111, ..z111.. IF I dont know the relative position of the columns, but I know I need the following variables: var<-c(a1, c1,a11,b11,f111) Can I use read.columns to read the data into R? I have tried the following codes, but it does not work

Hello all, I have two questions. One probably has a very simple answer but I have checked the FAQ, other websites and still have not found an answer. I am new to using R. My very simple question is how to do line breaks when creating an axis title (xlab) ? I have a few that are too long and they get cut off... I found on the internet the "/n" function but I think I must be using it

Dear list, i am using the following code to produce  a lattice xyplot, but the axis.break-function is seemingly not executed. Date<-seq(as.Date("2006-08-29"), as.Date("2007-08-28"), by="2 weeks") Period<- var1<-rnorm(27, 90000000, 30000000) var2<-rnorm(27, 500000000,250000000) var3<-rnorm(27, 1000000,500000) var4<-rnorm(27, 600000,300000)

Hello! Trying to find out some guides how to set up freebsd under xen, all i found is pointing to www.fsmware.com, <http://www.fsmware.com/xenofreebsd/7.0/download/> which is down! Maybe someone could share freebsd image and config? _______________________________________________ Xen-users mailing list Xen-users@lists.xensource.com http://lists.xensource.com/xen-users

HI, Dear R community, I have one data set like this, What I want to do is to calculate the cumulative coverage. The following codes works for small data set (#rows = 100), but when feed the whole data set, it still running after 24 hours. Can someone give some suggestions for long vector? id reads Contig79:1 4 Contig79:2 8 Contig79:3 13 Contig79:4 14 Contig79:5 17

From: Christopher Snow <cs at m0n0.co.uk> > If that's the case, shouldn't the credit vote go to Fedora? Disclaimer: This is 100% my _opinion_ (don't take it otherwise) 1. Credit At some point a lot of credit can go to a lot of projects. Yes, many are thanx to Red Hat -- by far the largest of any commercial entity. But then many are thanx to many maintainers in general.

Hi everybody, after the initial positive community review of the http://wiki.xen.org/wiki/Mirage_Incubation_Project_Proposal, it is time to have a formal vote. Who can vote? - Project leads of Mature Xen.org projects - Committers of Mature Xen.org projects - Others can register their vote in support of the project, but the vote is not binding. How to vote? - Vote via the voting form at

Hi everyone. I'm trying to break the y axis on a plot. For instance, I have 2 series (points and a loess). Since the loess is a "continuous" set of points, it passes in the break section. However, with gap.plot I cant plot the loess because of this (I got the message "some values of y will not be displayed"). Here's my code: library(plotrix); #generate some data x

Hi All, I am experimenting with using binary data (1, 0) in the kohonen package. I would like either 1 or 0 to be mapped in SOM, but the package gives me 1 or 0 (as expected), but something in between. I am not scaling the binary data, but still there are in-between values. Can anybody suggest a way to keep the data binary, even in the SOM plot? Thanks a lot. George

> first <- data.frame(a=1:3, b=4:6) > second <- data.frame(b=7:9, c=10:12) > third <- merge(first, second, by="b", all=TRUE) > third b a c 1 4 1 NA 2 5 2 NA 3 6 3 NA 4 7 NA 10 5 8 NA 11 6 9 NA 12 It's easy to replace the Nas with whatever value you want. No, merge() does not work with more than two data sets at a time. Andy From: mtb954 at gmail.com

In the name of Allah. hello, how are you? I sent a mail for you,yesterday but I didn't recieve to my aim. I need to some data frames in version 1.7.0 such as "florida","Barley","rainforest" and some others. I don't know these data frames are in which library. I typed "data(package = .packages(all.available = TRUE))" that was in "data()"

Me encuentro de vacaciones hasta el proximo 16/02. Por favor contactar a Pablo Minsteras (pablo.minsteras at utopixnetworks.com) o Gustavo Hernandez (gustavo.hernandez at utopixnetworks.com). Muchas Gracias, Gustavo Scheveloff

Hi all, I have a question about the SOM routine. You can either supply the initial representatives for the lattice yourself or else they are chosen randomly from the dataset. Is it possible to pass the random-seed as an argument somehow, when choosing the random initialisation of the lattice? As it is now, each time I run a SOM on a dataset with the same settings the resulting SOM will still