similar to: How to plot a blank plot

Hello everyone, I had a few questions that I have not been able to figure out despite a lot of reading. 1) Adding a title to a multiplot figure: If I plot multiple plots with >par(mfrow=c(2,6)) how do I add an overall title to the figure, not the individual plots? 2) Saving histories on the fly: Is there a .Rhistory equivalent to >save.image() ? I would like to be able to save the

Estimados Quisiera tener la sintaxis generar una serie de números entre 0 y 1 con distribución gamma de 25 puntos y que la suma de todos los puntos debe dar 1. Por ejemplo: c (0.000, 0.001, 0.012, 0.043, 0.078, 0.104, 0.117, 0.116, 0.108, 0.094, 0.078, 0.063, 0.049, 0.038, 0.028, 0.021, 0.015, 0.011, 0.008, 0.005, 0.004, 0.003, 0.002, 0.001, 0.001) No hemos podido

Estimado Jorge Perdóneme que lo moleste de nuevo, hay otra condición además de que sum(y)=1 y es que y[1] tiene que dar 0 en el ejemplo y<- c (0.0000000000, 0.6321985783, 0.2325728597, 0.0855587737, 0.0314753138, 0.0115791209, 0.0042597205, 0.0015670636, 0.0005764905, 0.0002120790) y[1]=0 sum(y)=1 esto se utiliza en el paquete EpiEstim donde se

Hi, Can anyone help to write a more elegant version of my code? I am sure this can be put into a loop but I am having trouble creating the objects b1,b2,b3,...,etc. b1 <- rigamma(50,1,1) theta1 <- rgamma(50,0.5,(1/b1)) sim1 <- rpois(50,theta1) b2 <- rigamma(50,1,1) theta2 <- rgamma(50,0.5,(1/b2)) sim2 <- rpois(50,theta2) b3 <- rigamma(50,1,1) theta3 <-

Hi, I'm trying to add points to 2 plots on the fly using par(mfg=vector) so switch between them. However, the appropriate scales aren't switched when changing from one plot to another, e.g. par(mfcol=c(2,1)) plot(1,1, col="blue") # blue plot plot(1.2,1.2, col="red") # red plot points(1.1,1.1) # appears to bottom left of red point par(mfg=c(1,1)) #

Hi, List, I used the following codes to generate ps plots but foo.ps contains nothing. Would someone please point out what is wrong with my codes? Thanks a million! postscript('foo.ps') par(mfrow=c(2,1)) par(mfg=c(1,1)) hist(rnorm(100),col='blue') par(mfrow=c(2,2)) par(mfg=c(2,1)) hist(rnorm(50),col='blue') par(mfg=c(2,2)) hist(rnorm(60),col='blue') dev.off()

Dear List, I have the following function where I plot on one page a histogram and a boxplot. I use option labels = TRUE in hist to show labels. However the label on the highest bar in the histogram is not showing entirely. How can I correct this? Thanks for any suggestions? histtst <- function(n = 100,...){ set.seed(15) # makes rnorm reproducible par(mfrow = c(2,1)) x <-

Hi, I'd like to use Latin Hypercube Sampling (LHC) in the the context of uncertainty / sensitivity analysis of a complex model with approximately 10 input variables. With the LHC approach I'd like to generate parameter combinations for my model input variables. Therefore I came across an simple example here on the mailing list (

Hola, LLevo buscando desde hace tiempo como hacer el Goodness of fit test en R. Es decir, me explico, intento hacer una cosa parecida que se hace en Minitab, por ejemplo, yo tengo un conjunto de datos, y lo que quiero es sabes que tipo de distibución es, en minitab se hace un histograma para ver si se ajusta bien o no a la campana de Gauss, luego vemos si aproximar la distribución de la muestra

Hi All, I am in the annoying position of having to present some data to someone who seems to be somewhat less than numerate. I need to label the y-axes of a multhist with the y-axis labeled not as counts but as percentage of a population. Plotting the standard histogram is in a way fine, all I need is to: -- have a left-handside y-axis labels for pop 1 and a right-handside y-axis labels for pop2

I am just trying to teach myself how to use the mle function in R because it is much better than what is provided in MATLAB. I am following tutorial material from the internet, however, it gives the following errors, does anybody know what is happening to cause such errors, or does anybody know any better tutorial material on this particular subject. >

Hi, Try either: set.seed(28) stats1<- as.data.frame(matrix(rnorm(5*10000),ncol=5)) pdf(paste("test",1,".pdf",sep="")) par(mfrow=c(2,1)) lst1<- lapply(names(stats1),function(i) {hist(stats1[,i],100,col="lightblue",main=paste0("Histogram of ",i),xlab=i );qqnorm(stats1[,i])}) dev.off() #or

Full_Name: Adrian Custer Version: 0.90.0 OS: Linux on Thinkpad (pentium) and desktop (K6) Submission from: (NULL) (128.32.251.234) When I create a boxplot, and then try to overlay a lowess fit or just the points, the points do not appear in the highest level and the lowess curve does not reach the highest level. However, if I add one to each of the models, the problem is solved. I tried this

Hi Hugo, I've been able to replicate your bug, including for other distributions (runif, rexp, rgamma, etc) which shouldn't be surprising since they're probably all drawing from the same pseudo-random number generator. ?Interestingly, it does not seem to depend on the choice of seed, I am not sure why that is the case. I'll point out first of all that the R-devel mailing list is

Hi Hugo, I've been able to replicate your bug, including for other distributions (runif, rexp, rgamma, etc) which shouldn't be surprising since they're probably all drawing from the same pseudo-random number generator. ?Interestingly, it does not seem to depend on the choice of seed, I am not sure why that is the case. I'll point out first of all that the R-devel mailing list is

Full_Name: Adrian Custer Version: 1.01 OS: gnu-linux (mandrake 7.1beta) Submission from: (NULL) (128.32.251.234) The character "*" displays too high in comparison to the graph space and other characters in the plot. Consider: > a_-0.15 > plot(c(0,100),c(0,1)) > for (i in 1:100){ points(i,(1-exp(a*i)),pch=".",col="blue")} > for (i in 1:100){

I am trying to run Samba 3.0.24-11.fc6 on Fedora Core 6. I've run across two problems that are keeping me from being successful. 1. IPTables. Even when I configure it through Webmin to not run at startup, it does anyway. I consider this a minor problem, at least for now, although I wouldn't mind some input on how to configure it to allow Samba. 2. Access to shares. We are running

Running FC6 and Samba 3.0.24-11.fc6. Workstation is XP Pro Initially, after disabling SELinux and turning off Iptables, I could access the shares that were set up. Came in this morning ready to start copying files and can't get in. Prompts me for a username and password which is rejected. I know it's not my workstation because I have another FC6/Samba server running and have no problem

Dear list, After running for a while, it crashes and gives the following error message: can anybody suggest how to deal with this? Error in if (ratio0[i] < log(runif(1))) { : missing value where TRUE/FALSE needed ################### original program ######## p2 <- function (Nsim=1000){ x<- c(0.301,0,-0.301,-0.602,-0.903,-1.208, -1.309,-1.807,-2.108,-2.71) # logdose

dear list, i have a problem using the par function. in one graphic device i want to have two plots so i tried to use par(mfror=c(1,2)). of course it worked out, but the height now is twice the length of the width for each single plot. what i actually wanted is something like par(mfrow=c(2,2)) where only the top (or bottom) two plots are drawn with entire length proportions of 1 to 2 ( height 1,