similar to: R-alpha: random number generator -- S-plus's

I have been using delay() quite successfully to 'load on reference' some large data files (earthquake catalogues). One minor difficulty I have not been able to overcome is that upon first reference, the data set is loaded into memory as expected, but also remains as a permanent object in the user's environment, instead of appearing to be only in the library from which it was loaded.

>>>>> Thomas Lumley writes: >> Two compatibility issues found while trying to convert a simulation >> from S to R. >> 1. set.seed() We don't have this function. According to Venables & >> Ripley it just picks a seed from a list of 1000 possibilities. How about >> >> "set.seed" <-function (i)

Hi, Is the below -- setSeed -- an okay way to set the seed in standalone applications making use of Rlib? It seems to work as expected. Is there a better way to do it? (I'm also looking at do_setseed but am unsure what to supply as op. findFun("set.seed", R_GlobalEnv)?) Thanks much. -- David Faden, dfaden at iastate.edu AIM: pitulx -- #include <assert.h> #include

In 0.65.0 lib.source() within library() works differently and fails to correctly "load" a library when it wants to change an object that it has defined previously: Below is a very short shar file that produces a "minimal" example [no attachment, just appended]. To reproduce {read on, before doing it!} do sh < libBug.shar <<<<< the file name you gave

Included below is a shar archive containing two patches that together: 1) make backup files get their setuid and setgid bits stripped by default 2) add a "-s" option that allows backup files to continue to have these privileges This means that if you update a collection of binaries with rsync, and one or more of them has a local-root security problem, the backup file(s) created when

Hi All, I'm new to R and am trying to run a unit root test on the vector "y" (a time series of inflation (i.e. changes in the Consumer Price Index quarter on quarter)). I've run the Augmented-Dickey-Fuller Test below (R's URCA package). It gives me an error that it cannot find the function ur.df unless I comment out the third last line of code (see below). I try to call

I am not a member of the samba listserv, but I wanted to contribute the following code to the samba effort. The following is clipped from some experimental changes I have made to my own copy of the samba source. I have not included all of the changes because I don't have the time to. But someone should be able to properly do the integration without too much effort. Summary, these

Dear contributors I have tried this experiment: x<-c() for (i in 1:12){ x[i]<-list(cbind(x1[i],x2[i])) #this is a list of 12 couples of time series I am using to perform a test } # that compares them 2 by 2 # ################# #trace statistic test<-data.frame() cval<-array( , dim=c(2,3,12)) for (i in 2:12){ for (k in 1:2){ for (j in 1:3){ result[k,j,i]<-

Hello everyone, I'm trying to estimate the parameters of the returns series attached using the GARCH code below, but I get the following error message: Error in solve.default(Hessian) : system is computationally singular: reciprocal condition number = 0 Error in diag(solve(Hessian)) : error in evaluating the argument 'x' in selecting a method for function 'diag' Can

If you found my subject heading to be confusing then I'm sure you'll enjoy the example I've included below. I find the apply type functions to be wonderful for avoiding loops but when I use them with existing functions, I end up using loops anyway to extract the vectors I want. I would appreciate it if someone could show me how to avoid these loops. Thanks. EXAMPLE:

There's a little source code (typo?) bug in "samba-2.0.4/source/lib/system.c" that hurts compilation at least on HPUX 9.0x (with HP's cc, not gcc!) !! In HPUX's "/usr/include/sys/time.h": .. # ifndef _STRUCT_TIMEVAL # define _STRUCT_TIMEVAL /* Structure returned by gettimeofday(2) system call and others */ struct timeval { unsigned long tv_sec;

Dear R users, I came up with some simple functions to give me the standard betas and tolerance values from a predefined lm() model. I have been trying to insert the results from these functions into the coefficients matrix in a modified summary.lm function that I'm calling summary2 (I'd never edit the summary.lm function directly!!). I managed to get the results inserted into the

I am attempting to lower the selectCC instruction to the instruction set of the backend I'm working on and I cannot seem to find a way to correctly implement this instruction. I know how this instruction should get implemented; I just have yet to find a way to do it. I want the select_cc instruction to be lowered into a comparison followed by a conditional move. I've attempted to use a

Hello all, I'm trying to create a heatmap using 2 matrices I have: z and v. Both matrices represent different correlations for the same independent variables. The problem I have is that I wish to have the values from matrix z to be represented by color intensity while having the values from matrix v to be represented by circle size. I currently have the following in front of me and an unsure

Hello, I am using smbspool and I seem to be getting an error that I can't find a solution to. I checked Google, Usenet and the samba archives. I seem to find people with the problem, but no answer. When I try to use the smbspool program, I get the following error: > root@HomeLinux:/home/steve# smbspool smb://192.168.1.4/LaserPrinter 0 0 0 0 0 README > ERROR: cli_session_request()

My question is about the "Signif. codes" , the output when I run matcoef =cbind(fit$par, se.coef,tval,2*(1-pnorm(abs(tval)))) dimnames(matcoef)=list(names(tval),c("Estimate","Std.Error","t value","pr(>|t|)")) cat("\nCoefficient(s):\n") printCoefmat(matcoef, digits=4, signif.stars = TRUE) Coefficient(s): Estimate

Hi, pretty good job you are doing with Samba... I have just downloaded Samba 2.0.4 and I have a problem trying to compile it. I have configured the "configure" file and ran it (./configure). When I type make, I can see all the different modules compiling but I have an error with the /source/lib/system.c file. It is concerning the timeval struct (which I concluded must be a system

I define the following function to convert a t-value with degrees of freedom DF to another t-value with different degrees of freedom fullDF: tConvert <- function(tval, DF, fullDF) ifelse(DF>=1, qt(pt(tval, DF), fullDF), 0) It works as expected with the following case: > tConvert(c(2,3), c(10,12), 12) [1] 1.961905 3.000000 However, it gives me warning for the example below although the

The following code causes opt to not terminate! With TOT this morning, and of a week ago: clang foo.c and clang -O1 foo.c work fine. clang -O2 foo.c and clang -O3 foo.c do not terminate. (At least after 10 minutes) If I generate the bit code (clang-cc -emit-llvmbc) and then run: opt -O3 foo.bc it does not terminate. //foo.c int get_id(int); typedef short

I really can't figure out how to get and set array indices from the c interface. so to get an element I'm calling tindex = *fn\SymbolTable(*index\name) index = LLVMBuildLoad(builder,tindex,"index") arr = *fn\SymbolTable(*array\name) arrptr = LLVMBuildLoad(Builder,arr,"arrayptr") tmp = LLVMBuildGEP(Builder,arrptr,index,0,"ptr") ptr =