search for: rcauchy

...and for the normal, but not for the cauchy. Any idea what I'm doing wrong? Thanks x <- rchisq(10000, df = 4) hist(x, freq = FALSE, breaks=100) curve(dchisq(x, df = 4), col = 2, add = TRUE) x <- rnorm(10000) hist(x, freq = FALSE, breaks=100) curve(dnorm(x), col = 2, add = TRUE) x <- rcauchy(10000) hist(x, freq = FALSE, breaks=100) curve(dcauchy(x), col = 2, add = TRUE) [[alternative HTML version deleted]]

Hi, I have the following R-code: pdf(file="test.pdf") hist(rcauchy(100)) dev.off() How can I set the size of the histogram in my pdf? I tried to do with par(...) but I didn't find the right parameter yet. Thanks, Max -- View this message in context: http://www.nabble.com/set-size-of-a-plot-tp22477083p22477083.html Sent from the R help mailing list archive a...

Hi R-users, Recently, I use maxNR function to find maximizer. I have error appears as follows Error in maxNRCompute(fn = fn, grad = grad, hess = hess, start = start, : NA in the initial gradient My code is mu=2 s=1 n=300 library(maxLik) set.seed(1004) x<-rcauchy(n,mu,s) loglik<-function(mu) { log(prod(dcauchy(x,mu,s))) } maxNR(loglik,start=median(x))$estimate Does anyone know how to solve this problem? Thanks, Kate [[alternative HTML version deleted]]

...lly in a way as simple as getting random numbers from a normal distribution with 'rnorm'. Is there any package to do this? Speed is an issue in this application. Or, as the Voigt distribution is a convolution of a Gaussian and a Lorentzian, can I simply combine random numbers from rnorm and rcauchy in some way (I'm not sure about the maths of this)? Best wishes, David. <DIV><FONT size="1" color="gray">This e-mail and any attachments may contain confidential, copyright and or privileged material, and are for the use of the intended addressee only. If you ar...

...ere's anything in R that could help me do that. Let's suppose I have a density function of a random variable, for example f(x) = (x^3)/4 0 < x < 2 and I would like to simulate it. For the common distributions (exponencial, gamma, cauchy) there are the r-functions (rgamma, rexp, runif, rcauchy, and so on).. But when the variable I want to simulate is not one of those, how should I procede? I read some references on the subject and saw that there are some algorithms that can do that, but I just wonder if there is any implemented in R? Thank you, --

...e is new to R, as there are too many hits for "Random" and they dive too deeply. My answer to this FAQ would be: In order to generate random numbers, there are "r" methods for every distribution. Distributions include Uniform (runif), Binomial (rbinom), Normal (rnorm), Cauchy (rcauchy), Chi Square (rchisq), Exponential (rexp), Gamma (rgamma), and many more. By the way, the R program itself is really badly in need of an "appropos()" to go alongside "help()". Perhaps "help()" can be renamed to "hunt()" or "huntAndPeck()".

Hello, It might be an easy question but if you have many variables to fit in the lm function, how do you take all without specifying var1+var2+...+var2100 in the terms parameter in response ~ terms? Cheers, Carol

...beta.o qbeta.o rbeta.o dunif.o punif.o qunif.o runif.o dnorm.o pnorm.o qnorm.o rnorm.o dlnorm.o plnorm.o qlnorm.o rlnorm.o df.o pf.o qf.o rf.o dnf.o dt.o pt.o qt.o rt.o dnt.o dchisq.o pchisq.o qchisq.o rchisq.o rnchisq.o dbinom.o pbinom.o qbinom.o rbinom.o rmultinom.o dcauchy.o pcauchy.o qcauchy.o rcauchy.o dexp.o pexp.o qexp.o rexp.o dgeom.o pgeom.o qgeom.o rgeom.o dhyper.o phyper.o qhyper.o rhyper.o dnbinom.o pnbinom.o qnbinom.o rnbinom.o dpois.o ppois.o qpois.o rpois.o dweibull.o pweibull.o qweibull.o rweibull.o dlogis.o plogis.o qlogis.o rlogis.o dnchisq.o pnchisq.o qnchisq.o dnbeta.o pnbeta....

I want to estimate the center of a distribution with lots of outliers in one tail, and thought I would use a function such as S-plus's location.m() with psi.fun=bisquare (as per MASS 3 p. 131). However, R seems not have such a function, so my questions are: 1) Is there an R equivalent to location.m()? 2) Would huber() give me results that are similar (i.e., close enough)? Thanks.

Folks, I wonder if anything could be said about the distribution of a random variate x, where x = N(0,1)/N(0,1) Obviously x is pathological because it could be 0/0. If we exclude this point, so the set is {x/(0/0)}, does x have a well defined distribution? or does it exist a distribution that approximates x. (The case could be generalized of course to N(mu1, sigma1)/N(mu2, sigma2) and one

...-I../include -c pbinom.c -o pbinom.o gcc -g -I../include -c qbinom.c -o qbinom.o gcc -g -I../include -c rbinom.c -o rbinom.o gcc -g -I../include -c dcauchy.c -o dcauchy.o gcc -g -I../include -c pcauchy.c -o pcauchy.o gcc -g -I../include -c qcauchy.c -o qcauchy.o gcc -g -I../include -c rcauchy.c -o rcauchy.o gcc -g -I../include -c dchisq.c -o dchisq.o gcc -g -I../include -c pchisq.c -o pchisq.o gcc -g -I../include -c qchisq.c -o qchisq.o gcc -g -I../include -c rchisq.c -o rchisq.o gcc -g -I../include -c dexp.c -o dexp.o gcc -g -I../include -c pexp.c -o pexp.o gcc -g -I../incl...

...seriesin: a numeric vector numlags: an integer, can be negative panelid: an optional panel id @DETAILS None @SEEALSO leadseries, chgseries, pchgseries, compoundseries @EXAMPLES x <- rnorm(10) xlag <- lagseries(x,2) lm( x ~ xlag ) d <- data.frame( x <- c( rnorm(20), runif(30), rcauchy(40) ), who= c( rep("firm1",20), rep("firm2", 30), rep("firm3",40)), year= c( 1961:1980, 1971:2000, 1971:2010 ) ) lagd <- data.frame( x=lagseries(d$x, panelid=who), who=d$who, year=lagseries(d$x, panelid=who) ) ', test = '...

Hello, Can you get eigenvalues in addition to eigevectors using prcomp? If so how? I am unable to use princomp due to small sample sizes. Thank you in advance for your help! Rebecca Young -- Rebecca Young Graduate Student Ecology & Evolutionary Biology, Badyaev Lab University of Arizona 1041 E Lowell Tucson, AZ 85721-0088 Office: 425BSW rlyoung at email.arizona.edu (520) 621-4005

Hi all, I am trying hard to do the following and have already spent a few hours in vain: I wanted to do the scatter plot. But given the high dispersion on those dots, I would like to bin the x-axis and then for each bin of the x-axis, plot the quantiles of the y-values of the data points in each bin: 1. Uniform bin size on the x-axis; 2. Equal number of observations in each bin; How to do

Dear R-users, Following is part of my data, where slide has 36 levels and block 48 levels. I have done boxplot for each slide on the same graph. There are outliers for each slide and I tried to use indentify functtion to identify outliers in such a way that when I click on an outlier or point, the points will be labelled by either their block or ID or by both but without success. How can I make

I would be very grateful for any help from members of this list for what might be a simple problem... We are trying to simulate the behaviour of a clinical measurement in a series of computer experiments. This is simple enough to do in R if we assume the measurements to be Gaussian, but their empirical distribution has a much higher peak at the mean and the distribution has much longer tails.

...g -I../include -I../../src/include -c dcauchy.c -o dcauchy.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pcauchy.c -o pcauchy.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qcauchy.c -o qcauchy.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rcauchy.c -o rcauchy.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dexp.c -o dexp.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pexp.c -o pexp.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qexp.c -o qexp.o gcc -ieee_with_inexact -g -I../include -I.....

...de -c qbinom.c cc -ieee_with_inexact -g -I../../src/include -c rbinom.c cc -ieee_with_inexact -g -I../../src/include -c dcauchy.c cc -ieee_with_inexact -g -I../../src/include -c pcauchy.c cc -ieee_with_inexact -g -I../../src/include -c qcauchy.c cc -ieee_with_inexact -g -I../../src/include -c rcauchy.c cc -ieee_with_inexact -g -I../../src/include -c dexp.c cc -ieee_with_inexact -g -I../../src/include -c pexp.c cc -ieee_with_inexact -g -I../../src/include -c qexp.c cc -ieee_with_inexact -g -I../../src/include -c rexp.c cc -ieee_with_inexact -g -I../../src/include -c dgeom.c cc -ieee_with_in...

...de -I../../src/include -c dcauchy.c = -o = dcauchy.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pcauchy.c = -o = pcauchy.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qcauchy.c = -o = qcauchy.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rcauchy.c = -o = rcauchy.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c dexp.c -o = dexp.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c pexp.c -o = pexp.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qexp.c -o = qexp.o gcc -ieee_with_inexact -g -I../in...

There is a problem with `rpois'. It does seem to take care about the order of the arguments. This is an example: > rpois(n=1,lambda=2) [1] 3 > rpois(lambda=2,n=1) [1] 2 0 It obviously uses the first argument as the number of samples to be drawn, which is wrong. I used Version 0.49 Beta (April 23, 1997). Fredrik