similar to: random number generation

Hello, I use newly R! I'd like to plot several data set together in one output window! How can I do that? Best regards Hadi --------------------------------- Never miss an email again! [[alternative HTML version deleted]]

Hi, Still playing around with Samba 2.0.5a running on HP-UX 10.20. Attempted to startup the smbd and nmbd daemon during system startup. smbd starts up but not the nmbd. Is there any hints. thanks in advance. best

Is this an interface bug? Using RGUI for windows I run into a "Not Responding" process (I "smartly" coded an infinite loop, yaiks!), I hit esc and the interpreter was stopped and I recovered the console functionality but the caption on the R icon in my windows taskbar (the individual icon shown for every software currently running in the session) was not updated so the

I was trying to generate random numbers with a gamma distribution. In R the function is: rgamma(n, shape, rate = 1, scale = 1/rate). My question is that if X~gamma(alpha, beta) and I want to generate one random number where do I plug alpha and beta in rgamma? and, what is the meaning and use of rate? Thanks for your attention, Jorge [[alternative HTML version deleted]]

Hi, I have the following vector which is created from 3 distinct distribution (three components) of gamma: x=c(rgamma(30,shape=.2,scale=14),rgamma(30,shape=12,scale=10),rgamma(30,shape=5,scale=6)) I want to plot the density curve of X, in a way that it shows a distinct 3 curves that represent each component. How can I do that? I tried this but doesn't work: lines(density(x)) Please

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 R Users This is a code I wrote and just want to confirm if the first 1000 values are raw gamma (z) and the next 1000 values are transformed gamma (k) or not. As I get 2000 rows once I import into excel, the p - values beyond 1000 dont look that good, they are very high. -- sink("a1.txt"); for (i in 1:1000) { x<-rgamma(10, 2.5, scale = 10) y<-rgamma(10, 2.5, scale = 10)

Hi all, I have this project that requires me to use 8 PSTN lines and possible more. I was thinking 2 TDM cards with FXO modules. The I got to read the "Qs about FXO/FXS cards" thread and that scared me. Can anybody recommend anything that is known to work ok with no mysterious problems? I was thinking OpenSwitch12 cards. What do you guys think? Any help is appreciated. Regards, Hadi

Hi, Has anyone encountered the problem of rgamma function in C? The following simplified program always dies for me, and I wonder if anyone can tell me the reason. #include <Rmath.h> #include <time.h> #include <Rinternals.h> SEXP generateGamma () { srand(time(NULL)); return (rgamma(5000,1)); } Has anyone encountered a similar problem before? Is there another way

Dear R people, This is a possibly silly question, but the rgamma function takes the shape and scale arguments and simulates gamma rvs corresponding to those values, right? But the shape and scale parameters have to be positive, right? However, rgamma quite happily returns to me values for negative values of shape and scale, and in some cases returns negative values eg. > rgamma(1, 1, -1) [1]

Hi all, I have the following x1<-paste("A", 1:6, sep = "") x2<- round(rgamma(6,2,1)) x3<-paste("B", 1:6, sep = "") x4<- round(rgamma(6,2,1)) data1 <- data.frame(x1,x2,x3,x4) I would like to get data2 <- c(A1=4, A2=1, A3=0,...) Is there any standard for such a case? Thank you very much in advance, Diego

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

Esta pregunta nos la hicimos muchos, alguien aproxima una respuesta a través de esta publicación: http://www.scielo.org.ar/pdf/ecoaus/v18n2/v18n2a07.pdf Saludos cordiales Teresa -- Ing. Agr. M.Sc.Teresa Boca Instituto de Clima y Agua INTA Castelar. Las Cabañas y Los Reseros s/n° 1712 Castelar Pcia de Buenos Aires, Argentina Cel: 01115-3572-9584 Tel: (54-11)-4621-0125/1684 Fax:(54-11)-4621-5663

Hi, I have asterisk installed on centos with phpagi. Also I have PRI card connect to it. it's possible to show the sip number when calling from sip number to external number thru the PRI, instead of showing the PRI number show the sip number ? Regards -Hadi.Salem

Dear All Can you please do me favor and let me know how can I convert *.wav files into 32 bit 44 KHz ? Please be informed that I have specific sound files in *.wav format that I converted them into *.gsm format with the aid of the following command : #sox FR00003.wav FR00003.gsm It got through but the voice quality is poor . I need to convert the original *.wav sound files (their file attribute is

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

Dear All, I would like to ask please if someone has a AGI script that select a value from a database and dial this value as a destination number Regards -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20090220/e2aa530c/attachment.htm

Hello everyone, In plotting several graphics, I'd like to be able to plot a blank plot as in: par(mfrow=c(2,1)) plot(BLANK) hist(rgamma(100000,6463.7,scale=0.015471),xlim=c(0,120),main="Emergence") I realize screen allows me to do this, but I figure the functionality must be there. Is there an equivalent to plot(BLANK)? thanks, adrian PS. please reply direct as I'm not

Hi: I am writing a C program which need a gamma random number generator. I download the source file of R and compile, make it myself. There is a "rgamma.c" function in the installing directory of R("/home/zhliu/Backup/R-2.0.1/src/nmath/rgamma.c"). My question is how to call this function in my own program which is in another directory. I can not copy this "rgamma.c"

Dear all, I am using the bigglm package to fit a few GLM's to a large dataset (3 million rows, 6 columns). While trying to fit a Poisson GLM I noticed that the coefficient estimates were very different from what I obtained when estimating the model on a smaller dataset using glm(), I wrote a very basic toy example to compare the results of bigglm() against a glm() call. Consider the