similar to: JSS Volume 18 (www.jstatsoft.org)

At www.jstatsoft.org we now have three special volumes Volume 22, Ecology and Ecological Modelling in R (Thomas Kneib and Thomas Petzoldt, eds.) Volume 20, Psychometrics in R (Jan de Leeuw and Patrick Mair, eds.) Volume 18, Spectroscopy and Chemometrics in R (Kate Mullen and Ivo van Stokkum, eds.) Coming soon -- Econometrics in R -- Political Methodology in R Gleam in my eye --

Hi I have some panel data for the 50 US states over about 25 years, and I would like to test a simple model via OLS, using this data. I know how to run OLS in R, and I think I can see how to create Panel Corrected Standard Errors using http://jackman.stanford.edu/classes/350C/pcse.r What I can't figure out is how to correct for autocorrelation over time. I have found a lot of R stuff on

in 2007 we will publish seven volumes of JSS (four of them special volumes). Volume 18 Spectroscopy and Chemometrics in R Volume 19 regular Volume 20 Psychometrics in R Volume 21 regular Volume 22 Ecology and Ecological Modeling in R Volume 23 is being started, it is a regular volume, and may not fill up completely. Volume 24 is 'Statistical Modeling of Social Networks with

Dear All, I am trying to find MLE by using "OPTIM" function. Difficult in differentiating some parameter in my objective function, I would like to use the returned hessian matrix to yield an estimate of Fisher's Information matrix. My question: Since the hessian is calculated by numerical differentiate, is it a reliable estimate? Otherwise I would have to do a lot of work to

Hello all, I am wondering if there is a way to do a spatial error probit/logit model in R? I can't seem to find it in any of the packages. I can do it in MATLAB with Gibbs sampling, but would like to confirm the results. Ideally I would like to use this model to predict probability of parcel conversion in a future time period. This seems especially difficult in a binary outcome model

Hello, I wonder why the control parameter REPORT is not supported by method SANN. Looking into optim.c I found an internal constant: #define STEPS 100 ... and decreasing this to 10 helped me fine-tuning the annealing parameters in an actual problem. Is there any reason why not passing nREPORT to samin and setting something like: STEPS = nREPORT / tmax Thomas P. -- Thomas Petzoldt

Dear R developers, I recently observed a NOTE on several MaxOS X package checks: sh: line 1: gs: command not found !!! Error: Closing Ghostscript (exit status: 127)! /usr/bin/texi2dvi: thumbpdf exited with bad status, quitting. See for details: http://www.r-project.org/nosvn/R.check/r-release-macosx-ix86/simecol-00check.html or

Hi All, I am currently weaving samba into an active directory domain. ntlm works fine. but wbinfo is not so good. wbinfo -r username returns Could not get groups for user username wbinfo -n username returns S-1-5-21-1482476501-343818398-682003330-6830 User (1) wbinfo -a username%password plaintext password authentication failed error code was NT_STATUS_NO_SUCH_USER (0xc0000064) error messsage

This could well be a package-specific issue, and so perhaps not for R-devel. The author of the relevant package is CC'd in any event. I'm seeing some odd behavior from Sweave with keep.source=TRUE when using driver=cacheSweaveDriver from the cacheSweave package. It would seem that code chunks with keep.source=TRUE aren't being parsed correctly; the corresponding Schunk starts at

Dear R developers, I tried to use nlminb instead of optim for a current problem (fitting parameters of a differential equation model). The PORT algorithm converged much better than any of optim's methods and the identified parameters are plausible. However, it took me a while before spotting the reason of a technical problem that nlminb, in contrast to optim, does not pass names of the

Dear all, The default plotting method for hclust trees looks just fine for few objects like in the example dataset. But when it comes to many features (eg some 1000 and more - I'm trying to visualize clustered microarray data) it renders a tree, that one cannot inspect, because of overlapping text and lines. My question is, is there a way or a plotting parameter for plotting a tree which is

PowerPoint Viewer 1997 work nicely, but the installation for PowerPoint Viewer 2003 seems incomplete. It emits a number of fixme:msi warnings. It would be a nice wishlist item. Cheers, Shaun $ wine --version Wine 20041201 $ md5sum ppviewer.exe 1facc59a11680d1f828ed5118a0270b1 ppviewer.exe $ wine ppviewer.exe fixme:msi:MsiInstallProductW L"ppviewer.msi" L""

Dear R developers, some of our packages come with additional programming examples in a directory called "/examples" which is created from "/inst/examples". This directory is linked from the docs (e.g. in inst/doc/index.html): <dl> <dt><a href="../examples/">examples</a>: <dd>Source code of examples </dl> Given, that we have a

Dear list - Does anyone have any ideas / comments about why I am receiving the following warning when I run lsoda: 1: lsoda-- at t (=r1), too much accuracy requested in: lsoda(start, times, model, parms) 2: for precision of machine.. see tolsf (=r2) in: lsoda(start, times, model, parms) I have tried changing both rtol and atol but without success. I saw the thread in the

Dear list The deSolve package allows you to specify the model code in C or Fortran. Thanks to the excellent vignette this works fine. However I have not yet managed to use forcing functions in C code. In pure R code this works very well with approxfun() specified outside the model: ############################################### #Model lvml <- function(t, x, parms) {

Dear all, I am attempting to explain patterns of arthropod family richness (count data) using a regression model. It seems to be able to do a pretty good job as an explanatory model (i.e. demonstrating relationships between dependent and independent variables), but it has systematic problems as a predictive model: It is biased high at low observed values of family richness and biased low at

homals-0.9.0 is on CRAN -- by Jan de Leeuw and Patrick Mair This package implements the methods discussed in Gifi, Nonlinear Multivariate Analysis, Wiley, 1990. In the Gifi terminology it covers homals, princals, canals, morals, criminals, and overals. The R implementation fills several gaps in Gifi, adding multiple ordinal, numerical, and polynomial data transformations. Differences with

homals-0.9.0 is on CRAN -- by Jan de Leeuw and Patrick Mair This package implements the methods discussed in Gifi, Nonlinear Multivariate Analysis, Wiley, 1990. In the Gifi terminology it covers homals, princals, canals, morals, criminals, and overals. The R implementation fills several gaps in Gifi, adding multiple ordinal, numerical, and polynomial data transformations. Differences with

http://idisk.mac.com/jdeleeuw-Public/jacobi This is paper/software for various techniques based on Jacobi plane rotations. There is R code for -- classical cyclical Jacobi Eigen diagonalization -- Jacobi-based SVD diagonalization -- approximate simultaneous diagonalization of symmetric matrices (De Leeuw/Pruzansky 1978) -- approximate simultaneous diagonalization of rectangular matrices

Hello One of my packages, untb_1.3-2, passes R CMD check under MacOSX (and apparently the systems used in the package check summary page on CRAN) but fails with the following message on R-2.6.0.tgz compiled last night on my (home) linux box. I hasten to add that I have never seen this error before on home-compiled pre-releases of R-2.6.0. Can anyone help me understand what is going on?