Re: [R] Interpolating / smoothing missing time series data

Francisco J. Zagmutt wrote: I don't have much experience in the subject but it seems that library(akima) should be useful for your problem. Try library(help=akima) to see a list of the functions available in the library. I hope this helps Francisco Yes, function aspline() of package

Re: [R] Creating very small plots (2.5 cm wide) in Sweave

Dear Francisco, thanks for your solution. It turns out that it's best for me to use \setkeys{Gin}{width=0.15\textwidth} directly before I call the plot - that seems to work just fine. Andrwe On Thu, Sep 08, 2005 at 05:44:59AM +, Francisco J. Zagmutt wrote: Others may propose more

[R] Text Size in Legend

Hello, I need to reduce the size of the text in a legend since the legend is overlapping with the curves in my plot. I've not been able to identify any way to achieve this in the documentation. Anyone have any suggestions on how to scale down the text or the overall legend? Thanks in advance

Re: [R] Text Size in Legend

Hi Chris, To change the scale of the whole legend you can use the argument 'cex' in legend. I hope this helps! Regards, Roula = Spyridoula Tsonaka Doctoral Student Biostatistical Centre Catholic University of Leuven Kapucijnenvoer 35 B-3000 Leuven Belgium Tel: +32/16/336899

[R] Using R map data to determine associated state for a coordinate?

Hi! I have no idea if this is maybe an easy task utilizing R since I read there is geographical map data in some package: I have a huge number of geographical points with their coordinates in Germany. Now I want to determine for each point in which Bundesland = state it is located. Can

[R] Too long to display problem

Dear list, I used read.xls in gdata package to read a worksheet in which certain field contains very long character strings (nucleotides sequence, nchar 10,000). Then, the values in these fields are automatically converted to TOO LONG TO DISPLAY. How can I get those original characters instead

[R] Survival model with cross-classified shared frailties

Dear All, The coxph function in the survival package allows multiple frailty terms. In all the examples I saw, however, the frailty terms are nested. What will happen if I have non-nested (that is, cross-classified) frailties in the model? Will the model still work? Do I need to take special

Re: [R] variables from command line

Omar == Omar Lakkis [EMAIL PROTECTED] on Wed, 7 Sep 2005 10:47:43 -0400 writes: Omar How can I pass parameters to an R script from the Omar command line. And how can I read them from within the Omar script? Omar This is how I want to invoke the script: R CMD BATCH Omar

[R] FW: Re: Doubt about nested aov output

Your response nicely clarifies a question that I've had for a long time, but which I've dealt with by giving each subject a unique label. Unless I'm missing something, both techniques should work as the toy example below gives exactly the same output in all 3 cases below (forgetting about the

Re: [R] Too long to display problem

Dear list, Please ignore this thread - the TOO LONG TO DISPLAY is brought by another tool when parsing data sets. Sorry for this ... Wuming On 9/8/05, Wuming Gong [EMAIL PROTECTED] wrote: Dear list, I used read.xls in gdata package to read a worksheet in which certain field contains very

Re: [R] Using R map data to determine associated state for a coordinate?

Werner Wernersen wrote: Hi! I have no idea if this is maybe an easy task utilizing R since I read there is geographical map data in some package: I have a huge number of geographical points with their coordinates in Germany. Now I want to determine for each point in which Bundesland

Re: [R] Interpolating / smoothing missing time series data

On 9/7/05 10:19 PM, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 9/7/05, David James [EMAIL PROTECTED] wrote: The purpose of this email is to ask for pre-built procedures or techniques for smoothing and interpolating missing time series data. I've made some headway on my problem in my

Re: [R] Using R map data to determine associated state for a coordinate?

On Thu, 8 Sep 2005, Thomas Petzoldt wrote: Werner Wernersen wrote: Hi! I have no idea if this is maybe an easy task utilizing R since I read there is geographical map data in some package: I have a huge number of geographical points with their coordinates in Germany. Now I

Re: [R] Prediction with multiple zeros in the dependent variable

On 08-Sep-05 John Sorkin wrote: I have a batch of data in each line of data contains three values, calcium score, age, and sex. I would like to predict calcium scores as a function of age and sex, i.e. calcium=f(age,sex). Unfortunately the calcium scorers have a very ugly distribution. There

[R] Time series ARIMAX and multivariate models

Dear List, The purpose of this e-mail is to ask about R time series procedures - as a biologist with only basic time series knowledge and about a year's experience in R. I have been using ARIMAX models with seasonal components on seasonal data. However I am now moving on to annual data (with

Re: [R] tcltk, X11 protocol error: Bug?

Dear Nicholas, This problem has been reported before (enter X11 protocol error on the R site search at http://finzi.psych.upenn.edu/search.html to see the previous threads), but as far as I know, there's no definitive explanation or solution. As well, things appear to work fine, despite the

[R] Converting a matrix to a dataframe: how to prevent conversion to factor

Colleages I am running R 2.1.0 on a Mac (same problem occurs in Linux). In some situations, I have mixed text/numeric data that is stored as characters in a matrix. If I convert this matrix to a dataframe, the numeric data becomes factors, not what I intend. TEXT- paste(Text,

Re: [R] Prediction with multiple zeros in the dependent variable

John Sorkin wrote: I have a batch of data in each line of data contains three values, calcium score, age, and sex. I would like to predict calcium scores as a function of age and sex, i.e. calcium=f(age,sex). Unfortunately the calcium scorers have a very ugly distribution. There are multiple

[R] Time Series Analysis: book?

There has been a few questions on the subject lately. Is there any book on the subject, if possible with a computer processing flavor, that you would highly recommend? Many thanks in advance, -- Jean-Luc __ R-help@stat.math.ethz.ch mailing list

Re: [R] Converting a matrix to a dataframe: how to prevent conversion to factor

Dennis Fisher [EMAIL PROTECTED] writes: Colleages I am running R 2.1.0 on a Mac (same problem occurs in Linux). In some situations, I have mixed text/numeric data that is stored as characters in a matrix. If I convert this matrix to a dataframe, the numeric data becomes factors,

Re: [R] FW: Re: Doubt about nested aov output

On 9/8/05, Ken Knoblauch [EMAIL PROTECTED] wrote: Your response nicely clarifies a question that I've had for a long time, but which I've dealt with by giving each subject a unique label. Unless I'm missing something, both techniques should work as the toy example below gives exactly the

Re: [R] FW: Re: Doubt about nested aov output

Thank you for your response. The single response/observer most probably explains the complaints that lmer was giving for my example. Maybe this small modification provides a better example and corrects a more serious error in my previous post: library(lme4) y-rnorm(30) cond - rep(gl(3,5,15), 2)

Re: [R] Prediction with multiple zeros in the dependent variable

On Thu, 8 Sep 2005, John Sorkin wrote: I have a batch of data in each line of data contains three values, calcium score, age, and sex. I would like to predict calcium scores as a function of age and sex, i.e. calcium=f(age,sex). Unfortunately the calcium scorers have a very ugly distribution.

Re: [R] Survival model with cross-classified shared frailties

On Thu, 8 Sep 2005, Shige Song wrote: Dear All, The coxph function in the survival package allows multiple frailty terms. Um, no, it doesn't. In all the examples I saw, however, the frailty terms are nested. What will happen if I have non-nested (that is, cross-classified) frailties in

[R] array indices in synced vectors

Let us start with the following definitions xxx-rep(c(1,2),times=5) yyy-rep(c(1,2),each=5) a-c(11,12) b-matrix(1:4,2,2) a[xxx] produces [1] 11 12 11 12 11 12 11 12 11 12 b[xxx,yyy] produces [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,]11111333

Re: [R] array indices in synced vectors

On Thu, 8 Sep 2005, Erich Neuwirth wrote: sapply(1:length(xxx),function(x)b[xxx[x],yyy[x]]) does what I need and produces [1] 1 2 1 2 1 4 3 4 3 4 Is there a function taking xxx,yyy, and b as arguments producing the same result? b[cbind(xxx,yyy)] Essentially, I am asking for a version of

Re: [R] Predicting responses using ace

I'm trying to run the print method, but according to the documentation it needs as a parameter an object created by |summary.areg.boot| . The thing is that |summary.areg.boot| gives me the following error: Error in bootj[, 1] : incorrect number of dimensions, when I do the simple call ---

Re: [R] Predicting responses using ace

I gave a quick read to the documentation again and noticed I misinterpreted it. It was print.summary.areg.boot the method I was referring to (although the summary error should still work). Sorry for the inconvenience Anyway, I used the print method on my |areg.boot| object and I got this:

Re: [R] Prediction with multiple zeros in the dependent variable

John: 1. As George Box long ago emphasized and proved, normality is **NOT** that important in regression, certainly not for estimation and not even for inference in balanced designs. Independence of the observations is far more important. 2. That said, it sounds like what you have here is a

[R] ROracle install problem

Hi, I am trying to install the ROracle package in a Linux-64 machine. I downloaded from Oracle's site their Instant Client bundle but it seems that ROracle needs some stuff not included in that kit in order to compile (in particuar, the 'proc' executable). I did not find any other linux client

[R] package installation error (LF versus CR)

Hello, I have the following problem in installing a package (in windows xp) rcmd install -c dlm [ ..stuff deleted ] ... DLL made installing DLL installing R files installing inst files installing data files installing man source files installing indices Errore in load(zfile,

[R] R API call from delphi

Hello, Has anyone tried to call R API from Delphi under windows ? How was it done if it was ? Or has anyone any idea about how it could be done ? Thanks for your answers Laurent Tessier [[alternative HTML version deleted]] __

Re: [R] Time Series Analysis: book?

TS is a huge topic. The book recomended by statisitcian might be different from the one recommended by econometrician. Finance guy might recommend another. Could you please be more specific? On 9/8/05, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: There has been a few questions on the subject

[R] data manipulation

Dear All, I would be grateful if you can help me. My problem is the following: I have a data set like: ID time X1 X2 11 x111 x211 12 x112 x212 21 x121 x221 22 x122 x222 23 x123 x223 where

Re: [R] package installation error (LF versus CR)

sonia wrote: Hello, I have the following problem in installing a package (in windows xp) rcmd install -c dlm [ ..stuff deleted ] ... DLL made installing DLL installing R files installing inst files installing data files installing man source files installing

Re: [R] Interpolating / smoothing missing time series data

(see inline) Sean Davis wrote: On 9/7/05 10:19 PM, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 9/7/05, David James [EMAIL PROTECTED] wrote: The purpose of this email is to ask for pre-built procedures or techniques for smoothing and interpolating missing time series data. I've made

[R] Tip: I() can designate constants in a regression

Just thought I would share a tip I learned: The function I() is useful for specifying constants to formulas and regressions. It will prevent nls (for example) from trying to treat the variable inside I() as something it needs to estimate. An example is below. -David P.S. This may be

Re: [R] data manipulation

Hi, This may not be the best solution, but at least it's easy to see what i'm doing, assume that your data set is called data: # remove the 4th column data1 = data[,-4] # remove the 3rd column data2 = data[,-3] # use cbind to add an extra column with only X1 #elements data1 = cbind(data1,

Re: [R] data manipulation

Marc Bernard [EMAIL PROTECTED] wrote: Dear All, I would be grateful if you can help me. My problem is the following: I have a data set like: ID time X1 X2 11 x111 x211 12 x112 x212 21 x121 x221 22 x122

Re: [R] Tip: I() can designate constants in a regression

David James [EMAIL PROTECTED] writes: Just thought I would share a tip I learned: The function I() is useful for specifying constants to formulas and regressions. It will prevent nls (for example) from trying to treat the variable inside I() as something it needs to estimate. An

[R] Multinomial Logit and p-values

Hi, I am trying to obtain p-values for coefficient estimates in a multinomial logit model. Although I am able to test for significance using other methods (e.g., Wald statistics), I can't seem to get R to give me simple p-values. I am sure there is a very simple solution to this, but the R

Re: [R] data manipulation

This is what reshape() does. -thomas On Thu, 8 Sep 2005, Marc Bernard wrote: Dear All, I would be grateful if you can help me. My problem is the following: I have a data set like: ID time X1 X2 11 x111 x211 12 x112 x212 21

[R] execute R expression from command line

Can I execute an R expression from the command line without having it in an infile, something like perl's -e flag. So it would look like: R {Rexpression;} outfile __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] data manipulation

Also see Hadley Wickham's reshape package for more bells whistles. -- HTH! Jim Porzak Loyalty Matrix Inc. On 9/8/05, Thomas Lumley [EMAIL PROTECTED] wrote: This is what reshape() does. -thomas On Thu, 8 Sep 2005, Marc Bernard wrote: Dear All, I would be grateful if

[R] Re-evaluating the tree in the random forest

Dear mailinglist members, I was wondering if there was a way to re-evaluate the instances of a tree (in the forest) again after I have manually changed a splitpoint (or split variable) of a decision node. Here's an illustration: library(randomForest) forest.rf - randomForest(formula = Species ~

Re: [R] R API call from delphi

On approach is to create a native/foreign interface to R by linking R as a library (libR.a and R.dll) file and calling the C routines in the library to i) initialize the R interpreter ii) call an R function We have done this with many languages and the procedure is well understood at this

Re: [R] Predicting responses using ace

Luis Pineda wrote: I gave a quick read to the documentation again and noticed I misinterpreted it. It was print.summary.areg.boot the method I was referring to (although the summary error should still work). Sorry for the inconvenience Anyway, I used the print method on my |areg.boot|

Re: [R] data manipulation

I am sure all this work but If you want exaclty the output to be the way you mentioned do this temp-read.table(yourfile, as.is=T, header=T) temp1-temp[, 1:3] temp2-temp[, c(1,2,4)] colnames(temp1)[3]-X colnames(temp2)[3]-X temp3-merge(temp1, temp2, all=T) temp3$type-toupper(substr(temp3$X, 1,2))

[R] writing data to sheet in excel workbook

Hi, I believe to remember there is a package that lets you write data from R to different sheets in a Excel workbook. I've been looking around on CRAN but could not find what I am looking for. Any help would be greatly appreciated. Cheers, Adi __

Re: [R] Leading in line-wrapped Lattice value and panel labels

On 9/7/05, Paul Murrell [EMAIL PROTECTED] wrote: Hi Deepayan Sarkar wrote: On 9/7/05, Tim Churches [EMAIL PROTECTED] wrote: Version 2.1.1 Platforms: all What is the trellis parameter (or is there a trellis parameter) to set the leading (the gap between lines) when long axis

Re: [R] execute R expression from command line

On 8 Sep 2005, [EMAIL PROTECTED] wrote: Can I execute an R expression from the command line without having it in an infile, something like perl's -e flag. So it would look like: R {Rexpression;} outfile With a bash-like shell, you can do: echo library(foo); somefunc(5) | R --slave HTH,

Re: [R] Survival model with cross-classified shared frailties

Hi Thomas, Thanks for the reply, coxme() seems to be the one I need. Best, Shige On 9/8/05, Thomas Lumley [EMAIL PROTECTED] wrote: On Thu, 8 Sep 2005, Shige Song wrote: Dear All, The coxph function in the survival package allows multiple frailty terms. Um, no, it doesn't.

[R] Setting width in batch mode

As instructed, I have spent a long time searching the web for an answer to this question. I am trying to use Sweave to produce lecture slides, and have the problem that I can't control the formatting of my R source. Setting options(width), as recommended in this forum, works fine on the R

Re: [R] writing data to sheet in excel workbook

On 9/8/05, adalbert duerrer [EMAIL PROTECTED] wrote: Hi, I believe to remember there is a package that lets you write data from R to different sheets in a Excel workbook. I've been looking around on CRAN but could not find what I am looking for. See

Re: [R] R API call from delphi

Follow this thread http://finzi.psych.upenn.edu/R/Rhelp02a/archive/50598.html Cheers Francisco From: Laurent TESSIER [EMAIL PROTECTED] Reply-To: [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Subject: [R] R API call from delphi Date: Thu, 8 Sep 2005 17:47:49 +0200 (CEST) Hello, Has anyone

[R] can't successfully use installed evir package

I'm next at installing packages. I seem to have successfully installed evir, but I can't use it. I'm wondering if I need to specify the installation to match my working directory, or something else. thx, G __ R-help@stat.math.ethz.ch mailing list

[R] clustering: Multivariate t mixtures

Hi, Before I write code to do it does anyone know of code for fitting mixtures of multivariate-t distributions. I can't use McLachan's EMMIX code because the license is For non commercial use only. I checked, mclust and flexmix but both only do Gaussian. Thanks Nicholas

[R] generating a vector from clusters of logicals

dear friends, I have a vector of clusters of TRUE and FALSE like c(TRUE,TRUE,TRUE...,FALSE,FALSE,FALSE,TRUE,TRUE...) and want to make that into a vector of c(1,1,1,1...2,2,2,2,.3,3,3,3) increasing the number assigned to each cluster as they change. How would I do that ? Best wishes

Re: [R] clustering: Multivariate t mixtures

On Thu, 08 Sep 2005 15:38:55 -0500 Nicholas Lewin-Koh wrote: Hi, Before I write code to do it does anyone know of code for fitting mixtures of multivariate-t distributions. I can't use McLachan's EMMIX code because the license is For non commercial use only. I checked, mclust and flexmix

Re: [R] generating a vector from clusters of logicals

On Thu, 08 Sep 2005 23:03:03 +0200 Troels Ring wrote: dear friends, I have a vector of clusters of TRUE and FALSE like c(TRUE,TRUE,TRUE...,FALSE,FALSE,FALSE,TRUE,TRUE...) and want to make that into a vector of c(1,1,1,1...2,2,2,2,.3,3,3,3) increasing the number assigned to each

Re: [R] clustering: Multivariate t mixtures

Hi, Actually that was my plan was to implement a new flexmix class. Thanks for the pointer to the jss paper, that will be helpful. Nicholas On Thu, 8 Sep 2005 23:07:13 +0200, Achim Zeileis [EMAIL PROTECTED] said: On Thu, 08 Sep 2005 15:38:55 -0500 Nicholas Lewin-Koh wrote: Hi, Before I

Re: [R] Time Series Analysis: book?

Wensui Liu wrote: TS is a huge topic. The book recomended by statisitcian might be different from the one recommended by econometrician. Finance guy might recommend another. Could you please be more specific? My software (http://moodss.sourceforge.net) collects, archives in a SQL database

Re: [R] generating a vector from clusters of logicals

From: Troels Ring [EMAIL PROTECTED] I have a vector of clusters of TRUE and FALSE like c(TRUE,TRUE,TRUE...,FALSE,FALSE,FALSE,TRUE,TRUE...) and want to make that into a vector of c(1,1,1,1...2,2,2,2,.3,3,3,3) increasing the number assigned to each cluster as they change. How

Re: [R] generating a vector from clusters of logicals

Try: x - c(TRUE,TRUE,TRUE,FALSE,FALSE,FALSE,TRUE,TRUE) tmp - rle(x) tmp$values - seq(along=tmp$lengths) new.x - inverse.rle(tmp) new.x Greg Snow, Ph.D. Statistical Data Center, LDS Hospital Intermountain Health Care [EMAIL PROTECTED] (801) 408-8111 Troels Ring [EMAIL PROTECTED] 09/08/05

[R] Data Expo 2006 (off-topic)

Hi This is to let R folks know about the Data Expo that is being run by the ASA Sections on Statistical Graphics, Statistical Computing, and Statistics and the Environment for JSM 2006. This competition provides a data set of geographic and atmospheric data from NASA and entrants are asked to

Re: [R] ROracle install problem

Ariel Chernomoretz [EMAIL PROTECTED] wrote: Hi, I am trying to install the ROracle package in a Linux-64 machine. I downloaded from Oracle's site their Instant Client bundle but it seems that ROracle needs some stuff not included in that kit in order to compile (in particuar, the 'proc'

[R] Re: General Matrix Inner Product?

Roger E. Khayat, Professor Department of Mechanical and Materials Engineering The University of Western Ontario London, Ontario, Canada N6A 5B9 Email: [EMAIL PROTECTED] Tel: (519) 661-2111 Ext 88253 Fax: (519) 661-3020 http://www.engga.uwo.ca/people/rkhayat/

[R] Coarsening Factors

It is not uncommon to want to coarsen a factor by grouping levels together. I have found one way to do this in R: sites [1] F A A D A A B F C F A D E E D C F A E D F C E D E F F D B C Levels: A B C D E F regions - list(I = c(A,B,C), II = D, III = c(E,F)) library(Epi) region -

Re: [R] Coarsening Factors

Murray Jorgensen [EMAIL PROTECTED] writes: It is not uncommon to want to coarsen a factor by grouping levels together. I have found one way to do this in R: sites [1] F A A D A A B F C F A D E E D C F A E D F C E D E F F D B C Levels: A B C D E F regions - list(I = c(A,B,C), II =

[R] change in read.spss, package foreing?

Dear All, it seems to me that the function read.spss of package foreign changed its behaviour regarding factors. I noted that in version 0.8-8 variables with value labels in SPSS were transformed in factors with the labels in alphabetic order. In version 0.8-10 they seem to be ordered preserving

Re: [R] Time Series Analysis: book?

1. Have you read the appropriate chapter in Venables and Ripley (2002) Modern Applied Statists with S (Springer)? If no, I suggest you start there. 2. Have you worked through the vignettes associated with the zoo package? If no, you might find that quite useful. [Are

[R] SPSS Dataset

How would one read SPSS data sets directly into R __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Time Series Analysis: book?

Quoting Spencer Graves [EMAIL PROTECTED]: 1. Have you read the appropriate chapter in Venables and Ripley (2002) Modern Applied Statists with S (Springer)? If no, I suggest you start there. 2. Have you worked through the vignettes associated with the zoo package? If no,

[R] strata in crr (cmprsk library)

Hi all, I am aware that crr lacks the friendly command structure of functions such as cph. All is clear to me about including covariates until I want to include a stratification term in the competing risk framework (no nice strat command). I am still a bit of a novice in R - I am looking for an

[R] Debugging R/Fortran in Windows

Hi, I'm trying to debug an R interface to a Fortran subroutine from Windows. (Yes, I know I should try Unix/Linux as well, but a quick attempt suggested that the (MinGW g77) Fortran compiler I have installed on my Windows laptop works better on this Fortran code.) I'm trying to follow the

Re: [R] SPSS Dataset

RSiteSearch(read spss data) -- library(foreign) ?read.spss Best, Matthias -Ursprüngliche Nachricht- Von: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Im Auftrag von ERICK YEGON Gesendet: Freitag, 09. September 2005 07:02 An: R-help@stat.math.ethz.ch Betreff: [R] SPSS Dataset