similar to: kernel density estimates

For the hist estimate >par(mex=1.3) >dens<-density(q) >options(scipen=4) > ylim<-range(dens$y) > h<-hist(q,breaks="scott",freq=FALSE,probability=TRUE, +? right=FALSE,xlim=c(9000,16000),ylim=ylim,main="Histogram of q(scott)") > lines(dens) >box() ? For the kernel estimate>options(scipen=4) > d <- density(q, bw =

Dear All, Attached are the codes of a histogram & a kernel density estimate and the output they produced. In fact the variable q is a vector of 1000 simulated values; that is I generated 1000 samples from the pareto distribution, from each sample I calculated the value of q ( a certain fn in the sample observations), and thus I was left with 1000 values of q and I don't know the

I am using the following commands for plotting kernel density for three kinds of crops density(s22$Net_income_Total.1, bw="nrd0",adjust=1, kernel=c("gaussian"))->t plot(t, xlim=c(-30000,40000), main="Net Income Distribution", axes=F, ylim=c(0,0.00035). xlab="Value in Rupees") par(new=T) density(s33$Net_income_Total.1, bw="nrd0",adjust=1,

I first answered to the email below in r-help, but as I did not see any response, and it looks like a bug/unwanted behavior, I am also posting here. I have observed this in RGui, whereas it seems not to happen in RStudio. Similar to OP, I sometimes have a problem with functions using the progress bar. Frequently, the console is cleared after x iterations when the progress bar is called in a

Hi all, I want the y-axis label to be ( in symbols) g(sigma given alpha) where given is the conditional sign. I've tried ylab=expression(g(sigma|alpha))) but it gave me g(|(sigma,alpha)) where the sigma and alpha are in greek but the conditional sign is misplaced (before the bracket) Any help would be appreciated. Maram. [[alternative HTML version deleted]]

Hi all, I've a question concerning the R 3.3.1 version. I have a long code that I used to run on versions earlier to the 3.3.1 version, and when I copied the code to the R console, I can still see the code while the loop is executing , along with the output printed after each iteration of the loop. Now, on the 3.3.1 version, after I copy the code to the console, it disappears and I only see

Hi all, I've a question concerning the R 3.3.1 version. I have a long code that I used to run on versions earlier to the 3.3.1 version, and when I copied the code to the R console, I can still see the code while the loop is executing , along with the output printed after each iteration of the loop. Now, on the 3.3.1 version, after I copy the code to the console, it disappears and I only see

Hi all, I've a linear equation of the form: 0.95=2 ( [3+ln(x/3)]^-13  + 4 [3+2ln(x/3)]^-13  + [3+4ln(x/3)]^-13 ) and I want to solve it for x, can I do this using R? Thanks in advance. Maram. [[alternative HTML version deleted]]

I have no idea to be honest. I have never used the package, I simply did a search for it. Hopefully a more experienced user can help --- On Sun, 7/12/09, maram salem <marammagdysalem at yahoo.com> wrote: > From: maram salem <marammagdysalem at yahoo.com> > Subject: Re: [R] (no subject) > To: "John Kane" <jrkrideau at yahoo.ca> > Received: Sunday, July 12,

Hi all, I have a question regarding the density function which gives the kernel density estimator. I want to decide the bandwidth when using gaussian kernel, given a set of observations. I am not familiar with different methods for bandwidth determination. Below are the different ways in R on deciding the bandwidth. Can anyone give an idea on which ones are preferred. Also, how can I take

Dear all, How can I use the histogram density estimate (hist) to find the value of the cdf at a certain point? Thanks Maram [[alternative HTML version deleted]]

I would like to ask once more if this is reproducible also for others? If yes, should I submit it as a bug-report? Best, Jon On 11/28/2016 11:26 AM, Jon Skoien wrote: > I first answered to the email below in r-help, but as I did not see > any response, and it looks like a bug/unwanted behavior, I am also > posting here. I have observed this in RGui, whereas it seems not to >

Hi all, I'm using a certain  procedure to calculate the value of some variable(Bayes risk),B. So I got the values B1, B2, ........, B1000, each under certain input values and using a long procedure. Now, I want to put the values I got in a nummerical vector and find their minimum value. I think c( ) should work.For example if I have only 10 values I could have used

Dear all, I have the cdf of the following power fuction distribution: F(y)=(y/350)^a               ,0<y<350, where " a " is some parameter with range a>0. I want to use it as the argument of the discretize function of the actuar package. So I think I need to define this function to R so that if I entered a=1, I get the following F(y)=(y/350) and if I entered a=4.5, I get the

Hi all, I've a trivial question. If (q) is a continous variable,actually a vector of 1000 values. how to calculate the probability that q is greater than a specific value, i.e. P(q>45)?? Thanks Maram [[alternative HTML version deleted]]

Dear All, (this msg is a statistics/computing question to the list) I'm trying to implement a modern-version of a (classic) "several-step protocol" in Fishery Biology (due to Bhattacharya, 1967): analysis of length-frequency distribution of fish larvae to id cohorts and later estimate growh rates! I've trouble with the 1st step: using kernel density estimation (KDE) function

--- On Tue, 6/30/09, maram salem <marammagdysalem at yahoo.com> wrote: > From: maram salem <marammagdysalem at yahoo.com> > Subject: [R] (no subject) > To: r-help at r-project.org > Received: Tuesday, June 30, 2009, 6:34 AM > Hi Group, > I've a vector of 1000 numeric values for which I want to > draw a histogram. I've read this vector into R with no >

Dear All, I have a matrix m of the form  m             [,1]             [,2]   [1,]  9072.302 0.00000004462366   [2,]  9086.811 0.00000005628169   [3,]  9101.320 0.00000007126347   [4,]  9115.830 0.00000008986976   [5,]  9130.339 0.00000011260000   [6,]  9144.848 0.00000014018405   [7,]  9159.357 0.00000017344229   [8,]  9173.866 0.00000021363563   [9,]  9188.375 0.00000026389996  [10,]  9202.884

Sorry, I didn't mean a linear equation, of course this equation is not linear, I meant an equation in one unknown. Hi all, I've a linear equation of the form: 0.95=2 ( [3+ln(x/3)]^-13? + 4 [3+2ln(x/3)]^-13? + [3+4ln(x/3)]^-13 ) and I want to solve it for x, can I do this using R? Thanks in advance. Maram. ? ? ? ??? [[alternative HTML version deleted]] -------------- next

The following is from version 1.6.2 of R under Windows, or 1.6.1 under Mac OSX/X11 > density(1) Error in if (!(lo <- min(hi, IQR(x)/1.34))) (lo <- hi) || (lo <- abs(x[1])) || : missing value where logical needed I am not sure how this should be handled. I encountered it in connection with densityplot(). In that connection, it might be enough to modify density() so that it