Displaying 20 results from an estimated 6000 matches similar to: "Generating by inverting function"
2012 Aug 29
2
Estimation parameters of lognormal censored data
Hi, I am trying to get the maximum likelihood estimator for lognormal distribution with censored data;when we have left, interval and right censord. I built my code in R, by writing the deriving of log likelihood function and using newton raphson method but my estimators were too high " overestimation", where the values exceed the 1000 in some runing of my code.
is there any one can
2012 Aug 31
3
fitting lognormal censored data
Hi ,
I am trying to get some estimator based on lognormal distribution when we have left,interval, and right censored data. Since, there is now avalible pakage in R can help me in this, I had to write my own code using Newton Raphson method which requires first and second derivative of log likelihood but my problem after runing the code is the estimators were too high. with this email ,I provide
2006 Apr 26
1
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 can handle these distributions?
2011 Mar 15
3
fitting a distribution to a ecdf plot
Dear all,
I need to plot an cumulative distribution plot of a variable and then to fit a distribution to that, probably a weibull or lognormal.
I have plotted the ecdf as
> plot(ecdf(x))
but I haven't managed to fit the distribution. I have as well attached the data.
I would appreciate if you could help me on that.
Thank you.
Kind regards
Maria
-------------- next part --------------
2006 Jul 07
6
parametric proportional hazard regression
Dear all,
I am trying to find a suitable R-function for
parametric proportional hazard regressions. The
package survival contains the coxph() function which
performs a Cox regression which leaves the base hazard
unspecified, i.e. it is a semi-parametric method. The
package Design contains the function pphsm() which is
good for parametric proportional hazard regressions
when the underlying base
2008 Jul 06
4
Method for checking automatically which distribtions fits a data
Hi,
Suppose I have a vector of data.
Is there a method in R to help us automatically
suggest which distributions fits to that data
(e.g. normal, gamma, multinomial etc) ?
- Gundala Viswanath
Jakarta - Indonesia
2008 Feb 22
1
fitting a lognormal distribution using cumulative probabilities
Dear all,
I'm trying to estimate the parameters of a lognormal distribution fitted
from some data.
The tricky thing is that my data represent the time at which I recorded
certain events. However, in many cases I don't really know when the event
happened. I' only know the time at which I recorded it as already happened.
Therefore I want to fit the lognormal from the cumulative
2010 May 26
2
Survival analysis extrapolation
Dear all,
I'm trying to fit a curve to some 1 year failure-time data, so that I can
extrapolate and predict failure rates up to 3 years. The data is in the
general form:
Treatment Time Status
Treatment A 28 0
Treatment B 28 0
Treatment B 28 0
Treatment A 28
2008 May 29
1
Bimodal Distribution
Hello R Users,
I am doing a Latin Hypercube type simulation. I have found the
improvedLHS function and have used it to generate a bunch of properly
distributed uniform probabilities. Now I am using functions like qlnorm
to transform that into the appropriately lognormal or triangularly
distributed parameters for my modes. However I have a parameter which I
believe is bimodally distributed,
2011 Jun 02
4
[Plea to the R Gods] Theoretical and Empirical CDFs
http://r.789695.n4.nabble.com/file/n3567636/ecdfs.jpg ecdfs.jpg
http://r.789695.n4.nabble.com/file/n3567636/ecdf_curve.gif ecdf_curve.gif
Hello,
I have generated a plot of two empirical CDFs (attachment 1). As a result,
they are stepwise when plotted. The following code was used:
> plot(ecdf(mut), do.points=FALSE, verticals=TRUE, xlim=range(mut, non),
> col="red")
>
2006 Sep 21
1
survival function with a Weibull dist
Hi
I am using R to fit a survival function to my data
(with a weibull distribution).
Data: Survival of individuals in relation to 4
treatments ('a','b','c','g')
syntax:
---- > survreg(Surv(date2)~males2, dist='weibull')
But I have some problems interpreting the outcome and
getting the parameters for each curve.
--------- Value Std.
2004 Feb 03
4
how to change one of the axis to normal probability scale
Hi,
I would like to plot a graph with one axis in log scale and the other in
normal probability scale using R. I cannot change the axis scale to
probability scale. What can be a solution? Thanks.
CT
1999 Dec 09
1
nlm() problem or MLE problem?
I am trying to do a MLE fit of the weibull to some data, which I attach.
fitweibull<-function()
{
rt<-scan("r/rt/data2/triam1.dat")
rt<-sort(rt)
plot(rt,ppoints(rt))
a<-9
b<-.27
fn<-function(p) -sum( log(dweibull(rt,p[1],p[2])) )
cat("starting -log like=",fn(c(a,b)),"\n")
out<-nlm(fn,p=c(a,b), hessian=TRUE)
2011 May 19
2
recursive function
Hi,
I created a function for obtaining the normal cumulative distribution (I
know all this already exists in R, I just wanted to verify my
understanding of it). below is the code I came up with.
cdf<-function(x) {
erf<-function(x) {
# approximation to the error function (erf) of the
# normal cumulative distribution function
# from Winitzki
2009 Jan 04
1
Bivarite Weibull Distribution
HI
Every one
Could some one provide me definitions of following bivariate distributions
gamma, exponencial, Weibull, half-normal , Rayleigh, Erlang,chi-square
thanks
A.S. Qureshi
2001 Jul 02
1
nls newbie: help approximating Weibull distribution
Hi folks,
I tried to retain the Weibull distribution using the `nls' function
and proceeding along the lines of the example provided in the
`SSweibull' help (at least I thought so):
t <- (1:200)/100
v <- pweibull(t, shape=3, scale=1)
df <- data.frame(Time=t, Value=v)
Asym <- 1.0; Drop <- 1.0; lrc <- 0; pwr <- 1
df.estimate <- nls(Value ~ SSweibull(Time,
2010 Jul 19
2
integral in R
Hello All,
I have to create a variable that is a function of another one (already
created), its cumulative distribution function and the integral of this
cumulative distribution, with limits: 0 and the value of the variable.
To be clear, I have the variable called “cip”. And its cdf called “cdfcip”
I need to create the variable:
bip = cip + ((1 – cdfcip)^4)*integral((1-cdf(u))^4*du, 0, value
2003 Sep 09
2
Computing a CDF or many quantiles
Given f, a pdf over a finite interval, is there any existing R function that
can efficiently tabulate the cumulative distribution function for f, or
produce all N+1 quantiles of the form i/N? "Efficiently" here means better
than doing repeated integrations for each point.
2008 Oct 07
3
Fitting weibull, exponential and lognormal distributions to left-truncated data.
Dear All,
I have two questions regarding distribution fitting.
I have several datasets, all left-truncated at x=1, that I am attempting
to fit distributions to (lognormal, weibull and exponential). I had
been using fitdistr in the MASS package as follows:
fitdistr<-(x,"weibull")
However, this does not take into consideration the truncation at x=1. I
read another posting in this
2005 Jul 07
1
CDF plot
Dear all,
I have define a discrete distribution P(y_i=x_i)=p_i, which I want to
plot a CDF plot. However, I can not find a function in R to draw it
for me after searching R and R-archive. I only find the one for the
sample CDF instead my theoretical one.
I find stepfun can do it for me, however, I want to plot some
different CDF with same support x in one plot. I can not manage how to
do it with