similar to: Sparse covariance estimation (via glasso) shrinking to a "nonzero" constant

Hi all, I am trying to compile a fortran routine originally written for R from within matlab. It all worked fine, but the newest version of the Fortran code has a call to rchkusr() and I now get a message that this is an undefined symbol. I am under the impression that rchkusr() just serves to let the processor take care of interruption messages (typically, the user frantically

I could swear this worked fine earlier.. Now facter, and hence puppet, do not see the domain name anymore: What could cause this? C:\Temp>facter --version 1.6.2 C:\Temp>facter | findstr kernel kernel => windows kernelmajversion => 6.1 kernelrelease => 6.1.7601 kernelversion => 6.1.7601 C:\Temp>facter fqdn NCS-VDI-05.none C:\Temp>facter hostname NCS-VDI-05

The following example in ?dump.frames options(error = quote({dump.frames(to.file = TRUE); q()})) is useful for teaching the user how to save a frame dump when R encounters an error during non-interactive sessions. This command however causes an additional change that on encountering an error R exits with a 0 error status. Although it's just an example, it's an important one as it's

This is useful for probing probing for cache files such as: external_command ?ignore_error:(Some true) ~prog "ls .cache/something.*" will return command output (matched files) on its success or empty list whenits exit code is other than 0 (there are no such files). --- mllib/common_utils.ml | 15 ++++++++++----- mllib/common_utils.mli | 2 +- 2 files changed, 11 insertions(+), 6

Hi, Is there a function func(x) where x is an array such that it returns a list of the indices of all non-zero elements of the array? for example: a: 1 0 0 0 0 1 0 1 0 func(a) returns two vectors r and c: r c 1 1 2 3 3 2 daver +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ |David Richmond It works on a | + Dept. of Sociology complex scientific

Dear Prof. Roger Koenker, On Tue, 22 Oct 2024 09:08:12 +0000 "Koenker, Roger W" <rkoenker at illinois.edu> wrote: > > fN <- rqss(y~qss(x,constraint="N")+z) > > *** caught segfault *** > address 0x0, cause 'invalid permissions? Given a freshly produced quantreg.Rcheck directory, I was able to reproduce this crash by running R -d gdb # make

Can someone verify for me if the for loop below is really calculating the nonzero min for each row of a matrix? I have a bug somewhere in the is section of code. My first guess is how I am find the the nonzero min of each row of my matrix. The overall idea is to make sure I am investing all of my money, i.e. new.set is a set of indicator variables for each stock for a particular portfolio, i.e.

Gurus: I have a new version of my quantreg package with minimal changes, mainly to fix some obscure fortran problems. It fails R CMD check ?as-cran with the error: Running examples in ?quantreg-Ex.R? failed The error most likely occurred in: > base::assign(".ptime", proc.time(), pos = "CheckExEnv") > ### Name: plot.rqss > ### Title: Plot Method for rqss Objects

Hi All, 1) Is it possible to set the options such that R opens a new script editor every time I start the R and 2) specify the size of windows. Thanks for the suggestion and Best regards, Krishna [[alternative HTML version deleted]]

Dear all, I am still trying to find out what goes wrong when I use ace (I get an 'NAs in foreign fuction call (arg3)' but when I go through it with debug there seem no data to be missing). I had some errors (see at end of message) while installing 'acepack'. Do you think it is possible that these lead to the fact, that simple problems run ok but larger problems result in an

Hello! I am performing coupling of chains in MCMC and I need the same value of seed for two chains. I will show demo of what I want: R code, which might show my example is: niter <- 3 nchain <- 2 tmpSeed <- 123 for (i in 1:niter) { # iterations for (j in 1:nchain) { # chains set.seed(tmpSeed) a <- runif(1) cat("iter:", i, "chain:", j,

Full_Name: Russell Senior Version: 0.64.0 OS: linux-2.2.x Submission from: (NULL) (198.107.133.11) It appears that the GNU readline history size is limited to 50 entries. This is apparently occurring because R is not bothering to override the readline default value (from readline-4.0): ./history.c:59:#define DEFAULT_HISTORY_GROW_SIZE 50 It appears that it ought to be possible to override

Hi, I'm trying to compile two functions into one function. the first funtion is called 'fs' which is self-made function, another function is from the built-in 'integration' function that is copy-paste-edited. If built separatey, these functions work well. However that is not the case if combines together, where certainly I made mistake somewhere when constructing the code.

Hi list, It appears there is a problem with the OpenMPI I/O library on CentOS 6.2 & 6.3 (package openmpi-1.5.4-1.el6.x86_64). When I compile the attached program it ends up in the error path since MPI_File_open returns 16. The corresponding (unhelpful) message is: MPI_ERR_OTHER: known error not in list I couldn't find any pointers on the net and the same program works with OpenMPI

The following patch adds a trivial family method for linear models. -k --- src/library/base/funs/lm.orig Tue May 27 02:53:58 1997 +++ src/library/base/funs/lm Sun Aug 10 22:00:46 1997 @@ -412,3 +412,5 @@ } formula.lm<-function(x)formula(x$terms) + +family.lm <- function(x) { gaussian() } =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- r-devel mailing list

Dear R users, I'm pretty new on using lmer package. My response is binary and I have fixed treatment effect (2 treatments) and random center effect (7 centers). I want to test the effect of treatment by fitting 2 models: Model 1: center effect (random) only Model 2: trt (fixed) + center (random) + trt*center interaction. Then, I want to compare these 2 models with Likelihood Ratio Test.

Hi all, I am currently doing a consistency analysis using an lmer model and trying to use parametric bootstrapping for the confidence intervals. My model is like this: model<-lmer(y~A+B+(1|C/D)+(1|E),binomial) where E is the individual level for consistency analysis, A-D are other fixed and random effects that I have to control for. Following Nakagawa and Scheilzeth I can work out the

Hello, I would appreciate *any* ideas on this problem. I'm the maintainer of a package ("subselect"), which on CRAN's Daily Package Checks is OK on all flavours of R, except r-devel Linux x86_64, where there is a "memory not mapped" segfault with the very first example that is tried out (output below). Additionally, a user with an AMD64 machine has just reported a

I don't know the details of pls (in the pls.pcr package, I assume), but if you use validation="CV", that says you want to use CV to select the best number of components. Then why would you specify ncomp as well? Andy > From: ryszard.czerminski at pharma.novartis.com > > When I try to use ncomp parameter in pls procedure I get > following error: > > >

Hello, I'm working on a problem using nested for-loops and I don't know if it's a problem with the order of the loops or something within the loop so any help with the problem would be appreciated. To briefly set up the problem. I have 259 trees (from 11 different species, of unequal count for each species) of which I am trying to predict biomass. For each tree species I have 10000