similar to: Analysis of Dispersion

I have an unexpected result with the functions lag() and diff() in the plm (panel data) package when used with transform(). These plm-specific functions are supposed to generate lags and first differences within each panel. lag() does not work properly the first time (it reproduces the same series--this is a common time series pitfall), BUT then it does work properly when it is run a second

I am trying to make R-1.8.1 on (SunOS shell1 5.8 Generic_108528-15 sun4u sparc SUNW,UltraAX-i2). I did ./configure make Configure output seems ok. The make proceeds until the following line appears, repeated indefinitely (until I break): ./config.status: ./confstat28489-19881/subs.frag: cannot overwrite existing file I suspect that this may involve write permissions (and maybe the umask set in

I have used earlier rsync on cygwin successfully for a long. I have now moved to cygwin 1.7.0, which required moving to rsync 3.0.6-1 at the same time. (The -1 is for minor cygwin-specific patches.) When I do almost any command with rsync, such as the one below, it breaks at a dup2 after a fork in pipe.c (routine piped_child). Windows gives error code 6, which is translated by cygwin to errno 9 =

This is to both the rsync and cygwin lists, and I regret that it is vaguer than I would like, but ... With cygwin 1.7 (latest beta) and the rsync released with it (3.0.6-1) I try to sync a couple of Gb of Thunderbird mail files. The flag settings were -avzuP plus a --chmod and some includes/excludes to control which files to move. I would repeatably get "hangs" on particular files,

Dear R-help, I have a four-part question about regression, matrices, and sandwich package. 1) In the sandwich package, I would like to better understand the meat() function. >From the bread() documentation, for a simple OLS regression, bread() returns (1/n * X'X)^(-1) That is, for a simple regression (per the documentation on bread()): MyLM <- lm(y ~ x) bread(MyLM)

Dear all, This seems to be working, but I'd like to make sure that I'm not doing anything wrong. I am using by() to construct a complicated summary statistic by several factors in my data (specifically, the 90-50 income ratio by city and race). cityrace.by <- by(microdata, list(microdata$city,microdata$race), function (x) quantile(x$income, probs=0.9) / quantile(x$income, probs=0.5)

R Version 2.12.1 (2010-10-15) vs 2.12.0 has slowed down 8 fold for dual core and 17 fold for dual-core-dual-processor Macs. I have checked this result on 3 different macs using the following R-script: Using Version 2.12.0 on a dual core dual processor Mac: > source("http://www.bio.umass.edu/biology/kunkel/pub/R/CuriousResult.R") matrix multiplication 43.543 1.308 14.788

Hi, We have software raid level 5 using 6 IDE disks amounting to about 500Gb. Its an ext3 filesystem which is NFS exported to a number of other machines. Recently, after running /etc/init.d/nfs-kernel-server stop /etc/init.d/nfs-common stop /etc/init.d/nfs-common start /etc/init.d/nfs-kernel-server start the whole raid array is gone. The file system is consistent and is the right size but no

I cannot get the "R" button to appear in RWinEdt. I'm running R 2.0.1 under Windows XP. I did a clean install of the latest WinEdt. Previously (for an earlier installation of WinEdt) I had RWinEdt running OK. Now, even though I have the distributed RWinEdt_1.6-2.zip extracted to the right place in the R directory, and when I execute library(RWinEdt) in the R Console,

I have a between-subjects factor, A, with three levels, and a within-subjects factor with 4 levels. There are 6 subjects at each A level. I get the correct F tests using the aov function with the stack form of the data frame but I want the epsilon adjustment (Greenhouse-Geisser, Huynh-Feldt). for nonspheriicty as well. I got the SPSS results using the standard broad data frame (a row for

Hi, I just installed version 2.11.1 on my Mac OS X 10.5.2. When I try running R from my xterm, I get the following error: Error in loadNamespace(name) : there is no package called 'lme4' Fatal error: unable to restore saved data in .RData I am returned to the terminal prompt and can not run R. I can open this version of R from the R console and tried installing the lme4 source

I'm a rank beginner with R, and I'm stumped at how to get it to work with Xemacs under Windows XP. I'm afraid there are some VERY basic questions I have here. Yes, I did read http://ess.r-project.org/Manual/readme.html. Both R and Xemacs are installed and working. The current ess files are installed in a subdirectory of my Xemacs directory. I added the R bin directory to my

Hello, I administer a student (Windows based) computer lab, which has R installed. I want to permit students to install R packages, but of course they don't have write permission to R root folder. Following FAQ suggestion, I've set up a folder that they can write to, and set the R_LIB parameter on the shortcut to this folder. So far, so good - they can now install packages.

I just installed R-2.3.1pat under Windows XP as well as the associated RWinEdt. If I start RGui from its shortcut (but do _not_ also start RWinEdt) and then try to execute demo() or help(), I get a RWinEdt pop-up error window with message: File "D:\WP\WinEdtData\WinEdt\D:/WP/WinEdtData/WinEdt/R.ini" does not exist! Qualifier -e/-E does not specify an existing file! First

Greetings, I'd be grateful if a good Samaritan helps me to approach this problem.... with my data, I've created the following model lm(formula = OUTCOME ~ VAR1 + VAR2) summary(model) Call: lm(formula = OUTCOME ~ VAR1 + VAR2) Residuals: Min 1Q Median 3Q Max -1.4341 -0.3621 0.1879 0.4994 0.7696 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.89020

stuck on this in spite of finding a few items in the archive that seemed pertinent. just configuring a new installation of 2.2.8 from SunFreeware on a Solaris 8 system. installation done, swat running, did simple basic configuration. then went to terminal window and got /usr/local/samba/bin# ./smbpasswd -j provost -r umapvstads1 -U libadmin Password: Error connecting to

Hello, I'm trying to plot a fitted lm() line on a plot when the one explanatory variable is log transformed and log="x". I get different lines using abline and predict.lm(). #Example x <- 1:100 y <- rnorm(100) plot(y ~ x, log="x") abline(lm(y ~ log(x))) lines(x, predict(lm(y ~ log(x))), lwd=2) I'm sure I'm missing something but could someone tell me which

Hello. I'm getting an unexpected result when running smooth.spline(). Here is a simple example that replicates the error I'm getting: > aa <- c(1, 2, 3, 8, 8, 8, 8, 8, 8, 8, 8, 8, 12, 13, 14) > bb <- 1:length(aa) > plot(aa, bb) > smooth.spline(aa, bb) Error in smooth.spline(aa, bb) : need at least four unique 'x' values As you can see from the example, my

#This is probably due to my incomplete understanding of grep, but the below code has been working for some time to #search for .R with anything in front of it and return a list of scripts to source. Likely, the syntax for the #grep statement has been wrong all along. scripts2source <- (c("/home/ssefick/R_scripts/Convert_package.R",

Hello, This is an (hopefully) improved question of one I posted several weeks ago. Does anyone know of a function for fitting "two-part" models? These models are designed to handle count data with so many zeroes that they can't be fit well with zero-inflated Poisson models or other 'typical' GLMs. My understanding is that they work by first fitting a binomial model to