Displaying 20 results from an estimated 38 matches for "distr".

Did you mean:
distro

2011 Sep 21

2

problem with function "Truncate" in package "distr"

Hello all,
Can someone tell me why the following mixture of two log-normal
distributions does not get truncated? What puzzles me is that the
function works almost always, but for certain combinations (like the
one below), it does not.
# R code example
library(distr)
mix<-UnivarMixingDistribution(Lnorm(3.2,0.5),Lnorm(5.4,0.6),mixCoeff=c(0.3,0.7))
mix.trunc<-Truncate(mix,l...

2008 May 23

1

Evaluation "conflict" in combination of replicate() and rexp()as variable inside a function

...enerating functions inside a self-defined "wrapper" function I encounterd
a puzzling behaviour.
The following are intentionally simple (and rather nonsense-) examples to
isolate the relevant aspects. Please, note the seemingly "inconsistent"
behaviour for the second call of rdistr() (with distr = rexp) for which I
have not found any explanation (yet):
> rdistr <- function( distr, ...) replicate( 1, distr( n = 1, ...))
> rdistr( distr = rnorm)
[1] -0.8889223
> rdistr( distr = rexp)
[1] NaN
Warning message:
In distr(n = 1, ...) : NAs produced
> rdistr( dist...

2011 Oct 06

2

apply and functions with many arguments

Dear all,
I would like to use the following function
fitdist(data, distr, method=c("mle", "mme", "qme", "mge"),
start=NULL, fix.arg=NULL, ...)
for many different distr values like distr=c("norm","lnorm","pois") (just a small example)
and take back into a list the parameter name which is what is insi...

2017 Jun 19

0

quantreg::rq.fit.hogg crashing at random

...em?")
list(coefficients = coefficients, nit = z$it.count, flag = z$info)
}
#----- Here is the main program that crashes
library(MASS)
library(quantreg)
M=1e3; n=1e3; p=15; type=8; K=10 # no. of replications, sample size,
dimension of beta, type of quantile estimator, number of quantiles
distributions=c("norm","t1","t2","t3","t5","logistic","
exponential","Weibull")
method.wqr="fn" # interior point method; if "fn" then roughly matches with
method.cqr="ip"
# Create the covariance...

2017 Nov 22

2

function pointers?

We have a project that calls for the creation of a list of many
distribution objects. Distributions can be of various types, with
various parameters, but we ran into some problems. I started testing
on a simple list of rnorm-based objects.
I was a little surprised at the RAM storage requirements, here's an example:
N <- 10000
closureList <- vector("...

2017 Nov 15

1

ks.test() with 2 samples vs. 1 sample an distr. function

In the first example you are performing a one-sample test against a continuous cumulative distribution (in this case a normal distribution). In the second case you are performing a two-sample test. You drew your values for x non-randomly by specifying fixed intervals along a normal distribution, but ks.test() just sees that you have provided two samples, not one sample and values along a cumu...

2017 Nov 15

1

ks.test() with 2 samples vs. 1 sample an distr. function

...the example given on the German wikipedia page.
xi <- c(9.41,9.92,11.55,11.6,11.73,12,12.06,13.3)
I get the right results when I calculate: ks.test(xi,pnorm,11,1)
Now the question: shouldn't I obtain the same or a very similar result if I commpare the sample and a calculated sample from the distribution?
p<- c(0.125, 0.250, 0.375, 0.500, 0.625, 0.750, 0.875, 0.9999)
x <- qnorm(p,11,1)
ks.test(xi,x)
Why don't I?
Thanks for helping me!
Tonja

2017 Aug 27

1

Fwd: Find maxima of a function

...will all depend on the user.
Then on the result of kmeans, I wanted to find the maxima so that I can use
Normalized probability formula further for my logic and in order to do that
I used Newton's method and that was implemented in the optimr package.
I was getting the wrong result so I used distr package in order to check
the problem and I figured out I was getting the wrong maxima. So I can to
the conclusion which I have posted as my query.
PS: I am able to use distr package without any problem.
And I want to use dnorm because I have a Gaussian mixture model, so I did
not want to use rno...

2017 Aug 26

1

Find maxima of a function

...; wrote:
>
> Hi,
>
> Thanks for your mail, and time
>
> It is not working for some arguments, when mean value is like >6.
>
>
> case
>
> mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0.
> 1134463,0.2074027)
>
> rv <-UnivarMixingDistribution(Norm(486.4255, 53.24133),
>
> Norm(664.0713, 3.674773),
>
> Norm(669.0484, 4.101381),
>
> Norm(677.1753, 4.869985),
>
> Norm(683.2635, 7.288175),
&g...

2017 Aug 27

1

Fwd: Find maxima of a function

...eeks on a very easy issue (it is very long time). But if you don?t want to loose anymore, please, isolate your question from the rest of the work, create minimal reproduciple example, and let?s focus on the problematic part. I assure you that your problem will be solved faster.
I could install the distr package at the end, but I?m getting another error while running your code. But I don?t believe the question is related to this package.
Let me explain about extremum points although most of you know about it. Let?s assume we have a y(x) function. An extremum (max/min) point is defined as where...

2004 Oct 17

11

Errors while compiling packages with namespace?

Hello,
I try to set up namespaces for packages. It is fine for several of them,
except one whose compilation fails (under Windows XP & R 2.0.0):
---------- Making package svViews ------------
adding build stamp to DESCRIPTION
installing NAMESPACE file and metadata
Error in parse(file, n, text, prompt) : syntax error on line 21
Execution halted
make[2]: *** [nmspace] Error 1
make[1]: ***

2017 Aug 27

1

Fwd: Find maxima of a function

...very easy issue (it is very long time). But if you don?t want to loose anymore, please, isolate your question from the rest of the work, create minimal reproduciple example, and let?s focus on the problematic part. I assure you that your problem will be solved faster.
>
> I could install the distr package at the end, but I?m getting another error while running your code. But I don?t believe the question is related to this package.
>
> Let me explain about extremum points although most of you know about it. Let?s assume we have a y(x) function. An extremum (max/min) point is define...

2004 Oct 13

4

Cleaning up R-2.0.0 installation on GNU/Linux ?

...ckages in library ''/usr/local/lib/R/site-library'':
accuracy ** No title available (pre-2.0.0 install?) **
ade4 ** No title available (pre-2.0.0 install?) **
debug ** No title available (pre-2.0.0 install?) **
distr ** No title available (pre-2.0.0 install?) **
eha Event History Analysis.
ellipse ** No title available (pre-2.0.0 install?) **
faraway ** No title available (pre-2.0.0 install?) **
fortunes ** No title available (...

2009 Dec 11

4

how can generate from trunceted gamma distribution in R ?

Hi, all
How can generate a sample from truncated inverse gamma distribution in R?
thanks

2006 Apr 26

2

cdf of weibull distribution

Hi,
I have a data set which is assumed to follow weibull distr''. How can I find of cdf for this data. For example, for normal data I used (package - lmomco)
>cdfnor(15,parnor(lmom.ub(c(df$V1))))
Also, lmomco package does not have functions for finding cdf for some of the distributions like lognormal. Is there any other package, which c...

2011 Aug 01

3

Beta fit returns NaNs

Hi,
sorry for repeating the question but this is kind of important to me and i
don''t know whom should i ask.
So as noted before when I do a parameter fit to the beta distr i get:
fitdist(vectNorm,"beta");
Fitting of the distribution '' beta '' by maximum likelihood
Parameters:
estimate Std. Error
shape1 2.148779 0.1458042
shape2 810.067515 61.8608126
Warning messages:
1: In dbeta(x, shape1, shape2, log) : NaNs produced
2: In d...

2005 Feb 24

1

Density of the Multivariate T Distribution

Hi,
I am looking for an efficient way to compute the
values of the density function of a multivariate T
distribution - something like "dmvnorm", but for T
distr. Does this exist somewhere?
Many thanks,
Jan Bulla
Goettingen University

2005 Jun 15

1

Chi square convolution?

Hi,
I want to determine the confidence interval on the sum of two sigma''s.
Is there an easy way to do this in R? I guess I have to use some sort of
chisquare convolution algorithm???
Thanx,
Roy
--
The information contained in this communication and any atta...{{dropped}}

2011 Nov 14

0

Samba, RDP

...= recycle full_audit
recycle:repository = /hard1/basket
recycle:exclude = *.tmp | *.TMP | ~$* | ~WRL*
recycle:keeptree = yes
recycle:versions = yes
veto files = /*.mp3/*.lnk/*.avi/
mix protocol = SMB2
[????????]
browseable = yes
writable = yes
path = /hard1/total
[distr]
browseable = no
writable = yes
path = /hard1/distr
include = /autostart/samba/%U
Everything works perfectly, as long as users went to the file server (FS) through windows XP and every one from his computer.
But needed to create a separate server to perform different tasks, it...

2017 Aug 26

2

Fwd: Find maxima of a function

Hi,
Thanks for your mail, and time
It is not working for some arguments, when mean value is like >6.
case
mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0.
1134463,0.2074027)
rv <-UnivarMixingDistribution(Norm(486.4255, 53.24133),
Norm(664.0713, 3.674773),
Norm(669.0484, 4.101381),
Norm(677.1753, 4.869985),
Norm(683.2635, 7.288175),
Norm(72...