similar to: Using R in C++

Hallo, I didn´t found any facilities for Halbergs cosinor analysis in R. This analysis is well known in the Chronobiology as the least square approximation of time series using cosine function of known period (in my case of 24hours-period). I tried to write a script but crashed... Can you give me some advices, please!? Thanks Anne Berger Institute of Zoo- and Wildlife Research, Berlin, Germany

Thanks a lot to Demitris for a prompt answer some minutes ago on another tread (see below). To avoid excess mails on the list, I move onto next question: I have another small plotting problem that confuses me. I want to plot results from a field trial series, using the numbers of the trials as symbols in the plot. pch = as.character(trial_no) works fine, but truncates the trial number to the

Dear R-users, Just like other users (as seen from previous posts on the list), I would like to use female and male signs in plots. I found B. Ripley's post about using Unicode characters. However, it doesn't works for me. > text(locator(1),"\u2640") produces the following error: Error: invalid \uxxxx sequence But I can specify other Unicode characters as long I

Hola! I am encountering the following problem, in a multilevel analysis, using glmmPQL from MASS. This occurs with bothj rw1062 and r-devel, respectively with nlme versions 3.1-38 and 3.1-39 (windows XP). > S817.mod1 <- glmmPQL( S817 ~ MIEMBROScat+S901+S902A+S923+URBRUR+REGION+ + S102+S103+S106A+S108+S110A+S109A+S202+S401+S557A+S557B+ + YHOGFcat,

I?m in the process of preparing a package for CRAN. The package is called ?exportR? and since it really just consists of one function, I found it natural to call it ?exportR? too. The function returns a function that does the actual job, but it is set up to work in different ways, depending on the arguments given to its creator. In short: library( exportR ) exporter <- exportR( the, arguments

Hello, I would appreciate if somebody could help me clear my mind about the below issues. I have a factorial experiment to study the effects of Grazing and Fire on Forest biomass production. The experimental unit (to which the treatment combinations are applied) are PLOTs. The measures were made repeatedly for 13 years. I am planning to use the linear mixed effect model function lme in R for this.

Dear all, W.XP I recently upgraded to R2.10.0 and did what I usually do: copied all "special" libraries from the old installation and then used "update.packages()" on the command line. Nothing happened. Then I clicked on the same command on the droplist which generated "update.packages(ask='graphics')", but still with no effect. When I tried to load the

Hello all, I posted a similar question recently but I suspect that it was not well enough formulated to trigger any answers. So I try again: Is there a way to produce boxplots (or something similar) that uses the estimated fixed effects of an lme{nlme} object? When I want to know the mean, I use estimable{gmodels} together with an appropriate matrix. I did look for a similar tool to estimable in

Hi, I have a question about special character: when I create a data frame including some special characters, like '?', it displayed as '<U+00F8>'. I understand that it's one encoding code for '?'. but I want to display the letter as '?' on my screen. And more, when I save the data frame to a local position, it was also save as the encoding code

Hello all! R2.0.1, W2k. All packages updated. I?m heavily dependant on using mixed models. Up til?now I have used lme() from nlme as I have been told to. Together with estimable() from gmodels it works smooth. I also often run Rcmdr, mostly for quick graphics. After using Rcmdr, on reopening the R workspace all help libraries for Rcmdr (22 !) loads, among them nlme, but not Rcmdr itself. Why?

Dear all, W2k, R 2.5.1 I am trying out the AMMI function in the agricolae package, to analyse the dependence of environment for a certain cultivar. The function responds to four basic variables: ENV Environment GEN Genotype REP Replication Y Response My question concerns the 'REP' term. When I normally do an analysis of variance, the replication would refer to repeated observations

Dear everybody, I'm a new user of R 2.4.1 and I'm searching for information on improving the output of regression tree graphs. In the terminal nodes I am up to now able to indicate the number of values (n) and the mean of all values in this terminal node by the command > text(tree, use.n=T, xpd=T) Yet I would like to indicate automatically in the output graph of the tree some

Hello all, R2.1.1, W2k I try to make a plot of a simple regression model in this way: > with(njfA_bcd, { + plot(TC_OS.G31,Prot,cex = 2, col = "red", xlab= "TC/OS at GS32", + ylab="Grain crude protein (CP)") + }) This part works well and produces the datapoints as red circles. When I try to add a line, using a fitted linear model in a way that works perfect

Hello all! System: R2.0.1, W2k. All packages apdated with update.packages(). I use the default graphic device for my system (no active choice). I?m trying to fine-tune scatterplot() graphs (package car), for publication. I want to control plotting symbols and colours to match plots from other, less flexible, systems. I also need to make all text and symbols bigger, to enable printing of the

Hi R-colleagues. In the sem-package i have a problem to introduce hidden variables. As a simple example I take an ordinary factor analysis. The program: cmat=c(0.14855886, 0.05774635, 0.08003300, 0.04900990, 0.05774635, 0.18042029, 0.11213013, 0.03752475, 0.08003300, 0.11213013, 0.24646337, 0.03609901, 0.04900990, 0.03752475, 0.03609901, 0.31702970)

Hello, I decided to try to compile R version 1.0.0 on my home-hacked linux system. This is not a conventional distribution, but has been upgraded countless numbers of times (I think it started as a 1.2.13). It's up to 2.2.13 (kernel version) and I still use gcc 2.7.2.3. I have managed to deal with the fortran part of R with f2c. As it turns out I can compile the program OK, but when it

Dear all! W2k, R 2.5.1 I am working with an ongoing malting barley variety evaluation within Sweden. The structure is 25 cultivars tested each year at four sites, in field trials with three replicates and 'lattice' structure (the replicates are divided into five sub blocks in a structured way). As we are normally keeping around 15 varieties from each year to the next, and take in 10 new

Hi all. After upgrading to R2.0.1 i get > "??" [1] "??" > "??" [1] "\370" > "??" [1] "??" Whereas under R1.9.1 i get > "??" [1] "??" > "??" [1] "??" > "??" [1] "??" Any hints apreciated. Steen Steen Ladelund, statistician +4543233275

Dear all, A question related to the following has been asked on R-help before, but I could not find any answer to it. Input will be much appreciated. I got an unexpected sign of the "slope" parameter associated with a covariate (diam) using zeroinfl(). It led me to compare the estimates given by zeroinfl() and hurdle(): The (significant) negative estimate here is surprising, given

Greetings, I am trying to install R (0.63.2) on a Digital Unix 4.0b Alpha Server 2100A using gcc 2.8.1 and f77 v 0.5.2.3. It seems to compile OK. However, when I try to run R I get the following message: R : Copyright 1999, The R Development Core Team Version 0.63.2 (January 12, 1999) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under