Displaying 20 results from an estimated 900 matches similar to: "defining and plotting functions thanks to equation"
2003 Jul 15
3
clearing some part of a graph
Hi R lovers
2 questions:
1) I'd like to know how to clean the title, the sub title or the labels of
a graph. I know how to redefine it with the function title() but it
overwrites the previous title and do not replace it
2) How could I clear a whole plot (for example in a multiple figure
environment)
thanks for your help
vincent
2003 Jul 25
5
named list 'start' in fitdistr
Hi R lovers!
I'd like to know how to use the parameter 'start' in the function
fitdistr()
obviously I have to provide the initial value of the parameter to optimize
except in the case of a certain set of given distribution
Indeed according to the help file for fitdistr
" For the following named distributions, reasonable starting values
will be computed if `start'
2003 Jul 15
1
function acf in package ts
Hi R lovers!
I'd like to know if the 2 blue lines in the graph produced by the function
acf in the package ts represents the level for the test of significance of
the autocorrelation
thanks for help
vincent
******************************************************************
The sender's email address has changed to
firstname.lastname@ sgcib.com. You may want to update your
2003 Aug 05
1
error message in fitdistr
Hi R lovers
Here is a numerical vector test
> test
[1] 206 53 124 112 92 77 118 75 48 176 90 74 107 126 99 84 114
147 99 114 99 84 99 99 99 99 99 104 1 159 100 53
[33] 132 82 85 106 136 99 110 82 99 99 89 107 99 68 130 99 99
110 99 95 153 93 136 51 103 95 99 72 99 50 110 37
[65] 102 104 92 90 94 99 76 81 109 91 98 96 104 104 93 99 125
89
2005 Dec 03
1
Fit Frechet Distribution
hello everybody
i want to use the maximum likelihood method to estimate FRECHET parameters of
my sample data.
Should it work with fitdistr in the package MASS?
I only find how to do it for GEV, Gumbel, and almost all other distributions,
but FRECHET?
I would be very happy if somebody can tell me how to do fit the FRECHET
distribution!
Thanks
Nadja Riedwyl
2008 Sep 08
1
Gumbell distribution - minimum case
Hello,
I would like to sample from a Gumbell (minimum) distribution. I have
installed package {evd} but the Gumbell functions there appear to refer
to the maximum case. Unfortunately, setting the scale parameter
negative does not appear to work.
Is there a separate package for the Gumbell minimum?
--
_____________________________
Rich Gwozdz
Fire and Mountain Ecology Lab
College of
2000 Nov 09
2
simple mixture
Dear All,
I am trying to do some simple mixture analyses. For instance, I have a
sample of n observations and I suspect they come from two different
exponential distributions with parameters rate1 and rate2, respectively.
So, I want to estimate rate1, rate2, and the proportions of both kinds of
individuals in the sample. I had a look at the packages mda and mclust, but
they do not seem to do this
2003 Nov 06
3
import data troubles
HI R lovers,
I have been facing a petty trouble with datas' import :
I have a plain txt file (see attached file or the copy below) that I cannot
read either with scan or read.table
> scan(file="F:/Alt/HDG/Stoliaroff/Data/test.txt")
Error in scan(file = "F:/Alt/HDG/Stoliaroff/Data/test.txt") :
"scan" expected a real, got "??6"
>
2003 Nov 24
2
merging variable and character
Hi R lovers! I'd like to know the trick to get a plot with a title
associating variable and text
I'd like to plot something like
for (i in 1:10) plot(x[,i],ylab="Serie" & i)
to get Serie 1 Serie 2 so on and so forth... of course it doesn't work.
If somebody could give me an idea, it would be very nice
thank you very much and long life to R
Vincent
2009 Aug 24
2
[] for R
I am assuming the variable out is the output parameter.
However, I don't understand what is out[1+xi*xx>0]?
Can someone explain this to me?
Thanks in advance,
Chen
dGEV <- function(x, xi, mu = 0, sigma = 1, logvalue=FALSE)
{
xx <- (x-mu)/sigma
#use the new dGumbel which passes mu and sigma:
#if (xi==0) out <- dGumbel(xx,logvalue=TRUE)-log(sigma)
if (xi==0) {
2003 Jun 06
3
writing comments within a function
Hi R lovers!
I would like to know how you can write comments inside the code of a
function
is it latex style
% ?
or any other language style?
thanks a lot
*************************************************************************
Ce message et toutes les pieces jointes (ci-apres le "message") sont
confidentiels et etablis a l'intention exclusive de ses destinataires.
Toute
2005 Jul 27
3
fitting extreme value distribution
hi,
rgev function gives me random deviates and I have a data
set which I am fitting to an EVD,IS there a way I can plot
both observed and ideal evd on the same plot
thankyou
Rangesh
2010 Jan 07
1
Return values in fExtremes package
Hi,
I was usuing the fExtemes package, and wanted to obtain some of the values returned from the function gumbelFit(). For example, in the following code, I would like to access 'mu' and 'beta' from the object 'para'. How should I go about doing this? Is there any generic method to access the object?
-----------------------------------
>library("fExtremes")
2011 Dec 13
1
tcplot documentation in evd package
Hello, and please advise regarding any errors/omissions on my part.
However, the documentation for R's tcplot function in the evd package
appears to contain an error. I am using evd version 2.2-4 in R 2.14.0
with Windows XP.
> data(portpirie)
> mrlplot(portpirie) ## No Error
> tlim <- c(3.6, 4.2)
> tcplot(portpirie, tlim) ## Line from documentation
Error in fpot(data, u[1],
2003 Feb 09
3
Clustering partition and memory
Dear R-help list members
i would like to use R to produce clustering or partitioning of a dataset.
I am trying to use the functions:
- hierclust() of the package multiv
-pam(), agnes() and fanny() of the package cluster
But I cannot get any result because of lack of memory. Would you know any
clustering function not to greedy in memory?
I have tried to expand my memory limit with memory.limit()
2013 Feb 12
2
standard error very high in maximum liklihood fitting
Dear all,
I have been trying to fit my data (only right censored) with gumbel distribution using fitdistrplus. I am getting very high standard error. I have been wondering why.
The followings are the outputs:
fit1=fitdistcens(dr0, "gumbel", start=list(a=99, b=0.6), optim.method= "L-BFGS-B", lower = 0.0, upper = Inf)
> summary(fit1)
FITTING OF THE DISTRIBUTION ' gumbel
2011 Sep 18
2
calculating VAR of a (Gumbel) copula
Hello,
I am a new user of R (2.13.1), my operational system is Windows Vista.
I have a problem with the attached file SFEVaRHAC.r, calculating the VAR of
a Gumbel copula, based on the attached GumHAC_VaR_PL_w250_n1000_s2500.txt
1. I had a Error in file(file, "rt") : cannot open the connection message.
I solved it by reading a post in nabble to use setwd(choose.dir()) and
2009 Aug 17
2
Newbie that don't understand R code
I got some R code that I don't understand.
Question as comment in code
//where is t comming from, what is phi inverse
rAC <- function(name, n, d, theta){
#generic function for Archimedean copula simulation
illegalpar <- switch(name,
clayton = (theta < 0),
gumbel = (theta < 1),
frank = (theta < 0),
BB9 = ((theta[1] < 1) | (theta[2] < 0)),
GIG = ((theta[2] < 0) |
2009 Mar 16
1
Fw: Fitting GUMBEL Distribution - CDF function and P P Plot
Dera R Helpers,
I am re-posting my query.
Please guide me.
Maithili
--- On Fri, 3/13/09, Maithili Shiva <maithili_shiva at yahoo.com> wrote:
I am trying to fit the Gumbel distribution to a data. I am
using lmom package. I am getting problem in Cumulative
Distribution Function of Gumbel distribution as I am getting
it as a series of 0's and 1's thereby affecting the
P P
2008 Apr 22
4
how to convert non numeric data into numeric?
I am having the following error in my function
function(theta,reqdIRR)
{
theta1<-theta[1]
theta2<-theta[2]
n<-length(reqdIRR)
constant<- n*(theta1+theta2)
sum1<-lapply(reqdIRR*exp(theta1),FUN = sum)
sum2<-lapply(exp(theta2 - reqdIRR*exp(theta1)),FUN = sum)
sum = sum1 + sum2
log.fcn = constant - as.numeric(sum)
result = - log.fcn
return(result)
}
*error :