search for: distr

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