similar to: change factor labels

hi, this may be trivial, but we can't seem to find anything adequate, (although there is a work around with match() ). We are looking for something along the lines of plot(table1[table1$var2==or("a","b","c","d"),"var1"]) would be handy, with the potential or() function leading to what plot(table1[table1$var2=="a" |

Dear list, My apologies if a solution / explanation to this already exists on the list, but it is difficult to assign it to a certain keyword. test<-c(1:3) testfct <- function(x) {test[1]<-100} test [1] 1 2 3 testfct(1) [1] 1 2 3 Basically, I would like to write data into an external table that the function does not know. Why is this not working / what alternatives exist? Thanks,

Hello A very short question: Using multiple plots with par(mfrow=c(3,3)), how can I get R to indicate one common y- and x-axis label? I tried to use text() in par, but this is then overwritten when I plot the graphs. -- Sebastian Leuzinger Institute of Botany, University of Basel Schönbeinstr. 6 CH-4056 Basel Ph. 0041 (0) 61 267 3511 fax 0041 (0) 61 2673504 email: Sebastian.Leuzinger@unibas.ch

Dear R-list I am (forced) to change from Linux to Mac and am now looking for a new editor for R. I would like one that features a split window (console + editor) as well as syntax highlighting. Can anyone help? Especially the split-window feature does not seem to be easily available in the editors desribed on the R-help site, except Emacs, which I am reluctant to start using. I am on a MacBook Air

Hello, anyone got an idea on how to use stl() so that the remainder eventually becomes white noise? i used stl repeatedly but there is autocorrelation in the remainder that i can't get rid of. os: linux suse9.3 ------------------------------------------------ Sebastian Leuzinger Institute of Botany, University of Basel Sch??nbeinstr. 6 CH-4056 Basel ph 0041 (0) 61 2673511 fax 0041 (0)

hello, has anyone come across the following rather mysterious problem: when pasting large bits of code (100 and more lines) into the R console with the central mouse button (under linux), only part of the code is pasted, and the text interrupts somewhere arbitrarily. It does not happen when smaller bits are pasted subsequently. I use linux suse 9.3 with the latest version of R

Hello Trying to install the MySQL package, I get the following error. The help archive contains something on this issue but did not help. I work on linux suse 9.3 Configuration error: Could not locate the library "libz" required by MySQL. The library libz however is not on any mirrors I checked. ------------------------------------------------ Sebastian Leuzinger web

hello, is R able to suppress tick labels (not tick marks)? i know there is a way around this with axes=F and then draw new axes, but it would be easier to suppress them in the first place. -- ------------------------------------------------ Sebastian Leuzinger Institute of Botany, University of Basel Sch??nbeinstr. 6 CH-4056 Basel ph 0041 (0) 61 2673511 fax 0041 (0) 61 2673504 email

dear list while installing the package rimage, i get the error message checking jpeglib.h usability... no checking jpeglib.h presence... no checking for jpeglib.h... no configure: error: Sorry, can't find jpeglib header ERROR: configuration failed for package 'rimage' although i have installed jpeglib as required. does R not find this library? but why did it find the fftw library

Hi I run suse linux 9.1 and I installed MySQL server, client, devel, bench. DBI is installed, when I try to install RMySQL I get an error saying, that libz is missing. (paths to libs were set:export PKG_CPPFLAGS="-I/usr/include/mysql/" export PKG_LIBS="-L/usr/lib/mysql/ -lmysqlclient") so my question: where do I get the libz files (are these mysql files? if yes, why were

Hi everyone I'm using Tukey HSD as post-hoc test following a lme analysis. I'm measuring hemicelluloses in different species treated with three different CO2 concentrations (l=low, m=medium, h=high). The whole experiment is a split-plot design and the Tukey-function from the package multcomp is suitable for lme-analysis with random factors. The analysis works fine but I get a non

Hi Folks, I think I have searched exhaustively, including, of course R-help (D. Bates, S. Graves, and others) and but I remain uncertain about testing fixed effects with lmer(..., family=binomial). I gather that mcmcsamp does not work with Do we rely exclusively on z values of model parameters, or could we use anova() with likelihood ratios, AIC and BIC, with (or without)

Hi everyone I have the following problem: My dataframe has 3 variables: ID, Year and and an outcome variable. The dataframe contains repeated measurements because the subjects filled out a questionnaire every year. The time span covers 2 years. Now I want to check if there is a significant change in the outcome over the 2 years with a paired wilcox.test. The problem: Not every subject

Dear all I imported a Stata .dta file with the read.dta-function from the foreign-package. The dataframe's dimensions are > dim(d.apc) [1] 15806 1300 Importing needs up to 15 min and calculations with these data are rather slow (although I subset the data before starting analyses). My questions are: 1. Has someone experiences importing Stata files (alternatives to read.dta) ? 2.

Dear R-list member, I'm using the function pmnorm() (-->library(mnormt)) within a for-loop. Certain parameter values leads to an error message: "(In sqrt(diag(S)) : NaNs produced, In sqrt(1/diag(V)) : NaNs produced, In cov2cor(S) : diag(.) had 0 or NA entries; non-finite result is doubtful)" obviously because "NaNs" were produced. Is it possible to tell R that it

Dear r-list, I am trying to visually seperate the two main clusters of a dendrogram. The idea is to use: 'edgePar=list(lty=3)' for 'dend1[[1]]' and 'edgePar=list(lty=1)' for 'dend1[[2]]' I have not found a way to solve this. Any suggestions? Patrick hc <- hclust(dist(USArrests), "ave") (dend1 <- as.dendrogram(hc)) par(mfrow=c(2,2)) plot(dend1)

Hi all Maybe a stupid question, but I'm not able to figure this out from the manpage nor from the HOWTOs... How can I force a user to change his password at next logon? I tried: # pdbedit -P "user must logon to change password" -C 0 -u username # pdbedit -u username -P "user must logon to change password" -C 0 # pdbedit -P "user must logon to change password"

Hello. Slope estimates in lme are shrinkage estimates which pull the OLS slope estimates towards the population estimates, the degree of which depends on the group sample size and the distance between the group-based estimate and the overall population estimate. Although these shrinkage estimates as said to be more precise with respect to the true values, they are also biased. So there is a

Dear list, I want to retrieve the p-value of a two-polynomial regression. For a one-polynomial lm I can easily do this with: summary(lm(b~a, data=c)[[4]][[8]]. But how do I find the final p-value in the two-polynomial regression? Under $coefficients I don't find it Any suggestions? Patrick alt <-(2260,2183,2189,1930,2435, 2000,2100,2050,2020,2470, 1700,2310,2090,1560,2060,

Hi list I just compiled 3.0rc1 on Solaris 9 (UltraSparc IIe) with gcc 3.3. Seems to work fine but I am surprised by the enormous size of the binaries: (ub-blade1):/usr/local/samba/sbin # ll total 232272 -rwxr-xr-x 1 root other 44719632 Aug 28 15:07 nmbd -rwxr-xr-x 1 root other 89726544 Aug 28 15:07 smbd -rwxr-xr-x 1 root other 44807356 Aug 28 15:07 swat -rwxr-xr-x