Displaying 20 results from an estimated 1000 matches similar to: "R's integrate function"
2008 Sep 27
3
Double integration - Gauss Quadrature
Hi,
I would like to solve a double integral of the form
\int_0^1 \int_0^1 x*y dx dy
using Gauss Quadrature.
I know that I can use R's integrate function to calculate it:
integrate(function(y) {
sapply(y, function(y) {
integrate(function(x) x*y, 0, 1)$value
})
}, 0, 1)
but I would like to use Gauss Quadrature to do it.
I have written the following code (using R's statmod package)
2008 Oct 19
1
multivariate integral with ADAPT when the parameter is close to boundary
Dear All,
There is one problem I encountered when I used ADAPT to compute some
2-D integral w.r.t beta density.
For example, when I try to run the following comments:
fun2<-function(theta){return(dbeta(theta[1],0.005,0.005)*dbeta(theta[2],0.005,0.005))}
int.fun2<-adapt(ndim=2,lo = c(0,0), up = c(1,1),functn = fun2,eps = 1e-4)
It seems it will take very long time to run. Acturally, I
2009 Oct 06
1
ldapsam, smbpasswd and posixAccount
Hi,
I have a setup with Samba (3.4) as PDC with ldapsam as backend. LDAP is
managed by Samba, no external helper scripts. When I add a new user with
"smbpasswd -a foo" it works fine, user is created and the openldap is
populated with samba and posix attributes.
Now, when I add a new user to the directory by running ldapadd against
a ldif file which contains only posix related
2011 Dec 01
3
Change the limits of a plot "a posteriori"
Hi all
How can I change the limits (xlim or ylim) in a plot that has been already
created?
For example, consider this naive example
curve(dbeta(x,2,4))
curve(dbeta(x,8,13),add=T,col=2)
When adding the second curve, it goes off the original limits computed by R
for the first graph, which are roughly, c(0,2.1)
I know two obvious solutions for this, which are:
1) passing a sufficiently large
2011 Aug 01
3
Beta fit returns NaNs
Hi,
sorry for repeating the question but this is kind of important to me and i
don't know whom should i ask.
So as noted before when I do a parameter fit to the beta distr i get:
fitdist(vectNorm,"beta");
Fitting of the distribution ' beta ' by maximum likelihood
Parameters:
estimate Std. Error
shape1 2.148779 0.1458042
shape2 810.067515 61.8608126
Warning
2001 Jun 06
3
error in dbeta (PR#970)
Full_Name: Hans Peter Wolf
Version: 1.2.1
OS: hpux10.20
Submission from: (NULL) (129.70.84.25)
dbeta computes a wrong result with parameters (1.3,1)
> version
platform hppa2.0-hp-hpux10.20
arch hppa2.0
os hpux10.20
system hppa2.0, hpux10.20
status
major 1
minor 2.1
year
2013 Sep 18
1
dbeta may hang R session for very large values of the shape parameters
Dear all,
we received a bug report for betareg, that in some cases the optim call in betareg.fit would hang the R session and the command cannot be interrupted by Ctrl-C?
We narrowed down the problem to the dbeta function which is used for the log likelihood evaluation in betareg.fit.
Particularly, the following command hangs the R session to a 100% CPU usage in all systems we tried it (OS X
2011 Jul 29
1
How to interpret Kolmogorov-Smirnov stats
Hi,
Interpretation problem ! so what i did is by using the:
>fit1 <- fitdist(vectNorm,"beta")
Warning messages:
1: In dbeta(x, shape1, shape2, log) : NaNs produced
2: In dbeta(x, shape1, shape2, log) : NaNs produced
3: In dbeta(x, shape1, shape2, log) : NaNs produced
4: In dbeta(x, shape1, shape2, log) : NaNs produced
5: In dbeta(x, shape1, shape2, log) : NaNs produced
6: In
2007 Nov 08
1
64-bit R-build on Mac OS X 10.4 - make check failures
Hi all,
I compiled 64-bit R on an Apple Mac G5 running OS X, but it failed
make check. Simon Urbanek suggested I post results to R-devel.
> On Nov 6, 2007, at 10:23 PM, Steven McKinney wrote:
>
> > Hi Simon,
> >
> > Would you be able to give more guidance on how to compile 64-bit
> > libiconv for Tiger,
>
> You can get the sources from Apple and compile
2007 Nov 24
2
how to compute highest density interval?
Suppose i want to compute a 95% highest density for a beta distribution
beta(a,b)
the two end points x1 and x2 shoudl satisfy the following two equations:
pbeta(x1,a,b)-pbeta(x2,a,b)=95%
dbeta(x1,a,b)=dbeta(x2,a,b)
Is there any fast way to compute x1 and x2 in R?
[[alternative HTML version deleted]]
2004 Nov 09
3
Strange results for Beta Distribution
Dear All,
I got these results from the example in the function "dbeta":
>x <- seq(0, 1, length=21)
> dbeta(x, 1, 1)
[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
Any Idea?
TIA
Giovanni
dr. Giovanni Parrinello
Section of Medical Statistics
Department of Biosciences
University of Brescia
25127 Viale Europa, 11
Brescia Italy
Tel: +390303717528
Fax: +390303701157
1997 Apr 15
1
R-alpha: Bug & Patch in dbeta.c (0.50 - PreR 7)
dbeta(1, a,b) would return 1 instead of 0.
Here is the patch for ..../src/math/dbeta.c :
--- dbeta.c~ Sun Nov 24 23:43:10 1996
+++ dbeta.c Tue Apr 15 21:25:30 1997
@@ -23,9 +23,7 @@
{
if (a <= 0.0 || b <= 0.0)
DOMAIN_ERROR;
- if (x <= 0)
+ if (x <= 0 || x >= 1.0)
return 0.0;
- if (x >= 1.0)
- return 1.0;
return MATH_CHECK(pow(x, a - 1) * pow(1.0 - x, b - 1.0) /
2012 Jan 03
6
calculate quantiles of a custom function
Hi,
I guess that my problem has an obvious answer, but I have not been able to
find it.
Suppose I create a custom function, consisting of two beta-distributions:
myfunction <- function(x) {
dbeta(x,2,6) + dbeta(x,6,2)
}
How can I calculate the quantiles of myfunction?
I have not seen any continous function treated in the docs, and applying the
"quantile function" gives me an
2009 Sep 04
2
plot positive predictive values
Hi,
I'm trying to fit a smooth line in a plot(y ~ x) graph.
x is continuous variable
y is a proportion of success in sub-samples, 0 <= y <= 1, from a Monte
Carlo simulation.
For each x there may be several y-values from different runs. Each run
produces several sub-samples, where "0" mean no success in any sub-
sample, "0.5" means success in half of the
2000 Aug 24
1
too large alpha or beta in dbeta ?
Dear friends.
Is this as expected ? Is alpha and beta too large simply ?
> dbeta(.1,534,646)
[1] NaN
Warning message:
NaNs produced in: dbeta(x, shape1, shape2, log)
Best wishes
Troels
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or
2000 Aug 25
1
Re: [R] too large alpha or beta in dbeta ? (PR#643)
>>>>> "TL" == Thomas Lumley <thomas@biostat.washington.edu> writes:
TL> On Thu, 24 Aug 2000, Troels Ring wrote:
>> Dear friends.
>>
>> Is this as expected ? Is alpha and beta too large simply ?
>>
>> > dbeta(.1,534,646)
>> [1] NaN
>> Warning message:
>> NaNs produced in:
2005 Jul 27
1
error message running R2WinBUGS
*Dear R-user,
*
I try to run Winbugs from R using bugs function in R2WinBUGS.My model works
well in Winbugs except that I can't get DIC. Since I don't need DIC, when I
try to run Winbugs from R , I set "DIC=FALSE". My model is as following:
model {
for (i in 1:N) {
for(j in 1 : T ) {
x[i, j] ~ dbin(p[i, j],n[i])
#Hier.prior
p[i, j] ~ dbeta(alpha[i, j], beta[i, j])
alpha[i, j]
2010 Jan 04
3
how to plot multiple density functions in one graph
Hello,
I am new to R and have two easy questions.
How can you plot multiple density functions in one graph? I have five beta
densities that I would like to plot in one graph. I understand how to plot
one beta density as a line:
plot (x,(dbeta(x,shape1=,shape2=,), type ="l")
Does the Pareto distribution need to be added to R with an additional
package?
thanks,
John
[[alternative
2013 Jan 22
2
Assistant
Good-day Sir,
I am R.Language users but am try to? estimate parameter of beta distribution particular dataset but give this error, which is not clear to me: (Initial value in "vmmin" is not finite)
beta.fit <- fitdistr(data,densfun=dbeta,shape1=value , shape2=value)
kindly assist.
expecting your reply:
2013 Jul 16
2
Problem following an R bug fix to integrate()
I have been told by the CRAN administrators that the following code generated
an error on 64-bit Fedora Linux (gcc, clang) and on Solaris machines (sparc,
x86), but runs well on all other systems):
> fn <- function(x, y) ifelse(x^2 + y^2 <= 1, 1 - x^2 - y^2, 0)
> tol <- 1.5e-8
> fy <- function(x) integrate(function(y) fn(x, y), 0, 1,