search for: cdf

Displaying 20 results from an estimated 360 matches for "cdf".

Did you mean: ecdf
2007 Jul 10
3
ECDF, distribution of Pareto, distribution of Normal
Hello all, I would like to plot the emperical CDF, normal CDF and pareto CDF in the same graph and I amusing the following codes. "z" is a vector and I just need the part when z between 1.6 and 3. plot(ecdf(z), do.points=FALSE, verticals=TRUE, xlim=c(1.6,3),ylim=c(1-sum(z>1.6)/length(z), 1)) x <- seq(1.6, 3, 0.1) lines(x,pgpd(x,...
2011 Feb 17
0
[BioC] Make.cdf.package error
Hi everybody, I tried to analyze a custom Affymetrix 3'-biased Array. So I wanted to make a cdf package. (My CDF file size is 1.12Go). I tried several methods but the same error occured Method 1 > #Set the working directory > setwd("D:/Analyse R/Cel files") > #library to create cdf env > library("makecdfenv") >#Create cdf environment >pkgpath &l...
2009 Dec 28
2
[BioC] make.cdf.package: Error: cannot allocate vector of size 1 Kb
My machine has 8GB memory. I had quit all other programs that might take a lot of memory when I try the script (before I post the first message in this thread). The cdf file is of only 741 MB. It is strange to me to see the error. On Mon, Dec 28, 2009 at 2:38 AM, Wolfgang Huber <whuber at embl.de> wrote: > Dear Peng Yu > > how big is the RAM of your computer? You could try with closing all other > applications before running this script. You cou...
2006 Dec 04
0
How to calculate area between ECDF and CDF?
...orov-Smirnov test, but of course this is invalid when the distribution parameters are estimated from the data. So I'm tinkering with an alternative method (excuse my naivet? if this is a bad idea, I'm a relative statistical novice) that calculates the area of the difference between the ECDF of the data and the CDF of the estimated function (somewhat like KS, which looks at the greatest distance between these). My thought is to compare this observed area to a distribution of simulated areas derived by monte carlo simulation (draw N random samples from the estimated function, calcul...
2010 Feb 09
1
how to adjust the output
Hi R-users,   I have this code below and I understand the error message but do not know how to correct it.  My question is how do I get rid of “with absolute error < 7.5e-06” attach to value of cdf so that I can carry out the calculation.   integrand <- function(z) { alp  <- 2.0165   rho  <- 0.868   # simplified expressions   a      <- alp-0.5   c1     <- sqrt(pi)/(gamma(alp)*(1-rho)^alp)   c2     <- sqrt(rho)/(1-rho)   t1     <- exp(-z/(1-rho))   t2     <- (z/(2*c2))...
2007 Sep 12
1
Integrate() error message, I am at a loss
.... DO NOT PASS it! alpha <- mu/sqrt((1+theta/n)) normconst <- 1 if(renorm==1 & rec==0){ term1 <- integrate(nctspa, -Inf, alpha, n=n, mu=mu, theta=theta)$value term2 <- integrate(nctspa, alpha, Inf, n=n, mu=mu, theta=theta)$value normconst <- 1/(term1+term2) } cdf <- numeric() pdf <- cdf c3 <- n^2+2*n*a^2+a^4 c2 <- (-2*mu*(a^3+n*a))/c3 c1 <- (-n^2-n*a^2-n*theta+a^2*mu^2)/c3 c0 <- (n*a*mu)/c3 q <- c1/3-(c2^2)/9 r <- 1/6*(c1*c2-3*c0)-1/27*c2^3 b0 <- sqrt(-4*q)*cos(acos(r/sqrt(-q^3))/3)-c2/3 t1 <- -mu+a*b0 t2 <- -a*t1/b0/n/2 nu...
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 d...
2012 Jun 14
2
plot cdf
Good Afternoon, I'm trying to create a cdf plot, with the following code. It works well, but I have little doubt, if you can help solve. When I create the plot, like the graph line would still not appear with point #cdf x<-table(Dataset$Apcode) View(s) hist(s) *plot(ecdf(x))* x<-1 37607 2 26625 3 5856 4 25992 5 30585 6 16064 7...
2012 Jul 11
2
Computing inverse cdf (quantile function) from a KDE
Hello, I wanted to know if there is a simple way of getting the inverse cdf for a KDE estimate of a density (using the ks or KernSmooth packages) in R ? The method I'm using now is to perform a numerical integration of the pdf to get the cdf and then doing a search for the desired probablity value, which is highly inefficient and very slow. Thanks, -fj [[alternativ...
2006 Sep 23
1
Fitdistr() versus nls()
Hello R-Users, I'm new to R so I apologize in advance for any big mistake I might be doing. I'm trying to fit a set of samples with some probabilistic curve, and I have an important question to ask; in particular I have some data, from which I calculate manually the CDF, and then I import them into R and try to fit: I have the x values (my original samples) and the y values (P(X<x)). To attempt the fit I've both fitdistr() and nls(), in the way you can read in the piece of code at the end of the email. Because the fit with all data doesn't work...
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.
2010 Sep 17
1
How to denote a line on the graph
Please give me some help, many thanks. I graphed a expected CDF line of a binomial distribution on a graph, And I have some observed points (observed CDF) from 4 groups fall on the smooth CDF line. I cannot really get the legend I want legend ('topleft', c('a, 'b', 'c', 'd', 'expected CDF'), col=c(1,2,3,4), pch=c(0...
2010 Nov 22
1
need smooth cdf lines
Hi, I would like to overlap the cdf curve for observed and generated data  Here is my code: plot(cdf,main ="CDF of the sum for winter season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", xaxs="i",yaxs="i",col=c("black","red"), lty=c(1,1),ylab="Cumulative probability&quot...
2003 May 08
2
approximation of CDF
Hi all, is there any package in R capable of smooth approximation of CDF basing on given sample? (Thus, I am not speaking about ecdf) In particular, I expect very much that the approximation should subject to the property: f(x0)<=f(x1) for x0<x1, where x0 and x1 belong to range of the sample given. Polynomial approximation could be OK f...
2005 Jul 06
0
download.file() yields incomplete files with method="internal" (PR#7991)
Summary: When I use method="wget" with download.file(), I consistently get a download of the entire file. When I use method="internal", I infrequently get the entire file, but usually get only part of the file. This behavior occurs with .cdf (a weather file format - basically binary) from a UCAR site. I am not sure this is a bug, since it could be some internet inconsistency, or I may be doing something wrong, but given that wget works and internal doesn't, it might be a bug. Appended to this message is an example of the output f...
2007 Sep 12
0
Problem with integrate()
...ent called rec. DO NOT PASS it! alpha <- mu/sqrt((1+theta/n)) normconst <- 1 if(renorm==1 & rec==0){ term1 <- integrate(nctspa, -Inf, alpha, n=n, mu=mu, theta=theta)$value term2 <- integrate(nctspa, alpha, Inf, n=n, mu=mu, theta=theta)$value normconst <- 1/(term1+term2) } cdf <- numeric() pdf <- cdf c3 <- n^2+2*n*a^2+a^4 c2 <- (-2*mu*(a^3+n*a))/c3 c1 <- (-n^2-n*a^2-n*theta+a^2*mu^2)/c3 c0 <- (n*a*mu)/c3 q <- c1/3-(c2^2)/9 r <- 1/6*(c1*c2-3*c0)-1/27*c2^3 b0 <- sqrt(-4*q)*cos(acos(r/sqrt(-q^3))/3)-c2/3 t1 <- -mu+a*b0 t2 <- -a*t1/b0/n/2 nu...
2008 Dec 11
2
how to get the CDF of a density() estimation?
Hi, I've estimated a simple kernel density of a univariate variable with density(), but after I would like to find out the CDF at specific values. How can I do it? thanks for your help, with it I am very close to finish my first little bit more serious work in R, Viktor
2008 Dec 16
1
How to make a smooth ( linear ) CDF plot?
This question might seem silly, because I felt that it MUST be in the mailing list archives or help files somewhere, but I simply couldn't find it. I want to make some simple CDF (cumulative distribution function) plots to check whether distributions are Gaussian / normal. But in order to check how "normal" the distribution is, I really need the y-axis to be Gaussian as well (meaning that the axis is not linearly scaled), so that if the distribution is indeed nor...
2017 Oct 07
1
beta binomial distribution
Hi, I need to write two inequalities depend on cumulative distribution (CDF) of beta binomial distribution where alpha and beta are unknown and need to find them. CDF of betabinomial(2,10,alpha,beta) <0.3<=CDF of betabinomial(3,10,alpha,beta) and CDF of betabinomial(5,10,alpha,beta) <0.8<=CDF of betabinomial(6,10,alpha,beta) How I can do that using r studio p...
2010 Nov 09
1
agrep pmatch recursive???
...ess - 64 bit windows 7, R 2.11.1 I am trying to match the character contents of one list, called 'exclude', to those of a second list, called 'dataset' dataset is a list of file names with folder locations, and looks like this when called: > dataset [1] "A/10-10-29a-13.cdf" "A/10-10-29a-14.cdf" "A/10-10-29a-15.cdf" "A/10-10-29a-16.cdf" "A/10-10-29a-17.cdf" .....(about 20 file location/names in this case)..... I also have the 'exclude' list, which looks like this when called: > exclude Row.names 1 10-10-2...