similar to: alarm() doesn't beep

Dear all I'm using Xubuntu Lucid and I keep getting the same random numbers whenever I start a new session of R. For example, I keep getting > sample(1:1000, 1) [1] 87 or > rnorm(1:10) [1] -1.3618103 0.4241701 1.0720076 0.2208145 -0.5375314 -0.4846588 [7] 0.7576768 0.6527407 -0.6868786 0.8718527 I expected that some set.seed() instruction woudl be present in a config file in

I'm creating a stacked bar chart using the likert command in the HH package. My data are in a data frame, with two numeric variables and a categorical variable, I can't get likert to use the column containing the categorical variable as a my y axis label. Here is a quick example: library(HH) #my data are:

Dear all I am working on unbalanced panel data and I can readily fit a "pooling" model using plm(), but not a "within" or "random" model. Reproducing the examples in vignette("plm") and in the AER package I encountered no such issues. ##unfortunately I cannot disclose the data, and it is too big anyway > dim(ibes.kld.exp.p[x.subs , ]) [1] 13189 34

Dear all This may be obvious, but I cannot get it working. I'm trying to subset & sort a data frame in one go. x <- iris x$Species1 <- as.character(x$Species) ##subsetting alone works fine with(x, x[Sepal.Length==6.7,]) ##sorting alone works fine with(x, x[order(Sepal.Length, rev(sort(Species1))),]) ##gets subsetted, but not sorted as expected with(x, x[(Sepal.Length==6.7) &

Dear R users I would like to do some spreadsheet style expansion of dates. For example, I would need to obtain a vector of months. I approached in an obviously wrong way: > paste(01:12) [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" > as.Date(paste(01:12),

Dear all, I installed the amd64 cran2deb [1] JGR binaries on a fresh Debian testing, and JGR is missbehaving. When JGR starts up, one core of the CPU goes to 100% and stays there, even though I am keeping it idle and not performing any operation. I tried JGR with both sun-java6-jdk and openjdk-6-jdk, debian-liv:/home/liviu# update-alternatives --config java There are 4 choices for the alternative

Dear all, Given the following vector: > (z <- c('R project', 'hello world', 'something Else')) [1] "R project" "hello world" "something Else" I know how to obtain all capitals or all lower case letters: > tolower(z) [1] "r project" "hello world" "something else" > toupper(z) [1] "R

Dear all Is there any R function that would perform currency conversion using up-to-date exchange rates? I would be looking for a function that allows to download recent exchange rates (say, from Yahoo!) and then use these in converting currencies (say, USD to EUR). I am not sure whether r-sig-finance would be more appropriate, but the (off-)topic feels general enough to me. Thank you Liviu --

Dear all How does one convert a "non-symmetric" list to a vector? See below: > x <- list() > x[[1]] <- letters[1:5] > x[[2]] <- letters[6:10] > x[[3]] <- letters[11:12] > x [[1]] [1] "a" "b" "c" "d" "e" [[2]] [1] "f" "g" "h" "i" "j" [[3]] [1] "k"

Dear all I'm unable to find an example of extracting the rotated scores of a principal components analysis. I can do this easily for the un-rotated version. data(mtcars) .PC <- princomp(~am+carb+cyl+disp+drat+gear+hp+mpg, cor=TRUE, data=mtcars) unclass(loadings(.PC)) # component loadings summary(.PC) # proportions of variance mtcars$PC1 <- .PC$scores[,1] # extract un-rotated scores of

Dear all I am getting this strange error when checking my package. Would you have an idea what causes it? Thank you Liviu * checking Rd cross-references ... WARNING Error in .find.package(package, lib.loc) : there is no package called 'KernSmooth' Calls: <Anonymous> -> lapply -> FUN -> .find.package Execution halted > sessionInfo () R version 2.10.0 (2009-10-26)

Dear all, I'm trying to gsub() "%" with "\%" with no obvious success. > temp1 <- c("mean", "sd", "0%", "25%", "50%", "75%", "100%") > temp1 [1] "mean" "sd" "0%" "25%" "50%" "75%" "100%" > gsub("%",

Dear all I would like to start R with Rcmdr from the cli, without tweaking Rprofile.site. This has been discussed in the past [1], but I don't see a solution that (1) could be used with any working directory and (2) would avoid starting Rcmdr on every R start-up. Personally I tried the following, which starts R but not Rcmdr liv at liv-laptop:~$ R --interactive -e 'require(Rcmdr)'

Dear all How does print.htest display the p-value in scientific notation? > (x <- cor.test(iris[[1]], iris[[3]])) Pearson's product-moment correlation data: iris[[1]] and iris[[3]] t = 21.65, df = 148, p-value < 2.2e-16 alternative hypothesis: true correlation is not equal to 0 95 percent confidence interval: 0.8270 0.9055 sample estimates: cor 0.8718 Above the p-value comes

Dear all I have a 'character' vector containing mixed formats (thanks Excel!) and I'd like to translate it into a default "%Y-%m-%d" Date vector. x <- c("1/3/2005", "13/04/2004", "2/5/2005", "2/5/2005", "7/5/2007", "22/04/2004", "21/04/2005", "20080430", "13/05/2003",

Dear all, Does this make any sense: var() = cov() != acf(lag.max=0, type="covariance")? I have daily data of IBM for May 2005, and I'm using the logarithmic return: > ibm200505$LRAdj.Close [1] NA 0.0203152 0.0005508 -0.0148397 -0.0025182 0.0092025 -0.0013889 [8] 0.0098196 -0.0103757 -0.0274917 0.0005716 -0.0159842 -0.0074306 0.0091710 [15] 0.0002898 0.0226306

Dear all Is there a clean plyr version of the following by() and do.call(rbind, ...) construct: > df<-data.frame(a=1:10,b=11:20,c=21:30,grp1=c("x","y"),grp2=c("x","y"),grp3=c("x","y")) > dfsum<-by(df[c("a","b","c")], df[c("grp1","grp2","grp3")], range) >

Dear all Is it possible to set globally the option .progress = "text" to all the apply functions in 'plyr'. For example, current default is daply(..., .progress = "none"). I would like to set it to daply(..., .progress = "text"), so as to avoid writing the argument every time I call such a function. I looked into ?daply and ?create_progress_bar without much

Dear all It is difficult to use round(..., digits=2) on a data frame since one has to first take care to remove non-numeric variables such as 'character' or 'factor': > head(round(iris, 2)) Error in Math.data.frame(list(Sepal.Length = c(5.1, 4.9, 4.7, 4.6, 5, : non-numeric variable in data frame: Species > head(round(iris[1:4], 2)) Sepal.Length Sepal.Width Petal.Length

Dear all What is the easy way to drop a variable by using its name (and not its number)? Example: > data(iris) > head(iris) Sepal.Length Sepal.Width Petal.Length Petal.Width Species 1 5.1 3.5 1.4 0.2 setosa 2 4.9 3.0 1.4 0.2 setosa 3 4.7 3.2 1.3 0.2 setosa 4 4.6 3.1