Displaying 20 results from an estimated 684 matches for "rho".
Did you mean:
rao
2017 Jun 27
0
Seg Fault memory violation
....7 and 1.8 respectively on the same system. We have multiple installations and they all exhibit the same behavior. Any help would be greatly appreciated. Thanks for your time all.
-------------- next part --------------
Program received signal SIGBUS, Bus error.
bcEval (body=body at entry=0x857990, rho=rho at entry=0x17f3d798, useCache=useCache at entry=TRUE) at eval.c:5873
5873 {
(gdb) bt
#0 bcEval (body=body at entry=0x857990, rho=rho at entry=0x17f3d798, useCache=useCache at entry=TRUE) at eval.c:5873
#1 0x00007ffff78d4138 in Rf_eval (e=0x857990, rho=rho at entry=0x17f3d798) at eval.c:624...
2006 May 09
1
Seg fault when installing package from bad repository
...6290 in do_download (call=0x18d0078, op=0x181c2a8,
args=0x18c3b50, env=0x18cfd64) at
../../../../R-2.3.0/src/main/internet.c:99
#8 0x010cefc4 in do_internal (call=0x18c3b50, op=0xbfff2a03,
args=0x0, env=0x18cfd64) at ../../../../R-2.3.0/src/main/names.c:1089
#9 0x0109d798 in Rf_eval (e=0x18d0040, rho=0x18cfd64) at
../../../../R-2.3.0/src/main/eval.c:391
#10 0x0109f87c in do_set (call=0x18d17b0, op=0x180a190,
args=0x18d0008, rho=0x18cfd64) at
../../../../R-2.3.0/src/main/eval.c:1303
#11 0x0109d798 in Rf_eval (e=0x18d17b0, rho=0x18cfd64) at
../../../../R-2.3.0/src/main/eval.c:391
#12 0x0109d798 i...
2003 Apr 16
0
vignette segfault bt
.../ext/R/cvs/R-devel/src/main/gram.y:1083
#9 0x081175d1 in do_parse (call=0x883d674, op=0x8205024,
args=0x8839138,
env=0x9207e84) at source.c:76
#10 0x080c99e4 in do_internal (call=0x883d5e8, op=0x81fb5e8,
args=0x92088bc,
env=0x9207e84) at names.c:1056
#11 0x080a3b28 in Rf_eval (e=0x883d5e8, rho=0x9207e84) at eval.c:418
#12 0x080a4f62 in do_begin (call=0x883cdc4, op=0x81fba9c,
args=0x883d5cc,
rho=0x9207e84) at eval.c:1072
#13 0x080a3b28 in Rf_eval (e=0x883cdc4, rho=0x9207e84) at eval.c:418
#14 0x080a3db2 in Rf_applyClosure (call=0x8954ae0, op=0x883d4ec,
arglist=0x9207f64, rho=0x920...
2003 Apr 21
0
sweave provoked segfault (PR#2809)
...0) at
memory.c:1534
#3 0x08097e35 in Rf_duplicate (s=0x8324ff0) at duplicate.c:108
#4 0x080c9abf in GetObject (cptr=0xbfffc2d0) at objects.c:40
#5 0x080ca663 in do_usemethod (call=0x82c5d28, op=0x8209874,
args=0x82c5d44,
env=0x8f1ecd4) at objects.c:401
#6 0x080a3b28 in Rf_eval (e=0x82c5d28, rho=0x8f1ecd4) at eval.c:418
#7 0x080a3db2 in Rf_applyClosure (call=0x8324ff0, op=0x82c5c48,
arglist=0x8d05e80, rho=0x8f385f0, suppliedenv=0x81f65a0) at
eval.c:609
#8 0x080a3937 in Rf_eval (e=0x8324ff0, rho=0x8f385f0) at eval.c:453
#9 0x080a586b in Rf_evalList (el=0x8324f9c, rho=0x8f385f0) at
ev...
2000 Dec 18
1
1.2.0 segfault
...n Rf_allocVector (type=14, length=2) at memory.c:1551
#1 0x5561c in Rf_duplicate (s=0xbcd6f0) at duplicate.c:140
#2 0x55518 in Rf_duplicate (s=0x99c508) at duplicate.c:133
#3 0x55518 in Rf_duplicate (s=0xc63468) at duplicate.c:133
#4 0x5b958 in do_set (call=0xb72804, op=0x1c6304, args=0xb727cc, rho=0x499a38)
at eval.c:1012
#5 0x5a114 in Rf_eval (e=0xb72804, rho=0x499a38) at eval.c:317
#6 0x5b300 in do_begin (call=0xb739ac, op=0x1c6294, args=0xb72820,
rho=0x499a38) at eval.c:803
#7 0x5a114 in Rf_eval (e=0xb739ac, rho=0x499a38) at eval.c:317
#8 0x5a730 in Rf_applyClosure (call=0x49a...
2010 Jan 26
1
newton method for single nonlinear equation
...ch for any help given.
newton.inputsingle <- function(pars,n)
{ runi <- runif(974, min=0, max=1)
lendt <- length(runi)
## Parameter to estimate
z <- vector(length=lendt, mode= "numeric")
z <- pars[1]
## Constant value
alp <- 2.0165 ; rho <- 0.868;
c <- sqrt(pi)/(gamma(alp)*(1-rho)^alp)
for (i in 1:n)
{ t1 <- exp(-pars[1]/(1-rho))
t2 <- (pars[1]*(1-rho)/(2*sqrt(rho)))^(alp-0.5)
bes1 <- besselI(pars[1]*sqrt(rho)/(1-rho),alp-0.5)
bes2 <- besselI(pars...
2004 Sep 13
6
Spare some CPU cycles for testing lme?
If anyone has a few extra CPU cycles to spare,
I'd appreciate it if you could verify a problem that I
have encountered. Run the code
below and tell me if it crashes your R before
completion.
library(lme4)
data(bdf)
dump<-sapply( 1:50000, function(i) {
fm <- lme(langPOST ~ IQ.ver.cen + avg.IQ.ver.cen, data = bdf,
random = ~ IQ.ver.cen | schoolNR);
cat("
2002 May 01
3
bivariate normal cdf and rho
Suppose F(x, y; rho) is the cdf of a bivariate normal distribution, with
standardized marginals and correlation parameter rho. For any fixed x and
y, I wonder if F(x, y; rho) is a monotone increasing function of rho,
i.e., there is a 1 to 1 map from rho to F(x, y; rho).
I explored it using the function pmvnorm in pac...
2001 Jan 23
0
1.2.1 segfault
...emory.c:1670
#2 0x72dbc in Rf_allocVector (type=13, length=4) at memory.c:1542
#3 0x48794 in CPtrToRObj (p=0xd9fb78, arg=0xd9f208, Fort=1) at dotcode.c:180
#4 0x511d0 in do_dotCode (call=0x8ad148, op=0xd9f668, args=0xda163c,
env=0xd9da48) at dotcode.c:1778
#5 0x56c44 in Rf_eval (e=0x3a9fb8, rho=0xd9da48) at eval.c:340
#6 0x58f18 in Rf_DispatchOrEval (call=0x3a94d0, generic=0x1632e8 "[",
args=0x3a94b4, rho=0xd9da48, ans=0xefffc1c4, dropmissing=0) at eval.c:1386
#7 0xbea40 in do_subset (call=0x3a94d0, op=0x1c3754, args=0x3a94b4,
rho=0xd9da48) at subset.c:471
#8 0x56aa0...
2011 Oct 19
1
Sparse covariance estimation (via glasso) shrinking to a "nonzero" constant
...hers have different
nonzero values", I'd like it to be "let's restrict some of these
elements are 0.3 and others are different or higher/lower than this
0.3". Does that make any sense?
Here is the glasso code if it helps:
> library(glasso)
> glasso
function (s, rho, zero = NULL, thr = 1e-04, maxit = 10000, approx = FALSE,
penalize.diagonal = TRUE, start = c("cold", "warm"), w.init = NULL,
wi.init = NULL, trace = FALSE)
{
n = nrow(s)
BIG = 1e+10
if (!is.matrix(rho) & length(rho) != 1 & length(rho) != nrow(s)...
2002 Apr 15
3
Greek in text()
I have gone over the examples and can't figure this out:
rho<-.77
text(x=.05,y=.5,paste(expression(rho),rho))
I was hoping to get this to print a Greek rho with 0.77 beside it.
Instead I get: rho 0.77 (i.e. Roman lettering)
The help on expression() is quite opaque so I don't understand how it
works.
Thanks for any help.
Bill Simpson
-.-.-.-.-.-.-....
2003 Apr 05
0
sweave/gc segfault bt
...xt/R/cvs/R-devel/src/main/gram.y:1147
#10 0x08117b61 in do_parse (call=0x883dc90, op=0x8205704,
args=0x88397c0,
env=0x94f5cf0) at source.c:76
#11 0x080c9e24 in do_internal (call=0x883dc04, op=0x81fbcc8,
args=0x94f52d4,
env=0x94f5cf0) at names.c:1056
#12 0x080a3d18 in Rf_eval (e=0x883dc04, rho=0x94f5cf0) at eval.c:418
#13 0x080a5172 in do_begin (call=0x883d3e0, op=0x81fc17c,
args=0x883dbe8,
rho=0x94f5cf0) at eval.c:1072
#14 0x080a3d18 in Rf_eval (e=0x883d3e0, rho=0x94f5cf0) at eval.c:418
#15 0x080a3fa5 in Rf_applyClosure (call=0x8955188, op=0x883db08,
arglist=0x94f5bbc, rho=0x9...
2008 Oct 03
1
Memory crash
...at memory.c:1317
#1 0x000000000041c31c in Rf_allocVector (type=14, length=3000)
at memory.c:1968
#2 0x000000000052a647 in duplicate1 (s=0x20cb4a0) at duplicate.c:214
#3 0x00000000004e1775 in do_makelist (call=<value optimized out>,
op=<value optimized out>, args=0x1e18238, rho=<value optimized out>)
at builtin.c:645
#4 0x000000000054697a in Rf_eval (e=0x10685a0, rho=0x1c557e0) at eval.c:492
#5 0x000000000054a7ef in do_set (call=0x10684f8, op=0x936188,
args=0x1068530,
rho=0x1c557e0) at eval.c:1424
#6 0x0000000000546896 in Rf_eval (e=0x10684f8, rho=0x1c...
2013 May 13
3
help: R GUI front-end has stopped working
...r (i in 1:n){
v=as.numeric(u<=u[i])
p=mean(v)
A=vector(length=k)
n1=n%/%k
A<-sapply(1:k,function(m) A<-sum(v[((m-1)*n1+1):(m*n1)]))
A[k]=sum(v[((k-1)*n1+1):n])
AA[[i]]<-A
B=c(rep(n1,k-1),length(v[((k-1)*n1+1):n]))
BB[[i]]<-B
l1=vector(length=(length(A)))
L1<-function(pe,rho){
for(j in 1:length(A)){
if(A[j]==0){l1[j]=0} else {l1[j]<-sum(log(pe+(-rho/(rho-1))*(0:(A[j]-1))))}
}
return(sum(l1))}
l2=vector(length=(length(B-A)))
L2<-function(pe,rho){
for(j in 1:length(B-A)){
if((B[j]-A[j])==0){l2[j]=0} else {l2[j]<-sum(log(1-pe+(-rho/(rho-1))*(0:(B[j]-A[j]-...
1999 Jul 02
2
segfault on Sparc Linux (RH6.0)
...es
Program received signal SIGSEGV, Segmentation fault.
0xb1dc8 in hqr2_ ()
(gdb)
bt gives:-
(gdb) bt
#0 0xb1dc8 in hqr2_ ()
#1 0xb4bb4 in rg_ ()
#2 0x3909c in do_dotCode (call=0x70421f64, op=0x70421b90, args=0x70679f28,
env=0x7067791c) at dotcode.c:386
#3 0x3fc7c in eval (e=0x70421f64, rho=0x7067791c) at eval.c:180
#4 0x41118 in do_set (call=0x70421fa0, op=0x702d4828, args=0x70421f8c,
rho=0x7067791c) at eval.c:761
#5 0x3fc7c in eval (e=0x70421fa0, rho=0x7067791c) at eval.c:180
#6 0x40b58 in do_begin (call=0x70421fb4, op=0x702d48c8, args=0x70421fc8,
rho=0x7067791c) at eva...
2007 Aug 01
0
a question about d F(x, y; rho) / d rho = f(x, y; rho)
Hello,
My name is IKJIN LEE and I am studying in University of Iowa.
Recently, I have been trying to prove the fact d F(x, y; rho) / d rho =
f(x, y; rho),
but I could not. In the mean time, I found that Professor S. Le Cessie
from Netherlands kindly gave you an elegant proof of
d F(x, y; rho) / d rho = f(x, y; rho).
If you don't mind, can I know the proof?
IKJIN LEE
2007 Mar 27
3
Use of 'defineVar' and 'install' in .Call
Dear all,
[system and version information below]
I am trying to modify a C function for finding the root of an
expression. The function is to be called from R as .Call with input
parameters:
f: expression for which we will find the root
guesses: interval for the solution
stol: tolerance
rho: environment
The original functions I use are:
SEXP mkans(double x) {
SEXP ans;
PROTECT(ans = allocVector(REALSXP, 1));
REAL(ans)[0] = x;
UNPROTECT(1);
return ans;
}
double feval(double x, SEXP f, SEXP rho) {
defineVar(install("x"), mkans(x), rho);
return(REAL(eval(f, rho)...
2007 Mar 27
3
Use of 'defineVar' and 'install' in .Call
Dear all,
[system and version information below]
I am trying to modify a C function for finding the root of an
expression. The function is to be called from R as .Call with input
parameters:
f: expression for which we will find the root
guesses: interval for the solution
stol: tolerance
rho: environment
The original functions I use are:
SEXP mkans(double x) {
SEXP ans;
PROTECT(ans = allocVector(REALSXP, 1));
REAL(ans)[0] = x;
UNPROTECT(1);
return ans;
}
double feval(double x, SEXP f, SEXP rho) {
defineVar(install("x"), mkans(x), rho);
return(REAL(eval(f, rho)...
2012 Apr 24
2
Some Help Needed
Dear all,
I need to do some calculation where the code used are below. I get
error message when I choose k to be large, say greater than 25.
The error message is
"Error in integrate(temp, lower = 0, upper = 1, k, x, rho, m) :
the integral is probably divergent".
Can anyone give some help on resolving this. Thanks.
Hannah
m <- 100
alpha <- 0.05
rho <- 0.1
F0 <- function(y,x,rho){
pnorm((qnorm(x, lower.tail = F)-sqrt(rho)*qnorm(y, lower.tail = F))/sqrt(1-
rho), lower.tail = F)
}...
2008 Feb 19
1
recursive function help
...spect I need a Vectorize somewhere,
but I can't suss it out. Any help would be appreciated. I've tried
traceback() and various debugging ideas to no avail (most likely due to
my inexperience with these tools.)
Here's what I have.
Nk <- function(m, C) {
if (length(m) > 1) {
rho <- C[1, -1]
Rmat <- C[-1, -1]
B <- diag(1/sqrt(1 - rho*rho)) %*%
(-rho %*% t(rho) + Rmat) %*%
diag(1/sqrt(1 - rho*rho))
integrate( function(x) dnorm(x) * Nk((m[-1] - rho*x)/sqrt(1 -
rho*rho), B), -10, m[1] )$value
} else {
pnorm(m[1])
}
}
my example i...