similar to: hist plot

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 All, I tried to plot a variety of lines(curves) on same figure. What I did is, plot(x=x1,y=y1) lines(x=x2,y=y2) lines(x=x3,y=y3) ... In my data, the maximum of y1 is much smaller than those maximums of other y vectors. So, in the figure I got, there are some curves which are not complete, I mean, they were cut off at the maximum of y1 at the y axis. Could anybody point out some right

Dear All, I am producing a figure with many curves on it. How do I make the legends for all those curves smaller so that it can fit the figure itself? The commands I used for ploting are: plot(x1,y1,col=1,lty=1) lines(x2,y2,col=2,lty=2) ... legend(0.3,0.4,c("name1","name2",...),col=1:20,lty=1:20) Any tips for making the legend fit the figure will very welcome! Thanks! Best,

Dear All, I am using wilcox.test to test two samples, data_a and data_b, earch sample has 3 replicates, suppose data_a and data_b are 20*3 matrix. Then I used the following to test the null hypothesis (they are from same distribution.): wilcox.test(x=data_a, y=data_b, alternative="g") I got pvalue = 1.90806170863311e-09. When I switched data_a and data_b by doing the following:

Dear All, I'm learning to train a neural network with my training data by using nnet package, then evaluate it with a evaluation set. My problem here is that, I need the trained network to be used in future, so, what should I store? and How? Any other options other than nnet package? Any example will be highly appreciated! Best, Baoqiang Cao

Hi Everyone, I'm trying to use cutree to get the clusters after hclust. What I used is: mycluster<-cutree(cnclust,h=0.5) Now, my problem is, how can I get the actual clusters? Thanks! Best, Baoqiang Cao

Dear All, Just start to use the long expected R, my focus will be doing clustering on microarray data, just wonder, anyone can show me any references to conquer the steep learning curve? Thanks! Best regards, Baoqiang Cao

Dear All, I was trying to format a numeric vector (100*1) by using outd <- format(x=m, sci=F, digits=2) > outd[1:10] [1] " 0.01787758" "-0.14760306" "-0.45806041" "-0.67858525" "-0.64591748" [6] "-0.05918100" "-0.25632276" "-0.15980138" "-0.08359873" "-0.37866688" >m[1:10] [1]

Dear All, I have such a two-class problem, one class is very large(~98% of total), and the other is just 2%. According to manual of nnet, I need setup "weights", so I intend to set 1 for class one, 49 for class 2. How do I do that? Just weights=49? Meanwhile I'd like to try svm(e1071), again, how do I setup "class.weights"? Thanks. BTW: Many thanks to Jake and Uwe for

Hi, I'd like to download some data files from a remote server, the problem here is that some of the files actually don't exist, which I don't know before try. Just wondering if a function in R could tell me if a file exists on a remote server? I searched this mailing list and after read severals mails, still clueless. Any help will be highly appreciated. B.C.

Dear All, I'm running the latest R on WinXP by using Rgui.exe --max-mem-size=30000Mb . After read in a huge data file to a data matrix, I tried to get a subset of the data matrix but failed with: Error: cannot allocate vector of size 12079 Kb Any tips to get out of it? If helps, > memory.limit() [1] 3145728000 > memory.size() [1] 842735624 Thanks! Best, Cao

Dear All, After followed the help of nnet, I could get the networks trained and, excitedly, get the prediction for other samples. It is a two classes data set, I used "N" and "P" to label the two. My question is, how do I get the predicted numerical value for each sample? Not just give me the label(either "N" or "P")? Thanks! FYI: The nnet example I

Hi, I have a question about using lm on matrix, have to admit it is very trivial but I just couldn't find the answer after searched the mailing list and other online tutorial. It would be great if you could help. I have a matrix "trainx" of 492(rows) by 220(columns) that is my x, and trainy is 492 by 1. Also, I have the newdata testx which is 240 (rows) by 220 (columns). Here is

Hi, I'm trying to plot something in the following way and would like if you could help: I'd like in a same plot window, two plots are shown, the left one is a bird-view plot of the whole data, the right half keep changing, i.e., different plots will be shown up on request, so that when I select/click on some where in the left plot, the right plot will be the corresponding plot. What I

Hi, I have the following two measurements stored in mat: > print(mat) [,1] [,2] [1,] -14.80976 -265.786 [2,] -14.92417 -54.724 [3,] -13.92087 -58.912 [4,] -9.11503 -115.580 [5,] -17.05970 -278.749 [6,] -25.23313 -219.513 [7,] -19.62465 -497.873 [8,] -13.92087 -659.486 [9,] -14.24629 -131.680 [10,] -20.81758 -604.961 [11,] -15.32194 -18.735 To calculate the ranking

Dear All, I'd like to know if there is anyway to download a certain structure file from http://www.pdb.org/. I tried the following but failed: tmp <- download.file(url="ftp://ftp.wwpdb.org/pub/pdb/data/structures/all/pdb", destfile="pdb100d.ent.gz") trying URL 'ftp://ftp.wwpdb.org/pub/pdb/data/structures/all/pdb' Error in download.file(url =

---- Original message ---- >Date: Thu, 14 Jul 2005 19:50:41 +0200 >From: Uwe Ligges <ligges at statistik.uni-dortmund.de> >Subject: Re: [R] East Asian language >To: Nan Lin <nlin at math.wustl.edu> >Cc: r-help at stat.math.ethz.ch > >Nan Lin wrote: > >> Dear all, >> >> I just installed R 2.1.1. The installation program automatically

Hi all, I'm trying to superimpose a rchisq density line over a histogram with something like: hist(alnlength) lines(density(rchisq(length(alnlength), 4)),col="red") But the rchisq line won't appear anywhere, Anyone knows what I am missing here? Thanks in advance, Albert.

Hi, I have a least square fitting problem with linear inequality constraints. pcls seems capable of solving it so I tried it, unfortunately, it is stuck with the following error: > M <- list() > M$y = Dmat[,1] > M$X = Cmat > M$Ain = as.matrix(Amat) > M$bin = rep(0, dim(Amat)[1]) > M$p=qr.solve(as.matrix(Cmat), Dmat[,1]) > M$w = rep(1, length(M$y)) > M$C = matrix(0,0,0)

Hi, I use a Debian Sarge with samba 3.0.23c (from testing) and cups 1.1.23. For unexpected reason, some printers are sometime stopped (pause), and I need to start then again. There is no trace of a "stop" in the log of cups. (and I'm the only one who can pause them, anyway). here is an example : cat 10.lpr.log | grep imp-cao | grep -E "(stop|start)" Nov 10 14:15:41