search for: 10,10

Dear all; Does anyone know how to add grid lines to a persp plot? I've tried using lines(trans3d..) but the lines of course are superimposed into the actual 3d surface and what I need is something like the plot shown in the following link: http://thermal.gg.utah.edu/tutorials/matlab/matlab_tutorial.html I'll appreciate any ideas Thanks PM

Dear all, I have a question about the plotting of the dynamical orbit by using rgl. My code is as follows: open3d() par3d(cex=2) bg3d("white") for(i in 1:out.length){ ind<-which plot3d(x[1:i],y[1:i],out[,"z"][1:i],xlim=c(-10,10),ylim=c(-10,10),zlim=c(-10,10),col="red",cex=2) Sys.sleep(0.0001) } ,where x,y,out[,"z"] are the position of the object in different time. I would like to show the orbit evoluting with time. However, I can just get animation with flash windows. I don't know how to use pla...

Full_Name: Craig A. McKinstry Version: 1.0.0 OS: WinNT4.0 and Win98 Submission from: (NULL) (192.101.100.130) When creating a multi-panel graphic, the command par(mfg=c(i,j,r,c)) is supposed to allow the user free movement between graphics panels to develop each panel separately. This works for the first invocation of par(mfg=c()) and for the first panel panel specified, but not for subsequen...

...ough (1,4) and (-2,-5) of the circle, but on Windows XP, R 2.20 the drawing was not good at all, and the known solutions were not shown in the graph. I guess I got it wrong? Is this use not intended ? Best wishes Troels Ring, MD Aalborg, Denmark symbols(x=1,y=-1,circles=5,inches=FALSE, xlim=c(-10,10),ylim=c(-10,10)) curve(3*x+1,-10,10,1000,add=T) abline(v=c(-2,1)) abline(h=c(-1,-5,4))

Hello all, Maybe I'm being thick, but I was trying to figure out a simple way to create a list with the same dimension of another list, but populated with NA values. masterlist = list( aa=list( a=matrix(rnorm(100),10,10), b=matrix(rnorm(130),13,10), c=matrix(rnorm(140),14,10)), bb=list( a=matrix(rnorm(100),10,10), b=matrix(rnorm(110),11,10), c=matrix(rnorm(120),12,10)), cc=list(...

...ment(s) (uv.coords(u, v)) > plot.new <- function() + {} > { } NULL > plot.new <- function() + { } > { } NULL > { for } Error: unexpected '}' in "{ for }" > area = 60*20 > lambda = 2 > N = rpois(1,lambda*area) > u = runif(N,20,40) > v = runif(N,-10,10) > plot(u,v,asp=1) > > plot(u,v,asp=1) > u = runif(N,20,40) > v = runif(N,-10,10) > plot(u,v,asp=1) > > plot(u,v,asp=1) > u = runif(N,20,40) > v = runif(N,-10,10) > plot(u,v,asp=1) > u = runif(N,20,40) > v = runif(N,-10,10) > plot(u,v,asp=1) > h = c...

...ug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@biostat.ku.dk # ###################################################### In version 0.99.0 for Windows NT I get the following: > .Random.seed [1] 0 21012 4461 28299 > rbinom(10,10,.3) [1] 4 5 5 1 5 2 4 2 6 2 > rbinom(10,9,.3) [1] 0 0 0 0 0 0 0 0 0 0 > rbinom(10,10,.3) [1] 3 3 2 2 3 1 3 3 3 2 > rbinom(10,11,.3) [1] 0 0 0 0 0 0 0 0 0 0 I tried similar inputs in version 0.65.1: > .Random.seed [1] 0 16049 70 30225 > rb...

...Cluster(8, type = "SOCK")) but that seems to defeat the point. library(foreach) library(doSNOW) registerDoSNOW(makeCluster(8, type = "SOCK")) getDoParWorkers() getDoParName() getDoParVersion() mySamples <- foreach (jj = 1:4, .combine=cbind) %dopar% { return(sample(1:10,10,replace=TRUE)) } mySamples ########## r 2.14.1 ########## > library(foreach) > library(doSNOW) > registerDoSNOW(makeCluster(8, type = "SOCK")) > getDoParWorkers() [1] 8 > getDoParName() [1] "doSNOW" > getDoParVersion() [1] "1.0.6" > >...

...h each line labeled with a different letter. But I'm not able to make the legend display the same kind of pattern '-a-', instead the letter is overwritten by the line. A simpler legend with only the letter is not very visible and the pt.bg does nothing with letters. Any idea? plot(1:10,10:1,lty=1,type='b', lwd=2,pch='a') legend("left", legend=c("ds1","ds2"), bty='n', col=1:2, lty=2,lwd=4,pch=letters) Thanks for your help! mario -- Ing. Mario Valle Data Analysis and Visualization Group | http://www.cscs.ch/~m...

Hello friends!! I have two questions, and I would like that you could answer me!!! I have created a plot as plot(range(10,-10),range(10,-10),col="blue",col.axis="blue",col.lab="blue",col.main="blue",col.sub="blue"); 1?) I want that the square of plot and the lines that indicates the value -10,-5,0,5,10 have a blue colour too. 2?)It?s different to the question one. I...

I have the following questions about the variance of the random effects in the survreg() function in the survival package: 1) How can I extract the variance of the random effects after fitting a model? For example: set.seed(1007) x <- runif(100) m <- rnorm(10, mean = 1, sd =2) mu <- rep(m, rep(10,10)) test1 <- data.frame(Time = qsurvreg(x, mean = mu, scale= 0.5, distribution = "weibull"), Status = rep(1, 100), Unit = gl(10,10) ) mod1 <...

...ument(s) (uv.coords(u, v)) > plot.new <- function() + {} > { } NULL > plot.new <- function() + { } > { } NULL > { for } Error: unexpected '}' in "{ for }" > area = 60*20 > lambda = 2 > N = rpois(1,lambda*area) > u = runif(N,20,40) > v = runif(N,-10,10) > plot(u,v,asp=1) > > plot(u,v,asp=1) > u = runif(N,20,40) > v = runif(N,-10,10) > plot(u,v,asp=1) > > plot(u,v,asp=1) > u = runif(N,20,40) > v = runif(N,-10,10) > plot(u,v,asp=1) > u = runif(N,20,40) > v = runif(N,-10,10) > plot(u,v,asp=1) > h = chu...

Hi, I am trying to use the optimize function to optimize a function. The results I am getting don't agree with what I compute on my own and when I look at the graph of f(x) = 100 + ((x-10)**2 + (x-10)) * cos(x-10), where -10 <= x <= 10 in gnuplot. I suspect I am making a mistake in the usage of the R optimize function, perhaps someone could point out where? > f<-function(x) 100 + ((x-10)**2 + (x-10)) * cos(x-10) to MAXIMIZE: > result=optimize(f, c(-1...

Hello R-help, I'm trying to make a fairly simple plot axis that goes something like this: plot(-10:10,-10:10, yaxt='n') axis(side=2, las=1, hadj=1, tck=-.01, cex.axis=.6) ...but as you can see, the labels are not close enough to the y-axis (where I want them... to save space for publication). Can anybody help me figure out how to move these labels over the the right a bit? Thanks, -...

Hi everyone This must be a FAQ but I can't find it anywhere... I want a postscript image of a contour() plot, with axes of equal length. Try R> postscript(file="~/f.ps") R> contour(matrix(rnorm(100),10,10)) R> dev.off() This isn't what I want: the plotting region is, as documented, quarter of an inch shy of the paper edge and the axes appear to be different lengths. contour() doesn't take a "asp" argument. postscript() does have a width and a length argument but I can...

Hi! I would like to post the following question: I was trying to figure out how to do the simulation shown in Fig 10.6 of John Verzani's book 'Using R for Intro Statistics'. It is on page 290, with a description on the previous page. It seems like a simple thing... Just needing to duplicate a procedure. (Perhaps I need to do it with a loop?) This is what I was trying: replicate(100, { abline(lm(y~x)...

Dear R users, suppose I have a matrix of observations for which I calculate all pair-wise correlations: m=matrix(sample(1:100,replace=T),10,10) w=cor(m,use="pairwise.complete.obs") How do I extract only those correlations that are >0.6? w[w>0.6] #obviously doesn?t work, and I can?t find a way around it. I would very much appreciate any help! Best wishes Christoph (using R 2.5.1 on Windows XP) -...

Hi, I plotted 12 graphs on a page and output to a png file. I wanted to have an overall title for all 12 graphs. What command can I use to do this? Below is the code that plotted the 12 graphs in one page. # FM10 by Month/ Export the plot to Wash2005FM10.png png(file="Wash2005FM10.png",bg="white") par(mfrow = c(3,4)) # Plot 12 Month of OFM10, FFM10 for(i in 1:12) { Temp <- subset (Wash2005, MM == i) plot (Temp$FM10.1, Temp$FM10,main=month.name[i],xlim=c(5,25),ylim=c(5,25)) }...

Maybe, the following is not a bug in the strict sense,.... The following code shows what the effect : p1 <- function(...) plot.default(...) p2 <- function(...) { n <- names(list(...)); plot.default(...) } par(mfrow=c(2,2)) p1(1:10, panel.first = grid(10,10)) # works okay (grid drawn) p2(1:10, panel.first = grid(10,10)) # doesn't draw the grid A consequence of this is that the panel.first argument (of plot.default) doesn't work with matplot() [well, matplot()'s doc doesn't claim it should, but still -----...

...add another context i.e 'internal' and the asterisk is complaining for not finding the required extension in the 'default' context. Asterisk schould point this to the internal contex and not default. here my simply config data : sip.conf : [general] port=5060 bindaddr=0.0.0.0 [10] type=friend secret=1234 host=dynamic context=internal [11] type=friend secret=1234 host=dynamic context=internal extensions.conf [default] exten =>2,1,Playback(digits/2) ; exten =>2,2,Goto(default,10,1) exten=>3,1,Playback(pbx-invalid) exten=3,2,Goto(default,4,1) exten=4,1,Playback(vm...