search for: nfac

R 1.4.1, Linux and Windows for(i in 1:1000) { X <- rnorm(i) XX <- fft(fft(X), inverse=T)/i if(max(Mod(XX-X)) > 1e-10) print(i) } [1] 392 [1] 588 [1] 968 [1] 980 and I then get a segfault during gc(). The answers are way off, with imaginary parts 1e10 or more. These numbers are all multiples of 7^2 or 11^2. (Based on a report to R-help Date: Thu, 28 Mar 2002 09:37:34

...causes a segmentation fault or gets stuck in a loop R> mvfft(xx) Program received signal SIGSEGV, Segmentation fault. 0x80d36cb in fftmx (a=0x40617eb0, b=0x40617eb8, ntot=30, n=30, nspan=30, isn=-2, m=3, kt=0, at=0x407c8f58, ck=0x407c8f80, bt=0x407c8fa8, sk=0x407c8fd0, np=0x407c8ee4, nfac=0x8115a78) at fft.c:558 558 bt[j] = b[k1] - b[k2]; (gdb) print bt[j] Cannot access memory at address 0x407cb000. (gdb) print kspan $4 = 0 We're in a loop here where kspan is the increment for k1 and decrement for k2. The loop breaks when k1 < k2, but somehow kspan is zero. How it g...

There is a problem with 'default argument evaluation' when I use an existing function name as argument name : sintest <- function(x, y = 2, sin= sin(pi/4)) { ## Purpose: Test of "default argument evaluation" ## -------- Fails for R-0.49. Martin Maechler, Date: 9 May 97. c(x=x, y=y, sin=sin) } ## R-0.49: R> sintest(1) ##> Error in sintest(1) : recursive

...," factor names ") dostop<-FALSE # Ensure the return value is defined for (i in seq(along=nn)) { message("Factor ",nn[i]," in estimation data") ofac <- object[["model"]][[nn[i]]] print(ofac) # diagnostic oldlev<-levels(factor(ofac)) nfac <- newdata[[nn[i]]] # ?? note need for double brackets. Why? message("Factor ",nn[i], " in new data") print(nfac) # diagnostic message("about to try levels on newfac") newlev<-levels(factor(nfac)) for (j in seq(along=newlev)) { if (!(newlev[j]...

Hi Developers, I am looking for another new project to help me get more up to speed on R and to learn something outside of R internals. One recent R issue I have run into is finding a fast implementations of the equivalent to the following SAS code: /* MDPC is an integer sort key made from two integer columns */ MDPC = (MD * 100000) + PCO; /* sort the dataset by the key */ PROC SORT;

...n", envir = parent.frame()) if (is.function(family)) family <- family() if (is.null(family$family)) { print(family) stop("'family' not recognized") } fltype <- mkFltype(family) FL <- lmerFactorList(formula, mf, fltype) nFacs <- length(FL) #My insert cnames <- with(FL, c(lapply(Ztl, rownames), list(.fixed = colnames(X)))) nc <- with(FL, sapply(Ztl, nrow)) Ztl <- with(FL, .Call(Ztl_sparse, fl, Ztl)) Ztl[[nFacs]] <- t(rand_matI)...

...causes a segmentation fault or gets stuck in a loop R> mvfft(xx) Program received signal SIGSEGV, Segmentation fault. 0x80d36cb in fftmx (a=0x40617eb0, b=0x40617eb8, ntot=30, n=30, nspan=30, isn=-2, m=3, kt=0, at=0x407c8f58, ck=0x407c8f80, bt=0x407c8fa8, sk=0x407c8fd0, np=0x407c8ee4, nfac=0x8115a78) at fft.c:558 558 bt[j] = b[k1] - b[k2]; (gdb) print bt[j] Cannot access memory at address 0x407cb000. (gdb) print kspan $4 = 0 We're in a loop here where kspan is the increment for k1 and decrement for k2. The loop breaks when k1 < k2, but somehow kspan is zero. How it g...