similar to: regression methods for circular(?) data.

Dear R-developers, Using .Call I pass a S4 class with e.g. the following class definition: setClass("mmatrix",representation( data="matrix") ) On the "C side" i do mat = GET_SLOT(vs,install("data")); and then: printf("%f\n",NUMERIC_POINTER(mat)[1]); The above print statement produces the correct output if xx<- new("mmatrix")

Dear R-helpers, The function optim implements algorithms that I would like to use. I have function implemented in R, which given the parameters of which minimization is to take place returns a scalar as well as the gradient. Unfortunately optim requires two function _fn_ and _gr_ where fn returns the function value and gr the gradient. Splitting my function in two functions would be easy,

Hi, How to add values to an array at any position. Asking because of the following: e.g. y<-c(0.1,NaN,0.2,NaN) #or data frame x<-na.omit(y) take some columns from x and do some computation with functions which do not allow NaN 's. After the computing add NaN's at positions stored in attr(x,"na.action") of the result vector. /E -- Dipl. bio-chem. Witold Eryk

Dear R-developers, I am working on an C interface to some c functions. I compiled the c file using R CMD SHLIB without any compilation errors. All libs used are added to the LD_LIBRARY_PATH. However dyn.load("interface.so") produces undefined symbol: N_VNew_Serial; which is defined in one of the libs I link against. Where, and how to start searching for the error? Eryk

Dear R-developers, I am getting the following error message (see bottom of e-mail) and can't make a clue out of it. The zzz.R file contains .First.lib <- function(lib, pkg) library.dynam("SBMLodeSolveR",pkg,lib) .Last.lib <- function(libpath) library.dynam.unload("SBMLodeSolveR", libpath) There is a SBMLodeSolveR.so file in the src directory which works if loaded

Dear R-developers, I am trying to reproduce the autoconf.ac example from R-ext and fail. My autoconf file looks like this [autoconf.ac] # original by Friedrich Leisch, much changed by BDR AC_INIT([SBMLodeSolveR]) dnl Select an optional include path, from a configure option dnl or from an environment variable. AC_ARG_WITH([sbmlode-include],

Hello, I am almost new in R, so perhaps my question will be silly. I try to use R for analyzing capture-recapture data in epidemiology. A cancer registry has different sources of patients. We know in each list, patients already known in all other list. The aim is to use capture-recapture models for estimating the number of patients unknow of all the sources. Because no order in sources, one

Normal 0 21 false false false MicrosoftInternetExplorer4 Dear R-helpers,   I was wondering whether there are any good books and/or website links that introduce mark-recapture analysis in R. In

Hi Does R have a proceedure/software for capture recapture? Thank you. Lawrence Lessner Please respond to LLessner at nycap.rr.com

Dear R helpers, I generate a qq plot using the following function call. qqmath(~val|ind,data=xx ,distribution=function(p) qt(p,df=19) ,ylab="Sample Quatinles" ,xlab="Theoretical Quantiles" ,aspect=1 ,prepanel = prepanel.qqmathline ,panel=function(x,y) { panel.qqmathline(y, distribution=function(p) qt(p,df=19),col=2)

Dear R-devlopers, I removed yesterday the R.2.1. installation from windows NT laptop and I am getting the following error when building a pacage with R.2.2 ---------- Making package msbase ------------ adding build stamp to DESCRIPTION making DLL ... make[3]: *** No rule to make target `d:/prog/r/rw2011/include/R.h', needed by `llogic.o'. Stop. .... I have no idea from where

Hello, I have a couple questions regarding generalized linear mixed models specifically around fitting the random effects terms correctly to account for any pseudo-replication. I am reading through and trying to follow examples from Zuur et al. Mixed Effects Models and Extensions in Ecology with R, but am still at bit unsure if I am specifying the models correctly. Background information: Our

Dear R gurus. Is there an R function equivalent to octaves / (Right division) withouth forming the inverse of Y' using solve ? [snip - from octave docu] Right division. This is conceptually equivalent to the expression (inverse (y') * x')' but it is computed without forming the inverse of Y'. If the system is not square, or if the coefficient matrix

Hi All, My kit : Centos4.4 + Postfix + SA + ClamAV (not that its relevant) General question - Does anyone know of a product which can take a folder / mbox of mail (spam, ham, virus infected stuff, scam mail and hoaxes) and "replay" it back against a test server to determine the servers ability to filter? I have looked at some low level (IP) packet replayers, but i cant really find

New function below is a bit faster due to more efficent memory handling. for-loop FTW! directionless_circular_permutations2 <- function( n ) { n1 <- n - 1L v <- seq.int( n1 ) ix <- combinations( n1, 2L ) jx <- permutations( n-3L, n-3L ) jxrows <- nrow( jx ) jxoffsets <- seq.int( jxrows ) result <- matrix( n, nrow = factorial( n1 )/2L, ncol = n ) k

Jeff, I wanted to let you know that your function is faster than generating the directional circular permutations and weeding. Here is the time for n = 10. I compared with just doing the permutations, there is no point in proceeding further with the weeding since it is slower at the start itself. system.time(directionless_circular_permutations(10)) user system elapsed 1.576 0.000

I don't know if this is more efficient than enumerating with distinct directions and weeding... it seems kind of heavyweight to me: ####### library(gtools) directionless_circular_permutations <- function( n ) { v <- seq.int( n-1 ) ix <- combinations( n-1, 2 ) jx <- permutations( n-3, n-3 ) x <- lapply( seq.int( nrow( ix ) ) , function( i ) {

Thanks! Yes, however, this seems a bit wasteful. Just wondering if there are other, more efficient options possible. Best wishes, Ranjan On Thu, 29 Mar 2018 22:20:19 -0400 Boris Steipe <boris.steipe at utoronto.ca> wrote: > If one is equal to the reverse of another, keep only one of the pair. > > B. > > > > > On Mar 29, 2018, at 9:48 PM, Ranjan Maitra

Hi, While checking the package I am getting all the time the following complain about the function: Undocumented S4 methods: generic coerce and siglist matrix,Massvector I have not defined the function nowhere in the R files. My gues is that this has something to do with the following definition of the setAs function. setAs("matrix","Massvector" ,def=

Hi! I am trying to calculate a dissimilarity matrix using daisy. The matrix vectver is binary as i test with: > levels(as.factor(vectver)) [1] "0" "1" But the call to daisy gives me the following error message.: > dfl1 <- daisy(vectver, type = list(asymm = c(1:length(vectver[,1])))) Error in daisy(vectver, type = list(asymm = c(1:length(vectver[, 1])))) : at least