similar to: incl.non.slopes=FALSE does not work at predict.lm

Dear all, I want to use the "rearrange" command which is based on Chernozhukov et al paper and is included in the quantreg package. So, I run a quantile regression in which I included dummy variables for state and years in order to estimate the respective fixed effects quantile regression. The problems are the followings: 1. At example that is stated in the help****, I don't

Hi, I am trying to use (i) as an index but R considers it as a function and not as text. To be more specific I would like for example to estimate some regressions named qrnox1, qrnox2, qrnox3,..... and so on. But when I am using qrnox(i) ot qrnox[i] it tries to find the ith element of vector qrnox. The thing is that I want to estimate the qrnoxi regression and not the qrnox(i) function or

Dear all, I am running a code using bootstraps for estimating standard errors but the mac crashes. When I use small number of bootstraps (100) it works fine but if I increase that number it crashes. Thanks in advance dimitris The code, the error and my mac characteristics are the following: ##############code##################### qr.1<- rq(y~factor(year)+factor(state)+x1+I(x^2)+I(x^3),

Hi, I cannot find a 'vectorized' solution to this 'for loop' kind of problem. Do you see a vectorized, fast-running solution? Objective: Take the value of X at each timepoint and calculate the corresponding value of Y. Leading 0's and all 1's for X are assigned to Y; otherwise Y is incremented by the number of 0's adjacent to the last 1. The frequency and

Dear R users I want to compute Kendall's Tau between two vectors x and y. But x and y may have zeros in the same position(s) and I wrote the following function to be sure to drop out those "double zeros" "cor.kendall" <- function(x,y) { nox <- c() noy <- c() # for (i in 1:length(x)) if (x[i]!= 0 | y[i] != 0) nox[length(nox)+1]<- x[i] for (i in

Greetings, I am working on a platform for analysis automation. I need to run several Virtual Environments concurrently and record information about their behavior. I wrote some months ago about the capability of reading the Memory during the Environment's execution (in paused state). What do I need is the complete linear memory image, byte per byte, nothing special; I will give this output

Hello I am trying to use all subsets regression on a test dataset consisting of 11 trails and 46 potential predictor variables. I would like to use Mallow's Cp as a selection criterion. The leaps function would provide the required output but does not work with this many variables (see below). The alternative function regsubsets should be used, but I am not able to define the function in

Dear list I'm plotting ( boxplot() and plot() ) some data for a publication. The editor would like the text labels on the plots in a larger font. I'm doing something like this: <snip> jpeg( filename = "D:/Martin/Work/CleanPath/RAF1%03d.jpg", width = 1000, height = 600, pointsize = 12, quality = 100, bg = "white", res = 96, restoreConsole = TRUE )

Hi all, I'm trying to find a way to get xyplot to produce a second set of axes outside the right hand side of the graph. This is my progress so far: EE <- equal.count(ethanol$E, number=9, overlap=1/4) xyplot(NOx ~ C | EE, data = ethanol, prepanel = function(x, y) prepanel.loess(x, y, span = 1), xlab = "Compression Ratio", ylab = "NOx (micrograms/J)",

So we got an Ubuntu 10.10 VPS, and go there with SSH putty, working with Xming as X server. This is what we did so far, to host a dedicated nox server. Code: apt-get update sudo apt-get install wine sudo apt-get install xorg wget *noxinstall url* unzip noxinstall.zip wine NoxInstall.exe (an installer and we put it in /root/Nox/) Wine wouldn't recognice a CD, so a crack gets installed

Dear all I am trying to fit Variable Coefficients Models on Unbalanced Panel Data. I managed to fit such models on balanced panel data (the example from the "plm" vignette), but I failed to do so on my real, unbalanced panel data. I can reproduce the error on a modified example from the vignette: > require(plm) > data("Hedonic") > Hed <- pvcm(mv ~ crim + zn + indus

I ran the following script from xyplot Examples using Tin-R on Windows and saw no plot produced. EE <- equal.count(ethanol$E, number=9, overlap=1/4) xyplot(NOx ~ C | EE, data=ethanol, prepanel = function(x,y) prepanel.loess(x, y, span=1), xlab="Compression Ratio", ylab="NOx (micrograms/J)", panel = function(x,y) { panel.grid()(h = -1, v=2)

I need to split a data frame into 3 columns. The column I want to split contains indices of lag (prefix L1 or L2 and suffix 01, 03, 04), station name (shown in the sample data as capitalized G, P and S) and pollutant name. Names with no ?L? prefix or 01/04 suffix are lag 0. Lag 01 is average of lag 0 and 1, and 04 is average of 0 to 4 days. How can one do that in R? I will ignore the other

Dear all, I'm having some problems with using the step() function inside another function. I think it is an environment problem but I do not know how to overcome it. Any suggestions are appreciated. I've prepared a simple example to illustrate my problem: > library(MASS) > data(Boston) > my.fun <- function(dataset) { + l <- lm(medv ~ .,data=dataset) + final.l <-

I have found a few issues with puppet and the up2date provider (version 0.22.4) that I''d like to run by those who use up2date (to see if others are experiencing the same issues). 1. Puppet doesn''t seem to find the latest updates of packages that are already installed. In looking at "/usr/lib/site_ruby/1.8/puppet/provider/package/up2date.rb", I found that

Hello, I am using the following iptables POSTROUTING rule to NAT some RFC 1918 addresses: iptables -t nat -A POSTROUTING -s 192.168.19.23 ! 192.168.0.0/255.255.0.0 -p tcp --dport 80 -j SNAT --to-source 10.32.4.2 (I am using SNAT instead of MASQUERADE for performance reasons). I have several addresses on the 192.168.0.0/16 subnet that I am SNAT''ing similarly. Problem is, ''tc

Hello, Any word on when extended attribute support will be added to OCFS2? What are the impediments to someone implementing this? Alternatively, has anyone implemented a patch to add extended attribute support to OCFS2? I've been watching the open source Lustre fs development, too. I noticed that they have added extended attribute support. However, OCFS2 is more attractive for some

Hi: (Tried to find a bug report about this issue, but was unable to find it, let me know if this is a known issue) I have been working on an interface to highlight xyplot panels on mouse overs in JavaGD but I have stumbled with what seems to be a bug in trellis.focus. I am using R 2.8 with lattice 0.17-15 *** To replicate the bug: 1.- display an xyplot. For example, from the xyplot help page:

Hello, When I kick off SVN1992 with storeconfigs=yes, I get the following error: err: Could not store configs: SQLite3::SQLException: SQL logic error or missing database: INSERT INTO resources ("exported", "line", "title", "host_id", "restype", "source_file_id") VALUES(NULL, 23, ''svnpkgadd-emacs-nox-etc'', NULL, NULL,

Hi all! I want to thanks all people here for their help, good job guys! :o) And nox, it's my question: I'm using smbldap-tools 0.8.2 from samba.idealx.org. In all the docs I read about it, I read that I must put these lines in smb.conf: passwd program = /usr/local/sbin/smbldap-passwd.pl -o %u passwd chat = *new*password* %n\n *new*password* %n\n *successfully* I look in the code of