search for: dcauchy

...ributions for variance parameters in hierarchical models" in Bayesian Analysis, 3, 515-533. I am only interested, however, in recreating the portion of the graph for the overlain prior density for the half-Cauchy with scale 25 and not the posterior distribution. However, when I try: curve(dcauchy, from=0, to=200, location=0, scale=25) the probabilities for the half-Cauchy values seem to approach zero almost immediately after 0 whereas in Gelman 2006 the tail appears much fatter giving non-zero probabilities out to 100. I am interested in replicating this because I want to use half-Cauch...

...ntly, 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]]

...term. x <- c(1.808,1.799,1.179,0.574,3.727,0.946,3.719,1.566,3.596,3.253) y <- c(1.816,1.281,-1.382,0.573,3.793,0.935,1.775,1.474,3.679,3.889) fn = function(x,a=0,b=1){ a+b*x } sample.ab <-function(x,y,a,b,s,da,db){ bstar = runif(1,b-db,b+db) astar = runif(1,a-da,a+da) logalpha = sum(dcauchy(y,location=fn(x,astar,bstar),scale=s,log=T) - dcauchy(y,location=fn(x,a,b),scale=s,log=T)) logu = log(runif(1,0,1)) acc = (logu < logalpha) b = acc*bstar + (1-acc)*b a = acc*astar + (1-acc)*a list(b=b,a=a,acc=acc) } samples = function(x,y,a,b,s,ds){ sstar = runif(1,s-ds,s+ds) while(sstar...

...hat 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]]

...er/SDKs/MacOSX10.4u.sdk -mmacosx- version-min=10.4 -std=gnu99 -dynamiclib -Wl,- headerpad_max_install_names -mmacosx-version-min=10.4 -undefined dynamic_lookup -single_module -multiply_defined suppress -L/usr/local/ lib -o pnmath.so bd0.o beta.o chebyshev.o choose.o d1mach.o dbeta.o dbinom.o dcauchy.o dchisq.o dexp.o df.o dgamma.o dgeom.o dhyper.o dlnorm.o dlogis.o dnbeta.o dnbinom.o dnchisq.o dnf.o dnorm.o dnt.o dpois.o dt.o dunif.o dweibull.o fmax2.o fmin2.o ftrunc.o gamma.o gamma_cody.o i1mach.o imax2.o imin2.o lbeta.o lgamma.o lgammacor.o mlutils.o pbeta.o pbinom.o pcauchy.o pchisq...

Hi everybody, I''m using R on a 64-bit Ubuntu 9.04 (Jaunty). I prefer to install R packages from source, even if they are available in Synaptic. The problem is that I can''t install/compile kernlab. Everything works fine until it gets to the lazy loading part: ** preparing package for lazy loading Creating a new generic function for "terms" in "kernlab"

....o qgamma.o rgamma.o dbeta.o pbeta.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 pnchis...

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

....o gcc -g -I../include -c qbeta.c -o qbeta.o gcc -g -I../include -c rbeta.c -o rbeta.o gcc -g -I../include -c dbinom.c -o dbinom.o gcc -g -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 gc...

...=4) > mtext("log(dbinom(k))", col="red", adj=1) > > > > cleanEx(); ..nameEx <- "Cauchy" > > ### * Cauchy > > flush(stderr()); flush(stdout()) > > ### Name: Cauchy > ### Title: The Cauchy Distribution > ### Aliases: Cauchy dcauchy pcauchy qcauchy rcauchy > ### Keywords: distribution > > ### ** Examples > > dcauchy(-1:4) [1] 0.15915494 0.31830989 0.15915494 0.06366198 0.03183099 0.01872411 > > > > cleanEx(); ..nameEx <- "ChickWeight" > > ### * ChickWeight > > flush(st...

...xact -g -I../include -I../../src/include -c pbinom.c -o pbinom.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qbinom.c -o qbinom.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rbinom.c -o rbinom.o gcc -ieee_with_inexact -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...

...clude -c rchisq.c cc -ieee_with_inexact -g -I../../src/include -c dbinom.c cc -ieee_with_inexact -g -I../../src/include -c pbinom.c cc -ieee_with_inexact -g -I../../src/include -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...

.../include -I../../src/include -c pbinom.c -= o = pbinom.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c qbinom.c -= o = qbinom.o gcc -ieee_with_inexact -g -I../include -I../../src/include -c rbinom.c -= o = rbinom.o gcc -ieee_with_inexact -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 -i...

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

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

...-I../../src/include -c dbinom.c -o dbinom.o gcc -g -O2 -I../include -I../../src/include -c pbinom.c -o pbinom.o gcc -g -O2 -I../include -I../../src/include -c qbinom.c -o qbinom.o gcc -g -O2 -I../include -I../../src/include -c rbinom.c -o rbinom.o gcc -g -O2 -I../include -I../../src/include -c dcauchy.c -o dcauchy.o gcc -g -O2 -I../include -I../../src/include -c pcauchy.c -o pcauchy.o gcc -g -O2 -I../include -I../../src/include -c qcauchy.c -o qcauchy.o gcc -g -O2 -I../include -I../../src/include -c rcauchy.c -o rcauchy.o gcc -g -O2 -I../include -I../../src/include -c dexp.c -o dexp.o gcc -g...

...-I../../src/include -c dbinom.c -o dbinom.o gcc -g -O2 -I../include -I../../src/include -c pbinom.c -o pbinom.o gcc -g -O2 -I../include -I../../src/include -c qbinom.c -o qbinom.o gcc -g -O2 -I../include -I../../src/include -c rbinom.c -o rbinom.o gcc -g -O2 -I../include -I../../src/include -c dcauchy.c -o dcauchy.o gcc -g -O2 -I../include -I../../src/include -c pcauchy.c -o pcauchy.o gcc -g -O2 -I../include -I../../src/include -c qcauchy.c -o qcauchy.o gcc -g -O2 -I../include -I../../src/include -c rcauchy.c -o rcauchy.o gcc -g -O2 -I../include -I../../src/include -c dexp.c -o dexp.o gcc -g...

...rom pt.c making qt.d from qt.c making rt.d from rt.c making dchisq.d from dchisq.c making pchisq.d from pchisq.c making qchisq.d from qchisq.c making rchisq.d from rchisq.c making dbinom.d from dbinom.c making pbinom.d from pbinom.c making qbinom.d from qbinom.c making rbinom.d from rbinom.c making dcauchy.d from dcauchy.c making pcauchy.d from pcauchy.c making qcauchy.d from qcauchy.c making rcauchy.d from rcauchy.c making dexp.d from dexp.c making pexp.d from pexp.c making qexp.d from qexp.c making rexp.d from rexp.c making dgeom.d from dgeom.c making pgeom.d from pgeom.c making qgeom.d from qgeom....

...a.o qgamma.o rgamma.o dbeta.o pbeta.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...

...rom pt.c making qt.d from qt.c making rt.d from rt.c making dchisq.d from dchisq.c making pchisq.d from pchisq.c making qchisq.d from qchisq.c making rchisq.d from rchisq.c making dbinom.d from dbinom.c making pbinom.d from pbinom.c making qbinom.d from qbinom.c making rbinom.d from rbinom.c making dcauchy.d from dcauchy.c making pcauchy.d from pcauchy.c making qcauchy.d from qcauchy.c making rcauchy.d from rcauchy.c making dexp.d from dexp.c making pexp.d from pexp.c making qexp.d from qexp.c making rexp.d from rexp.c making dgeom.d from dgeom.c making pgeom.d from pgeom.c making qgeom.d from qgeom....