similar to: bus error macosx/off-topic

Hi all, I am working on Macos x 10.3 (Panther) to build a package consisting of C/C++ code that is called from R. In the C/C++-sources I use several commands to print info to the console: std::cout << "info" << endl; and: Rprintf("info\n"); Both work fine when R is run on the command line but neither works when running Raqau (I did check the preference boxes

Hello All, In some package I use a c-routine which calls a fortran routine which expects a char-string as input. As per the writing R-extensions manual, the Fortran routine is declared in C as: void F77_NAME (setoptions) (char **option); and then it is calles as follows: char **option; option = new char*[1]; option[0] = new char[256]; option[0] = strcpy(option[0],"Iteration

Hi All, Does someone have an idea of how to cleverly convert a categorical timeseries into a transition matrix? Ie, I have something like: x<- c(1,1,2,1,1,2,2,2,1,2), And I want a matrix with counts and/or probabilities: > tr <- matrix(c(2,3,2,2),2,2) > tr [,1] [,2] [1,] 2 2 [2,] 3 2 Meaning that there are two transitions from 1 to 1, two from 1 to 2, three from 2 to 1

-- R Help List -- I am looking for a bivariate normal cdf routine in R. I have some fortran routines for this, which appear to be based on 15-point quadrature. Any guidance/suggestions on making these in loadable R-functions would be appreciated. Thanks, Dan =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Daniel A. Powers, Ph.D. Department of Sociology University of Texas at Austin

Dear useRs, I'm trying to create a barplot like so: x=matrix(1:10,2,5) barplot(x,leg=c("left","right"),besid=T) The legend is placed in default position topright, however the data are plotted there too. I tried controlling the legend position by adding x="topleft" but this results in an error that x matches multiple formal arguments. Leaving out the legend

Is there a package in R that handles general linear (in-)equality + box constrained optimization? If it is not there, could anyone advise me which way to go? And/or point me to packages that solve these problems partially? best, ingmar -- Ingmar Visser Department of Psychology, University of Amsterdam Roetersstraat 15, 1018 WB Amsterdam The Netherlands http://users.fmg.uva.nl/ivisser/ tel:

Hello all, I am using an existing Fortran routine that takes a single character string as argument. The routine echoes the argument that I provide. When working on OS X 3.9 there seems to be no problem, ie the Fortran routine nicely echoes my argument. However, I compiled the same package on a PC (using all the tools provided in the R for windows faq), and the routine only echoes the first letter

Dear R-help, I'm trying to load the fGarch package and keep running into problems with tcltk: After succesfully instaling fGarch (and dependencies) I get: >library(fGarch) Loading required package: fBasics Loading required package: fImport Loading required package: fSeries Loading required package: robustbase Loading required package: fCalendar Loading required package: MASS Loading

Dear All, When installing the Ruuid package (from Bioconductor) from sources on my MAC (OS X 10.3.5, R version 2.0.0) I get the following errors: m00245:~ ivisser$ R CMD INSTALL -l /Users/ivisser/Library/R/library/ /Users/ivisser/Desktop/Ruuid * Installing *source* package 'Ruuid' ... loading cache ./config.cache checking for glib-config... no checking how to run the C preprocessor...

Hello, I tried to fit a poisson distribution but looking at the function fitdistr() it does not optimize lambda but simply estimates the mean of the data and returns it as lambda. I'm a bit confused because I was expecting an optimization of this parameter to gain a good fit... If I would use mle() of stats4 package or mle2() of bbmle package, I would have to write the function by myself

Dear R-help, I'm facing a problem with defining a repeated measures anova with weighted data. Here's the code to reproduce the problem: # generate some data seed=11 rtrep <- data.frame(rt=rnorm(100),ti=rep(1:5,20),subj=gl (20,5,100),we=runif(100)) # model with within factor for subjects/repeated measurements, no problem aov(rt~ti + Error(subj/ti),data=rtrep) #model with weights

Hello All, Does anyone know of a package that performs constraint-based clusters? Ideally the package could perform "Time-Ordered Clustering", a technique applied in a recent journal article by Runger, Nelson, Harnish (using MS Excel). Quote, "in our specific implementation of constrained clustering, the clustering algorithm remains agglomerative and hierarchical, but observations

Hi all, I'm writing a vignette for my package, and I would like to include some of the package help files in there as well. Is there an easy way of doing so? I tried using R CMD Rdconv to generate latex files from .Rd files but I am not sure how to include these into a .Rnw file (ie the vignette source). The resulting file from Rdconv do not readily compile using latex ... The other option I

I use the lme function from the nlme library (or alternatively from the Matrix library) to estimate a random effects model. Both functions return the covariance matrix of the estimated parameters. I have the following question: Is it possible to retrieve the information matrix of such a model (ie from the fitted object)? In particular, the information matrix can be computed as a sum of individual

Hello list, I am dealing with a noisy function (gradient,hessian not available) with simple boundary constraints (x_i>0). I've tried constrOptim() using nelder mead to minimize it but it is way too slow and the returned results are not satisfying. simulated annealing is so hard to tune and it always crashes R program in my case. I wonder if there are any packages or functions can do

Hi My data looks like this startDate="2008-06-01" dateRange =c( "2008-10-01","2008-12-01") Is there any method to find the week number from the startDate range ----- Thanks in Advance Arun -- View this message in context: http://r.789695.n4.nabble.com/Week-number-from-a-date-tp4410223p4410223.html Sent from the R help mailing list archive at Nabble.com.

Can anybody, please, explain me how many parameter are estimated using randomLCA? For examples, model "dentistry.lca2random" estimate 1 scale (or variance, b_j) parameter and 2 position parameters (a_cj)? Doesn't it? Do I need at least 4 diagnostic tests for such a model? What happens if I specify options blocksize and byclass? How many diagnostic tests (or rater) I need?

Dear List, We are using RHmm to cluster data through HMM. We would like to restrict the transition matrix of HMM, to get hierarchical connections between clusters. But, RHmm doesn't seem to support these restrictions. Can any one suggest a library to do that. Thanks, Kishor

I'm trying to get the this Fortran array bounds check to run on my code as my package has been rejected for failing this test. I read writing R extensions and it appears to run this test I should use a Makevars file with one of the following two flags depending on my compiler: FCFLAGS = -g -O2 -mtune=native -fbounds-check FFLAGS = -g -O2 -mtune=native -fbounds-check and store it in the src

I had written a vignette and included a \usepackage{Rd} command to make it possible to include latex'ed Rd files in the vignette. However, when loading Rd.sty texShop produces the following error: l. 180 ...d}[1]{\ifmmode\bm{#1}\else\textbf{#1}\fi} This is on Max OS X 3.9 and R 2.0.1 Has anyone seen this before and/or is it problematic? I'm not sure whether the output suffers from this