Displaying 20 results from an estimated 3000 matches similar to: "64-bit R-build on Mac OS X 10.4 - make check failures"
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
2008 Jul 27
1
64-bit R on Mac OS X 10.5.4
Hi Matt
Your method is the easiest way for me to install the 64-bit R. I followed the directions on your web site and then did the following:
R --arch=x86_64
source("http://bioconductor.org/biocLite.R")
biocLite(type = "source",lib = "/Library/Frameworks/R.framework/Versions/2.8/Resources/RLib64")
I got many errors and warnings which I copied to the attached file.
2008 Jul 26
1
64-bit R on Mac OS X 10.4.5
Hello
I have a Mac OS X 10.4.5. I am trying to build a 64-bit R by following the directions on this page: http://r.research.att.com/building.html
r_arch=x86_64 \
CC="gcc -arch x86_64 -std=gnu99" \
CXX="g++ -arch x86_64" \
OBJC="gcc -arch x86_64" \
F77="gfortran -arch x86_64" \
FC="gfortran -arch x86_64"
PATH=/usr/X11/bin:/usr/local/bin:$PATH
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
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
2006 Dec 19
3
Bug in rt() ? (PR#9422)
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
<<insert bug report here>>
Reproduced on Debian and Windows ...
On 2.4.x if you execute
set.seed(12345)
t.1 <- rt(n = 1000, df = 20)
set.seed(12345)
t.2 <- rt(n = 1000, df = 20, ncp = 0)
all.equal(t.1, t.2) ## Not close to true
This appears to be due to the fact that in 2.4.x rt is now
rt
function (n, df, ncp = 0)
{
if
2008 Jun 14
1
qt with ncp>37.62
help(qt) states that:
"ncp non-centrality parameter delta; currently except for rt(), only for
abs(ncp) <= 37.62"
so I would expect that calling qt with non-centrality parameter exceeding
37.62 should fail, instead e.g. calling
> mapply(function(x) qt(p = 0.9, df = 55, ncp = x),35:45)
gives:
[1] 40.21448 41.35293 42.49164 43.68862 44.82945 45.97048 47.11170 48.25310
[9]
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]]
2008 Sep 30
2
R's integrate function
Hello,
I am trying to use R's integrate function to calculate the following
integral for z=423:
integrate(function(y,z){
sapply(y, function(y,z){
integrate(function(x,z)
1/x*dbeta(0.01,x/(0.005/1.005),(1-x)/(0.005/1.005))*dbeta(y,x/(0.005/1.005),(1-x)/(0.005/1.005))*(1-y)^z,0,1,423)$value
})
},0,1,423)$value
but I receive an error message saying that the maximum number of
subdivisions is
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) /
2007 Jan 31
2
Bug in 'pchisq' for x=0.0 (PR#9485)
The function 'pchisq' from the 'stats' library gives a wrong result if the
argument equals exactly zero:
# Upper tail of central 1-df chi^2 distribution
> pchisq(1 , 1, ncp=0, lower.tail = F, log.p = FALSE)
[1] 0.3173105
> pchisq(0.5 , 1, ncp=0, lower.tail = F, log.p = FALSE)
[1] 0.4795001
> pchisq(0.01 , 1, ncp=0, lower.tail = F, log.p = FALSE)
[1]
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
2006 Dec 10
1
Noncentral t & F distributions
Dear List:
The square of the noncentral t-statistic with noncentrality parameter
\delta is a noncentral F with noncentrality parameter \lambda=\delta^2.
So, t^2_{\nu,\delta} = F_{1,\nu,\lambda=\delta^2}. Consequently, it
should follow that t^2_{1-\alpha/2,\nu,\delta} =
f_{1-alpha,1,\vu,\lambda=\delta^2}. However, this is not what is
happening with the following code. The central
2002 Oct 17
3
Non-central distributions
Hi Folks,
I note that, while the "chisq" functions
dchisq(x, df, ncp=0, log = FALSE)
pchisq(q, df, ncp=0, lower.tail = TRUE, log.p = FALSE)
qchisq(p, df, ncp=0, lower.tail = TRUE, log.p = FALSE)
rchisq(n, df, ncp=0)
all have a slot for the non-centrality parameter "ncp", of
the functions for the t and F distributions:
dt(x, df, log = FALSE)
2000 Nov 28
1
non-centrality parameter in pf() (PR#752)
Bug Description:
Problem with the function pf() when the non-centrality
parameter is large. Here is a sample command. You should
see a smooth line from 0 to about 55, and then the values
of pf() go crazy from 55 to 100.
############################
ncp <- seq(0,100,length=200)
plot(ncp,pf(5,7,2,ncp=ncp))
############################
Version:
platform = i686-pc-linux-gnu
arch = i686
os =