Search Archives

Dear all, Sometimes, during an R session, my computer hangs and I loose all the objects created during this session. Is there a way to automatically type save.image() and savehistory() every 5 minutes? Best regards, Remigijus

Dear All, A very simple question: > library(survival) > fit <- coxph(Surv(time, status) ~ x, data=aml) > survfit(fit) Call: survfit.coxph(object = fit) n events median 0.95LCL 0.95UCL 23 18 30 18 45 I believe I know what...

Dear all, I have two vectors and connect the points by line segments: x <- c(1990,1994,1995,1997) y <- c(30,40,80,20) plot(x,y,type="l") I want to get the values of y's on...

Hello, Could someone help me to explain the VERY big difference in applying two tests on multivariate normality: library(mvnormtest) data(EuStockMarkets) mshapiro.test(t(EuStockMarkets[15:29,1:4])) Shapiro-Wilk normality test data: Z W = 0.8161, p...

Hello R-help, My vector aa has some NA's: > aa_c(NA,1) > aa=="NA" [1] NA FALSE > aa==NA [1] NA NA How can I detect NA's? In other words, I would like my answer to be...

Hello, In a toy example below, there are some outliers which I would like to mark with their rownames. However, it seems that identify(...) does not work together with boxplot? set.seed(1) x<-data.frame(rcauchy(20)) rownames(x...

Dear All, I am on WindowsXP with 512 MB of RAM, R 2.4.0, and I want to read in a big file mln100.txt. The file is 390MB big, it contains a column of 100 millions integers. > mln100...

Dear List, There has just been this question (as well as the answer to it) in the S list: > I would like to put some text in plots... but in french. How can I > write accented letters? > I suppose that...

Hello, I want to test if the Weibull distribution is appropriate for the failure time. When trying to reproduce an example from MASS (the book, Ch. 13.2), I type library(survival) library(MASS) leuk.wei <- survreg( Surv(time)~ag...

Dear All, R 1.9.1, Windows When copying and pasting a few lines from the 'predict.Arima' help, I get an error message: > data(lh) > predict(arima(lh, order = c(3,0,0)), n.ahead = 12) Error in eval...

Dear All, There is an example for the priceIts function (the its package) which does not work for me as expected. > ?priceIts > x1 <- priceIts(instrument = c("^ftse"), start = "1998-01-01", + quote = "Close") Error in validObject(.Object) : invalid class "its...

Dear All, I have the test results in the following form > test <- data.frame(Name=letters[1:5],Score=c(7,NA,5,NA,10)) > test Name Score 1 a 7 2 b NA 3 c 5 4 d NA...

If you define Pareto density as p(x)=c*x^(-(c+1)) for x>1, then dpareto <- function(x,c){ if(c<=0)stop("c must be positive") # Diagnostic step ifelse(x<1,0,c/x^(c+1))} ppareto <- function...

Many thanks to all who replied to my e-mail. My problem was that I had not known about the approx function. By the way, if I have x <- c(1990,1994,1995,1997), is there an automated way to...

Try the H. Bengtsson's function plot.histogram from http://www.maths.lth.se/matstat/staff/hb/mypackages/R/plot.histogram.R Remigijus Saturday, January 25, 2003, 10:21:30 PM, you wrote: FHFPdR> Hi, I am having some trouble...

Dear all, I have a function MYFUN which depends on 3 positive parameters TETA[1], TETA[2], and TETA[3]; x belongs to [0,1]. I integrate the function over [0,0.1], [0.1,0.2] and [0.2...

Sunday, March 02, 2003, 9:48:26 PM, you wrote: MK> ripley at stats.ox.ac.uk wrote: >> You are not solving the same problem, which is more than a little `odd'. MK> Why? MK> >>> want to choose the three...

Dear all, I met with some problems when dealing with a time series with serial correlation. FIRST, I generate a series with correlated errors set.seed(1) x=1:50 y=x+arima.sim(n = 50, list(ar = c(0...

bindec (at least, on my computer) works only for moderate values of b: bindec <- function(b){ as.i <- as.integer(unlist(strsplit(b,""))) print(as.i) fl <- 2^(floor(round(log10(b),10)):0) print(fl) dec <- sum(as.i...

I use the following two function for a lagged regression: lm.lag=function(y,lag=1) summary(lm(embed(y,lag+1)[,1]~embed(y,lag+1)[,2:(lag+1)])) lm.lag.x=function(y,x,lag=1) summary(lm...