similar to: joining two points in rgl

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 [[alternative HTML version deleted]]

If my sample size is small is there a particular switch option that I need to use with t.test so that it calculates the t ratio correctly? Here is a dummy example? á =0.05 Mean pain reduction for A =27; B =31 and SD are SDA=9 SDB=12 drgA.p<-rnorm(5,27,9); drgB.p<-rnorm(5,31,12) t.test(drgA.p,drgB.p) # what do I need to give as additional parameter here? I can do it manually but

Dear useRs, I'm trying to simply fill in the area under a curve using RGL. Here' the set up: x <- c(0.75,75.75,150.75,225.75,300.75,375.75,450.75,525.75,600.75,675.75, 0.5,50.5,100.5,150.5,200.5,250.5,300.5,350.5,400.5,450.5, 0.25,25.25,50.25,75.25,100.25,125.25,150.25,175.25,200.25,225.25) y <- c(0.05,4.91,9.78,14.64,19.51,24.38,29.24,34.11,38.97,43.84,

Hi all, How can I tell RGL to set the center for the rotation to the origin of the coordinate system (0,0,0). It seems that the default is to use the center of the display not the origin of the coordinate system. open3d() lines3d(c(0, 1), c(0,0), c(0,0)) lines3d(c(0,0), c(0, 1), c(0,0)) lines3d(c(0,0), c(0,0), c(0, 1)) TIA Mark ??????????????????????????????????????? Mark Heckmann Blog:

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[1]> cos(direction[i]) [1] 6.123032e-17 Browse[1]> cos(-1.570796) [1] 3.267949e-07

Can I insert an element in an array at a particular position without destroying the already existing element? X<-c(1,2,3,4,5) I want to insert an element between 1 and 2. Thanks ../Murli [[alternative HTML version deleted]]

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

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 the help archive but did

Is there a function that calculates the inter-rater agreement index (kappa) in R? Thanks ../Murli [[alternative HTML version deleted]]

Hi, I am having an issue when creating a postscript file from RGL window. It seems to cut off some of the axis labels. Here is the code I am using. I created a 3D plot using RGL_0.71 with R 2.5 on Windows XP. z1<-c(5,4,1,4.5,2,3,2,1,1) z2<-c(6,8,7,7.5,5,3.5,4,1,1) z3<-c(3,2,4,7,3,4.5,6,2,3) x1<-seq(1,9) x2<-seq(1,9) x3<-seq(10,18) y1<-seq(8,0) y2<--1*y1

Which is the parameter that is used to decrease the size of ylabs plotted in biplot? I tried playing with cex and cex.lab I am not getting it right pc <- princomp(USArrests) biplot(pc, xlabs = rep("", nrow(USArrests)),ylabs=(colnames(USArrests))) Thanks../Murli [[alternative HTML version deleted]]

I used the multcomp package sometime back for doing multiple comparisons. I see that it has been updated and the methods like simint are no longer supported. When I run the program it prompts to me to use glht. How do I get the lower and upper conf int and the pValues using glht? Does anyone have an example? Thanks ../Murli [[alternative HTML version deleted]]

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 fine in this case gene.rbf <-

I am interested in plotting histograms for the following data Isoform Tumor_65_198 Tumor_50_192 Tumor_80_167 Tumor_80_204 Tumor_95_197 Tumor_70_189 Tumor_90_202 Tumor_40_177 Tumor_60_21 Tumor_70_174 Tumor_70_147 Tumor_50_5 ABCC4-2007 1 1 1 6 1 9 10 1 2 0 10 1 ABCC4-2008 5 8 7 5 3 10 5 5 7 3 10 3 ABCC4-2009 0 0 0 0 0 0 0 0 0 0 0 0 ABCC4-2010

For a paper dealing with generalized ellipsoids, I want to illustrate in 3D an ellipsoid that is unbounded in one dimension, having the shape of an infinite cylinder along, say, z, but whose cross-section in (x,y) is an ellipse, say, given by the 2x2 matrix cov(x,y). I've looked at rgl:::cylinder3d, but don't see any way to make it accomplish this. Does anyone have any ideas? thx,

I am trying to add text at specific location on my graph. I know this can be done in R but I can't recollect. I was trying using locator() to identify the position and use identify() but I can get it to work. Can someone jog my memory? Thanks ../Murli [[alternative HTML version deleted]]

I am running a recursive feature selection wrapper on svm for a large data set. The routine gets progressively slower. I am deleting the objects not required. Is there anything else you advise me look for? I tried running it from the command line but it does not print anything on the screen. I did not specify any output file and was expecting the screen to be the default. Can any one help me

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, -0.25,1) Uninformative.data.class2<-rnorm(900,0,1) The above will give me one set of data

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

Here it is again, hope this is more clear I am using the following data (only a small subset is given): Habitat Fungus.yield Birch 20.83829053 Birch 22.9718181 Birch 22.28216829 Birch 24.23136797 Birch 22.32147961 Birch 20.30783598 Oak 27.24047258 Oak 29.7730014 Oak 30.12608508 Oak 25.76088669 Oak 30.14750974 Hornbeam 17.05307949 Hornbeam 15.32805111 Hornbeam 18.26920177 Hornbeam 21.30987049