Re: [R] error while plotting

2009-09-29 Thread Nair, Murlidharan T
-Original Message- From: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Sent: Tuesday, September 29, 2009 4:52 AM To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] error while plotting Nair, Murlidharan T wrote: > > -Original Message- > From: U

Re: [R] 3D to 2D projection

2009-09-28 Thread Nair, Murlidharan T
[mailto:dwinsem...@comcast.net] Sent: Monday, September 28, 2009 2:57 PM To: Nair, Murlidharan T Cc: r-h...@stat.math.ethz.ch Subject: Re: [R] 3D to 2D projection On Sep 28, 2009, at 1:41 PM, Nair, Murlidharan T wrote: > David, > Have you used persp or trans3d before? Yes, both. persp require

Re: [R] 3D to 2D projection

2009-09-28 Thread Nair, Murlidharan T
Thank you for this info. I think this has what I need. Cheers../murli -Original Message- From: hadley wickham [mailto:h.wick...@gmail.com] Sent: Monday, September 28, 2009 2:27 PM To: Nair, Murlidharan T Cc: David Winsemius; r-h...@stat.math.ethz.ch Subject: Re: [R] 3D to 2D projection

Re: [R] error while plotting

2009-09-28 Thread Nair, Murlidharan T
-Original Message- From: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Sent: Sunday, September 27, 2009 1:17 PM To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] error while plotting Nair, Murlidharan T wrote: > I am getting the following errors when I am try

Re: [R] 3D to 2D projection

2009-09-28 Thread Nair, Murlidharan T
r$z[(endPt1+1):(endPt1*2)],type="l") s3d$points3d(dnaStr$x[(endPt1*2+1):(endPtlength)],dnaStr$y[(endPt1*2+1):(endPtlength)],dnaStr$z[(endPt1*2+1):(endPtlength)],type="l") -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Saturday, September 26, 2009

Re: [R] error while plotting

2009-09-28 Thread Nair, Murlidharan T
, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] error while plotting Nair, Murlidharan T wrote: > I am getting the following errors when I am trying to plot the data below. I > cannot figure out the error. > Error in plot.window(...) : need finite 'xlim' values > In addi

[R] 3D to 2D projection

2009-09-26 Thread Nair, Murlidharan T
Is there a method that I can use to convert 3D coordinates into 2D? I was looking at persp and trans3d. Are those the ones I should be looking at ? Thanks ../Murli __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] error while plotting

2009-09-25 Thread Nair, Murlidharan T
I am getting the following errors when I am trying to plot the data below. I cannot figure out the error. Error in plot.window(...) : need finite 'xlim' values In addition: Warning messages: 1: In min(x) : no non-missing arguments to min; returning Inf 2: In max(x) : no non-missing arguments to ma

Re: [R] multicomp plotting

2009-09-25 Thread Nair, Murlidharan T
suggestion. Thanks .../Murli From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of Nair, Murlidharan T [mn...@iusb.edu] Sent: Thursday, September 24, 2009 2:06 PM To: r-help@r-project.org Subject: [R] multicomp plotting I

[R] multicomp plotting

2009-09-24 Thread Nair, Murlidharan T
I am trying to plot my multiple comparison data. Can anyone give me some input of the error I am getting. The data and code is appended below. Thanks ../Murli library(multcomp) sig.data<-structure(list(X = 1:63, Cell.lines = structure(c(1L, 6L, 13L, 25L, 33L, 42L, 2L, 7L, 14L, 26L, 34L, 43L

Re: [R] dotchart to barplots

2009-09-24 Thread Nair, Murlidharan T
___ From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of Nair, Murlidharan T [mn...@iusb.edu] Sent: Wednesday, September 23, 2009 5:11 PM To: Greg Snow; Henrique Dallazuanna Cc: r-help@r-project.org Subject: Re: [R] dotchart to barplots

Re: [R] dotchart to barplots

2009-09-23 Thread Nair, Murlidharan T
Thanks Greg. I was also thinking about it after I saw my plots. Cheers../Murli -Original Message- From: Greg Snow [mailto:greg.s...@imail.org] Sent: Wednesday, September 23, 2009 3:49 PM To: Nair, Murlidharan T; Henrique Dallazuanna Cc: r-help@r-project.org Subject: RE: [R] dotchart to

Re: [R] dotchart to barplots

2009-09-23 Thread Nair, Murlidharan T
Guess, I miss the argument when I ?barplot. Cheers../Murli -Original Message- From: Peter Ehlers [mailto:ehl...@ucalgary.ca] Sent: Wednesday, September 23, 2009 3:56 PM To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] dotchart to barplots Murli, Two points: 1. I

Re: [R] dotchart to barplots

2009-09-23 Thread Nair, Murlidharan T
Thanks Peter. Where did you find that option? It's really cool Cheers../Murli -Original Message- From: Peter Ehlers [mailto:ehl...@ucalgary.ca] Sent: Wednesday, September 23, 2009 3:56 PM To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] dotchart to barplots

Re: [R] dotchart to barplots

2009-09-23 Thread Nair, Murlidharan T
The journal wants black and white only :) -Original Message- From: Henrique Dallazuanna [mailto:www...@gmail.com] Sent: Wednesday, September 23, 2009 3:41 PM To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] dotchart to barplots Perhaps a white border: barplot(t

Re: [R] dotchart to barplots

2009-09-23 Thread Nair, Murlidharan T
ct.data[,1], cex.axis = 0.7, cex.names = 0.7, las=2) On Wed, Sep 23, 2009 at 4:26 PM, Nair, Murlidharan T wrote: > Ok, I could make it perpendicular by specifying las=2 > > barplot(t(as.matrix(intersect.data[,2:5])), >          beside = T, horiz = T, >          names.arg = intersect

Re: [R] dotchart to barplots

2009-09-23 Thread Nair, Murlidharan T
- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Nair, Murlidharan T Sent: Wednesday, September 23, 2009 3:21 PM To: Henrique Dallazuanna Cc: r-help@r-project.org Subject: Re: [R] dotchart to barplots I had tried names.arg=c(intersect.data[,1]) so that was the

Re: [R] dotchart to barplots

2009-09-23 Thread Nair, Murlidharan T
: Henrique Dallazuanna [mailto:www...@gmail.com] Sent: Wednesday, September 23, 2009 3:09 PM To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] dotchart to barplots Try this: barplot(t(as.matrix(intersect.data[,2:5])), beside = T, horiz = T, names.arg

[R] dotchart to barplots

2009-09-23 Thread Nair, Murlidharan T
Hi, I am trying to plot the following data so that it can be visually represented well. I tried the dotchart but I felt it was too spread out. Then I tried the barplot which is good enough for me. Is there a way to give the labels for the y-axis as in the dot chart? Also, I feel the grey level

Re: [R] joining two points in rgl

2009-08-20 Thread Nair, Murlidharan T
pts[2,],text=c("t"), adj = c(0.5, -1)) shade3d(cylinder3d(pts, radius=0.1, e2=rbind(c(1,0,0),c(1,0,0))), col="red") Thanks for adj. The rgl package is absolutely wonderful. It is extremely useful to me. Cheers../Murli ____ From: Duncan Murdo

Re: [R] joining two points in rgl

2009-08-19 Thread Nair, Murlidharan T
so that if I write the text at the exact same coordinate of the sphere the text is above the sphere and can be seen? Thanks ../Murli -----Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Nair, Murlidharan T Sent: Friday, August 14, 200

Re: [R] joining two points in rgl

2009-08-14 Thread Nair, Murlidharan T
: Thursday, August 13, 2009 8:24 PM To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] joining two points in rgl Nair, Murlidharan T wrote: > Hi!! > > I need to draw a cylinder/tube joining two points. I am trying to make > something presentable, I have been able to do it u

[R] joining two points in rgl

2009-08-13 Thread Nair, Murlidharan T
Hi!! I need to draw a cylinder/tube joining two points. I am trying to make something presentable, I have been able to do it using lines3d. But is it possible to increase the thickness of the lines? The size parameter does not work. Does any one have any suggestions? Thanks ../Murli library

Re: [R] some help with plotting a beautiful structure using rgl

2009-08-13 Thread Nair, Murlidharan T
Ok, I was able to figure out one part by using rgl.spheres. I am still trying to insert text using rgl.texts. Does any one have a simple example? -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Nair, Murlidharan T Sent: Thursday

[R] some help with plotting a beautiful structure using rgl

2009-08-13 Thread Nair, Murlidharan T
Hi!! I need some help in using the correct required symbol when plotting my DNA structure using rgl. I need to plot one strand using empty circles and the other using filled circles and if possible have sequence in the center. I tried using pch and lty but have not been able to get it to do wha

Re: [R] computing the radius of an arc

2009-07-27 Thread Nair, Murlidharan T
Greg Snow [mailto:greg.s...@imail.org] Sent: Friday, July 24, 2009 12:41 PM To: Nair, Murlidharan T; Hans W Borchers; r-h...@stat.math.ethz.ch; Bert Gunter; 'Gabor Grothendieck' Subject: RE: [R] computing the radius of an arc There is a function "rsaga.local.morphometry" in the RSAGA

Re: [R] probability on a barplot

2009-07-27 Thread Nair, Murlidharan T
Kindly, post the solution to the problem, so that it will benefit others. An example could would be great. Cheers../Murli -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Erin Hodgess Sent: Monday, July 27, 2009 1:33 PM To: R help

Re: [R] computing the radius of an arc

2009-07-24 Thread Nair, Murlidharan T
n(0.46365) proj2d.xy<-cbind(proj2D.x,proj2D.y) colnames(proj2d.xy)<-c("x","y") -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Hans W Borchers Sent: Friday, July 24, 2009 5:41 AM To: r-h...@stat.math.ethz.ch Subje

[R] computing the radius of an arc

2009-07-23 Thread Nair, Murlidharan T
Hi!! I am interesting in computing the radius of an arc that best approximates a curve. Is there an R function that I can use to draw an arc? Nothing useful came up when I searched help.search. Does anyone have any suggestion to do this? Thanks ../Murli ___

Re: [R] Matrix multiplication precision

2009-07-15 Thread Nair, Murlidharan T
n expert in the area. Cheers../Murli -Original Message- From: dmba...@gmail.com [mailto:dmba...@gmail.com] On Behalf Of Douglas Bates Sent: Wednesday, July 15, 2009 7:29 AM To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] Matrix multiplication precision On Wed,

[R] Matrix multiplication precision

2009-07-14 Thread Nair, Murlidharan T
Hi!! I am trying to multiply 5 matrices and then using the inverse of that matrix for further computation. I read about precision problems from the archives and the suggestion was to use as.numeric while computing the products. I am still having problems with the results. Here is how I am using

Re: [R] Trig functions strange results

2009-07-14 Thread Nair, Murlidharan T
To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] Trig functions strange results FAQ 7.31 This is what happens with floating point number and you are only printing out 7 digits of precision; look at the results > cos(-1.570796) [1] 0.003267949 > cos(-1.5707961)

[R] Trig functions strange results

2009-07-14 Thread Nair, Murlidharan T
I am trying to calculate coordinate transformations and in the process of debugging my code using debug I found the following Browse[1]> direction[i] [1] -1.570796 Browse[1]> cos(direction[i]) [1] 6.123032e-17 Browse[1]> cos(-1.570796) [1] 3.267949e-07 Browse[1]> direction[i] [1] -1.570796 Browse

[R] strsplit

2009-05-29 Thread Nair, Murlidharan T
I am trying split a string and use one part of it to label graphs. I am using strsplit for that. While I am able to split it, how do I access the separated parts. > filName<-"chrI_2223" >part<- strsplit(filName,"\\_") > part [[1]] [1] "chrI" "2223" > part[1] [[1]] [1] "chrI" "2223" I looked up

Re: [R] scatterplot3d

2009-05-14 Thread Nair, Murlidharan T
-Original Message- From: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Sent: Thursday, May 14, 2009 4:44 AM To: Nair, Murlidharan T Cc: r-help@r-project.org Subject: Re: [R] scatterplot3d Nair, Murlidharan T wrote: > Does anyone have any suggestion on this? I tried par(new

Re: [R] scatterplot3d

2009-05-13 Thread Nair, Murlidharan T
Does anyone have any suggestion on this? I tried par(new="T") It does not seem to work with scatterplot3d. Any suggestion is appreciated. Cheers../Murli -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Nair, Murlidha

[R] scatterplot3d

2009-05-13 Thread Nair, Murlidharan T
Hi!! I am interested in plotting 3 different (x,y,z) coordinate sets on the same plot. How do I achieve using scatterplot3d? Is there an easy method to call to change the projection or rotate the figure along a particular axis? The following is the code that plots if one at a time. Any help i

[R] boxplot labelling levels

2008-09-16 Thread Nair, Murlidharan T
I want the levels to appear in the boxplot instead of 1 and 2. What do I need to do for that? Here is the dummy code. x<-runif(100,50,80) x1<-runif(100,70,80) True.positives<-c(x,x1) splice<-factor(c(rep("Human.AA.200",100),rep("Human.AA.100",100))) splice<-factor(splice,levels=c("Human.AA.200

[R] ksvm accessing the slots of S4 object

2008-09-14 Thread Nair, Murlidharan T
I am using kernlab to build svm models. I am not sure how to access the different slots of the object. For instance if I want to get the nuber of support vectors for each of model I am building and store it in a vector. >ksvm.model <- ksvm(Class ~ ., data = somedata,kernel = "vanilladot", cross

Re: [R] how to draw a vertical line from points to x-axis

2008-09-08 Thread Nair, Murlidharan T
Anny, You can also do the following plot(0:10, 0:10, pch=16, type="h") Cheers../Murli -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Peter Alspach Sent: Sunday, September 07, 2008 5:10 PM To: Anny Huang; r-help@r-project.org Subject: Re: [R] how to dra

[R] data frame question

2008-09-03 Thread Nair, Murlidharan T
I have a data frame containing sequences and I am interested in changing a few sequences in a window and the swapping the original sequence back after I have completed my analysis. My temporary data frame that I am creating seq.in.window does not like the way I am making me assignment. The vari

Re: [R] Text

2008-09-02 Thread Nair, Murlidharan T
Dani, You can put text on your graphs using the locator function. Here is an example hist(rnorm(1000,0,1) #plots a histogram for you xy.cord<-locator(1) #Once you give this click on the graph to get the coordinates text(xy.cord, "Text you want to write") Hope this helps. Cheers../Murli -

[R] Error in .local(object, ...) : test vector does not match model !

2008-09-02 Thread Nair, Murlidharan T
I am getting a really strange error when I am using predict on an ksvm model. The error is "Error in .local(object, ...) : test vector does not match model !". I do understand that this happens when the test vectors do not match the Model. But in this case it is not so. I am attaching a portion

Re: [R] heatmap

2008-07-02 Thread Nair, Murlidharan T
___ From: [EMAIL PROTECTED] [EMAIL PROTECTED] On Behalf Of Nair, Murlidharan T [EMAIL PROTECTED] Sent: Wednesday, July 02, 2008 2:02 PM To: r-help@r-project.org Subject: [R] heatmap I am interested in plotting a matrix of data using heatmap to reflect frequency of positional occu

[R] heatmap

2008-07-02 Thread Nair, Murlidharan T
I am interested in plotting a matrix of data using heatmap to reflect frequency of positional occurrence of certain strings. I am not interested in clustering them. I also want to make a key for the colors.Can someone help me with that. I am appending a portion of my data here. Thanks ../Murli

Re: [R] stringdot

2008-06-25 Thread Nair, Murlidharan T
e vector of size 1.2 Gb. Is my usage wrong ?. I am running this under linux serve 6850 which has about 8Gb memory Any would be appreciated. Cheers../Murli From: [EMAIL PROTECTED] [EMAIL PROTECTED] On Behalf Of Nair, Murlidharan T [EMAIL PROTECTED] Sent: We

[R] stringdot

2008-06-25 Thread Nair, Murlidharan T
Hi!! I am trying to figure out how to use the string kernel "stringdot" in kernlab. k <- function(x,y) { (sum(x*y) +1)*exp(-0.001*sum((x-y)^2)) } class(k) <- "kernel" data(promotergene) ## train svm using custom kernel gene.k <- ksvm(Class~.,data=promotergene,kernel=k,C=10,cross=5) # works

Re: [R] dotchart

2008-06-04 Thread Nair, Murlidharan T
Thanks! I put the labels as follows dotchart(as.matrix(intersect.data[-1]), labels=intersect.data[,1], cex=0.8) Murli From: Henrique Dallazuanna [EMAIL PROTECTED] Sent: Wednesday, June 04, 2008 1:49 PM To: Nair, Murlidharan T Cc: r-help@r-project.org

[R] dotchart

2008-06-04 Thread Nair, Murlidharan T
I am trying to plot the following data using dotchart intersect.data<-structure(list(X = structure(c(1L, 3L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 2L, 4L, 5L, 6L, 7L), .Label = c("1-100", "1001-1100", "101-200", "1101-1200", "1201-1300", "1301-1400", "1401-1500", "201-300", "301-400", "401-500", "

[R] for loop step

2008-05-20 Thread Nair, Murlidharan T
How do I define the incremental step in a "for" loop? for (j in 1:10){ cat(j, "\n") } In the above example, if I want to increment j by 2 where do I specify that? Thanks [[alternative HTML version deleted]] __ R-help@r-project.org mail

[R] listing components of an object

2008-03-02 Thread Nair, Murlidharan T
Is there a method to list the components of an object, instead of looking at the help for that method? Let me be more clear with an example data(iris) ## tune `svm' for classification with RBF-kernel (default in svm), ## using one split for training/validation set obj <- tune(svm, Species

Re: [R] filling data into objects

2008-02-06 Thread Nair, Murlidharan T
Oh yes, that's what I was looking for. I can then do a Class1<-rbind(Informative.data, Uninformative.data) Thanks ../Murli -Original Message- From: Henrique Dallazuanna [mailto:[EMAIL PROTECTED] Sent: Wednesday, February 06, 2008 3:06 PM To: Nair, Murlidharan T Cc: [EMAIL P

[R] filling data into objects

2008-02-06 Thread Nair, Murlidharan T
I am trying to generate artificial data for feature selection. Basically trying to generate a total of 1000 features with 100 that are informative and rest are uninformative. Informative.data.class1<-rnorm(100,0.25,1) Uninformative.data.class1<-rnorm(900,0,1) Informative.data.class2<-rnorm(100,

[R] Non-independent paired measures

2007-10-29 Thread Nair, Murlidharan T
I am interested in generating simulated data for demonstrating sampling distribution of the mean difference between non-independent paired measures. Is there an elegant way do generate non-independent data sets? Thanks ../Murli [[alternative HTML version deleted]] ___

Re: [R] multicomp plotting

2007-10-24 Thread Nair, Murlidharan T
ay no adjustments to be applied or does it have to do with attributes typeof? I would appreciate if you could clarify this for me. Thanks. Cheers../Murli -Original Message- From: Richard M. Heiberger [mailto:[EMAIL PROTECTED] Sent: Wednesday, October 24, 2007 12:34 AM To: Nair,

[R] multicomp plotting

2007-10-23 Thread Nair, Murlidharan T
I am trying to plot a large set of multiple comparison data from a csv file. I used to do in the following manner, but its giving me an error now, and I am not sure what correct. I have put the data in the code using the output from dput for completeness. Can anyone help me with my mistakes? Tha

[R] correlated data

2007-10-17 Thread Nair, Murlidharan T
Hi!! I am trying to generate data with specific correlation using corgen method from the ecodist package. It does not seem to generate the data properly. I am listing the code below. Secondly, when I put the code given below all in one file and source it, it is unable to complete the task. But

Re: [R] combining vectors on unequal length

2007-10-05 Thread Nair, Murlidharan T
Thanks to all for their suggestions!! Cheers../Murli -Original Message- From: Gabor Grothendieck [mailto:[EMAIL PROTECTED] Sent: Thursday, October 04, 2007 3:30 PM To: Nair, Murlidharan T Cc: [EMAIL PROTECTED] Subject: Re: [R] combining vectors on unequal length Assuming they are

[R] combining vectors on unequal length

2007-10-04 Thread Nair, Murlidharan T
If I have two vectors X<-1:10 Y<-1:5 When I combine them using cbind, the shorter one is repeated and both are made of the same length. Is there a methods that does this without duplicating the shorter one. I want to use this to store the data back to a file. Thanks ../Murli [[alterna

Re: [R] MAD

2007-09-18 Thread Nair, Murlidharan T
-Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Monday, September 17, 2007 5:10 PM To: Nair, Murlidharan T Cc: [EMAIL PROTECTED] Subject: Re: [R] MAD On 9/17/07, Nair, Murlidharan T <[EMAIL PROTECTED]> wrote: > > I am calculating the median absolu

[R] MAD

2007-09-17 Thread Nair, Murlidharan T
I am calculating the median absolute deviation using mad function, and it tends to ignore the parameter constant=1, when I am calculating it for x=seq(1:5). Am I missing something here? x<-seq(1:5) mad(x)# gives [1] 1.4826 mad(x, constant=1)# gives [1] 1 #Here is the long form dev.from.median<-ab