Displaying 20 results from an estimated 3000 matches similar to: "question on integrate function & request for consultant"
2009 Sep 29
1
Probability of data values form empirical distribution
Hello,
Could someone help me please and to tell how to get the probability from empirical distribution (not parametric) for each data value (R function).
For example, for normal distribution there is such a function like:
“pnorm(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)”
I need the same function only for the empirical distribution (which does not correspond to any typical
2009 Sep 29
1
Probability of data values form empirical distribution
Hello,
Could someone help me please and to tell how to get the probability from empirical distribution (not parametric) for each data value (R function).
For example, for normal distribution there is such a function like:
“pnorm(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)”
I need the same function only for the empirical distribution (which does not correspond to any typical
2010 May 06
1
question about rolling regressions
Hi All,
I am using R 2.11.0 on a Ubuntu machine. I have a time series data set and
want to run rolling regressions with it. Any suggestions would be useful.
Here are the details:
(1) I convert relevant variables into time series objects and compute first
differences:
vad <- ts(data$ALLGVA/data$GDPDEF, start=1948, frequency=1)
emp <- ts(data$ALLEMP, start=1948, frequency=1)
vad.dif1 <-
2009 Nov 27
1
problem with "dynformula" from "plm" package [RE-POST]
Hello list,
I'm following the paper (http://www.jstatsoft.org/v27/i02/paper) on
how to use "plm" to run panel regressions, and am having trouble with
what I believe should be something very basic.
When I run the command (p.9 in the paper):
R>
dynformula(emp~wage+capital,log=list(capital=FALSE,TRUE),lag=list(emp=2,c(2,3)),diff=list(FALSE,capital=TRUE))
I see:
emp ~ wage +
2010 Jan 25
3
Issue using tapply
Hello all,
I am trying to use the tapply function to sum some values and change the
column names of the resulting vector.
I input
Emp Et
1 10565 ACC
2 7515 ADM
3 625 AGF
4 6243 CNS
5 12721 EDU
6 3924 FIN
7 18140 HLH
8 3686 INF
9 15841 MFG
10 243 MIN
11 1864 MNG
12 4664 OSV
13 5496 PRF
14 4988 PUB
15 2166 REC
16 2153 REL
17 16082 RTL
18 3582 TRN
19 757 UTL
20
2005 Mar 08
1
To convert an adjacency list model into a nested set model
Dear R-help
I am wondering if somebody wrote some code to convert an adjacency list
model into a nested set model.
In principal I want to do the same as John Celko mentioned it here with
SQL:
http://groups.google.co.uk/groups?hl=en&lr=lang_en&selm=8j0n05%24n31%241
%40nnrp1.deja.com
Assume you have a tree structure like this
Albert
/ \
/
2011 Feb 02
2
Efficient way to determine if a data frame has missing observations
I have a data set covering a large number of cities with values for characteristics such as land area, population, and employment. The problem I have is that some cities lack observations for some of the characteristics and I'd like a quick way to determine which cities have missing data. For example:
2011 Oct 16
2
ecdf
Hi,
Newbie here. I read the R for Beginners but i still don't get this.
I have the following data (this is just an example) in a CSV file:
courseid numstudents
101 209
141 13
246 140
263 8
321 10
361 10
364 28
365 25
366 23
367
2011 Jan 26
2
2 functions with same name - what to do to get the one I want
There seems to be 2 functions call ecdf...
http://lib.stat.cmu.edu/S/Harrell/help/Hmisc/html/ecdf.html
http://127.0.0.1:11885/library/stats/html/ecdf.html
How do I get the one ecdf {Hmisc} to run instead of the ecdf {stats}
A pointer in the right direction would be greatly appreciated.
Tried to instal Hmisc but got this message, so I assume I have it
> utils:::menuInstallPkgs()
Warning:
2011 Feb 10
1
How to determine the quantile boundary from an ECDF?
Given a dataset x, the ecdf is ecdf(x). Then I can use ecdf(x)(y) to find
the percentile of y. Given the ecdf is there a way to determine what is the
value of y that is the boundary of let's say 95 percentile? In other words,
is there a function I can call on the ecdf like:
fomeFunc( ecdf( x ), 0.95 )
Which will return the highest value of y, for which ecdf( y ) < 0.95?
The only solution
2003 Jun 18
2
Forward stepwise procedure w/ stepAIC
I'm attempting to select a model using stepAIC. I want to use a forward
selection procedure. I have specified a "scope" option, but must not be
understanding how this works. My results indicate that the procedure begins
and ends with the "full" model (i.e., all 17 independent variables)...not
what I expected. Could someone please point out what I'm not
2012 May 02
1
Two ecdf with log-scales
Hi,
i want to plot empirical cumulative density functions for two variables in
one plot. For better visualizing the differences in the two cumulative curves I'd like to log-scale the axis.
So far I found 3 possible functions to plot ecdf:
1) ecdf() from the package 'stats'. I don't know how to successfully set the log.scales? Combining two plots is not a problem:
plot(ecdf(x1))
2006 Apr 20
2
smooth the ecdf plots
Hi All,
I have codes as follows to get the ecdf plots:
>
day.hos2<-c(6,4,6,6,4,6,5,4,7,5,6,6,8,6,17,9,8,4,6,3,5,8,7,12,5,10,6,4,6
,13,7,6,6,25,4,9,96,6,6,6,6,9,4,5,5,4,10,5,7,6)
>
day.hos3<-c(5,6,7,6,4,5,6,6,6,6,19,7,5,9,8,8,7,5,6,20,40,5,8,7,7,5,6,13,
11,9,4,6,9,16,6,7,6)
> f<-ecdf(day.hos2)
> plot(f,col.p='red',col.h='red')
>
2010 Apr 19
2
ecdf
Hello,
I'd like to plot an empirical cumulative distribution function, except
instead of the fraction of values < x, I'd like the fraction of values > x.
I think this can be done using the ecdf function in {Hmisc}. I installed
the package and loaded it. However, when following the example given in the
documentation, I get an error:
x <- rnorm(100)
ecdf(x,what='1-F')
2012 May 26
1
Kolmogorov-Smirnov test and the plot of max distance between two ecdf curves
Hi all,
given this example
#start
a<-c(0,70,50,100,70,650,1300,6900,1780,4930,1120,700,190,940,
760,100,300,36270,5610,249680,1760,4040,164890,17230,75140,1870,22380,5890,2430)
length(a)
b<-c(0,0,10,30,50,440,1000,140,70,90,60,60,20,90,180,30,90,
3220,490,20790,290,740,5350,940,3910,0,640,850,260)
length(b)
out<-ks.test(log10(a+1),log10(b+1))
# max distance D
2011 Jul 16
2
ecdf() to nls() - how to transform data?
Hi,
I am using ecdf-function and want to use the ecdf()-data-points for nls() as
data-parameter.
nls() expects 'list' or 'environment' as a data-type, knots(ecdf(mydata))
gives me 'numeric'.
What should I do now?
Thanks in advance - Jochen
Here is the code:
#################################################
# --- Fit ---
# Gumbel-Dist-Function, cumulative,
2010 Nov 25
1
overlap cdf plots and add colors and etc
Hi r-users,
I would like to overlap 2 ecdf plots.
I tried this below and it gives me two plots of ecdf but just both just in
black.
par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i")
plot(ecdf(datobs))
lines(ecdf(gam_sum_gen))
Then I try to add colors etc and also the legend but fail.
par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i")
2004 Oct 26
2
Empirical P Value
I am trying to return the p value for a stat from the ECDF. That is the
index of the first occurrence,
on an ordered vector, of a value either greater than or equal to a given
value.
Ideally I would not have to order the vector beforehand.
Currently I use:
PValue<-function(stat, ECDF){
###Get the length of the ECFD
L<-length(ECDF)
###Loop through the ECDF until the p value is
2016 Jun 24
3
Ayuda ggplot2
Hola a todos!
Soy nueva en R y necesito hacer unos gráficos para una investigación, he
explorado un poco y estoy intentando usar ggplot2 ya que hace gráficos de
muy buena calidad...tengo los datos de varios años para diferentes grupos
de empresas y los pretendo graficar tanto en un solo grafico como en varios
(facet_wrap) pero tengo problemas con el eje de las X, ya que necesito que
aparezcan los
2012 Apr 24
2
Function from ecdf to ccdf
Dear all,
I would like to calculate the complementary cumulative distribution function. As it is known, the ccdf is the 1-ecdf(X)==1-F(x). (From ?ecdf help is shown that ecdf returns a function which is the F(x)
I would like to use that function inside my function and after call curve (which accepts function as input). My code looks like that
onVector<-seq(1:100)