[R] Decompose an irregular daily time series with missing values

Dear all, I have a daily time series for the months April to September (183 says) for a 40 years period. It contain missing values. I would like to extract a seasonal component, trend component and irregular component using an adaptive model. #Commands for making a time series:

[R] Heatmap clustering

I have constructed a heatmap comparing distances between 30 elements. I get the heatmap plot with red dots when distance is close, and yellow when two elements differ. My question is...how can I get clusters of similarity from the headplot (I don´t know the total number of clusters, otherwise

[R] correlation structure in lmer

Hi how can I specify a correlation structure in the lmer-function as it is possible in lme(formula, ..., corr=corAR1(form=...))? Thanks Fränzi __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] Simple Graph

Hi Sir We use the function mpg() for ticks on axes. If we define more number of ticks then each tick is not labeled. How to label each tick on axes? Regards -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. Email: [EMAIL PROTECTED] [EMAIL PROTECTED] [EMAIL PROTECTED]

[R] Single plot multiple levels in x?

Plotting with 2 x axis? One axis inside another, for example salary within state, 1-50 | 50 – 100 | 100+ | 1- 50 | 50 -100 | 100+ | … repeated bins for salary AL ! AR …… more states The values are

[R] rpart's loss matrix in ipred

Dear R users, I have been using the rpart procedure to predict the occurrence of depression in a large data file. Since the prevalence is very low (5%), I have been using classification trees with a loss matrix that penalized false negatives more than false positives. I have become

[R] Converting into time series object

Hi, I have a dataframe of trading dates along with the corresponding prices. I need to convert this into a time series object. How do I do this with my price values being the time series object and the dates/time being the trading dates. BR, Shubha [[alternative HTML version

[R] Q: Mean, median and confidence intervals with functions summary boxplot.stats

Dear R ussers, My question is, How can my mean be outside the confidence intervals ?! I think i have the answer for it, but i would like to hear some other ideas on it. First my data is not continuose but categorical, it is a titre calculated on a dilution serie. It is stored as a column of

[R] Reading lines from file

Hi. I have a text file containing a few hundred lines of numbers, each line has a different length. For example: 1 4 1 1 7 3 11 1 1 1 1 1 1 2 2 4 1 2 And so on. I need to do a simple plot function for each line, and then to save each plot into a separate file. Is there any way doing it from

Re: [R] Month end calculations

A simple example (avoiding using dates, just to show the principle) - this assumes that your data are already sorted (?order). temp - data.frame(subject = rep(1:2, each = 5), response = 1:10) print(temp) last - do.call(rbind, by(temp, temp$subject, function(x) tail(x, 1))) print(last) By

[R] How to signal the end of the table?

I am using a for loop to read a table row by row and I have to specify how many records are there in the table. I need to read row by row because the table is huge and the memory not large enough for the whole table.: number.of.records=100 fp=file(abc.csv,r) pos=seek(fp, rw=read) for (i in

Re: [R] customizing the color and point shape for each line drawn using lattice's xyplot

Does this help common.without.Method4 - subset(common, Method!=4) xyplot(MeanBxg ~ PercentVarExplained | bdg.f * bdx.f, data=common.without.Method4, groups=Method.f, type=l, auto.key=T) Ross Darnell -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gen

[R] Incomplete output with `sn' library package

Dear R users: I have a question regarding the output for two of the functions in the `sn' package, which deals with the mle fitting of skew normal curves to linear regressions. I'm using the examples and the dataset given as an example in the online documentation for this package, for the

Re: [R] Reading lines from file

uv wrote: Hi. I have a text file containing a few hundred lines of numbers, each line has a different length. For example: 1 4 1 1 7 3 11 1 1 1 1 1 1 2 2 4 1 2 And so on. I need to do a simple plot function for each line, and then to save each plot into a separate file. Is there any way

Re: [R] Excel

On 30/8/07 6:42 AM, Erich Neuwirth wrote: There is one feature in Excel which is extremely convenient, Pivot tables. Anybody doing any work as statistical consultant really ought to know about Pivot tables, and I am still surprised how many statisticians do not know about it. Neither Gnumeric

Re: [R] Incomplete output with `sn' library package

MANASI VYDYANATH wrote: Dear R users: I have a question regarding the output for two of the functions in the `sn' package, which deals with the mle fitting of skew normal curves to linear regressions. I'm using the examples and the dataset given as an example in the online

Re: [R] Strage result with an append/strptime combination

Hi again, Just a quick post to propose another solution (without chron) : strptime function on the whole object and not on each namefile. namefile-070707050642.dat namefile-append(namefile,namefile) namefile-append(namefile,namefile) namefile [1] 070707050642.dat 070707050642.dat

Re: [R] Simple Graph

amna khan wrote: Hi Sir We use the function mpg() for ticks on axes. If we define more number of ticks then each tick is not labeled. How to label each tick on axes? What is mpg()? Which version of R? Which device? Which OS? I guess R cannot print more labels because the ticks are too

Re: [R] Q: Mean, median and confidence intervals with functions summary boxplot.stats

Tom Willems wrote: Dear R ussers, My question is, How can my mean be outside the confidence intervals ?! I think i have the answer for it, but i would like to hear some other ideas on it. First my data is not continuose but categorical, it is a titre calculated on a dilution

Re: [R] How to signal the end of the table?

What about trying to use a database system and make queries? If you are at the end of a file, you will read elements of length 0... Uwe Ligges Yuchen Luo wrote: I am using a for loop to read a table row by row and I have to specify how many records are there in the table. I need to read row

Re: [R] Single plot multiple levels in x?

Richard Yanicky wrote: Plotting with 2 x axis? One axis inside another, for example salary within state, 1-50 | 50 – 100 | 100+ | 1- 50 | 50 -100 | 100+ | … repeated bins for salary AL ! AR

[R] Simple Graph

Hi Sir If I want 10 ticks on x-axis and 10 ticks on y-axis then I use lab=c(10,10, 20) Here first 10 is for 10 ticks on x-axis, second 10 is for 10 ticks on y-axis and 20 is the lenght of x and y labels in characters. Now my question is that how I can label all 10 ticks on x-aixs and all 10

[R] (no subject)

Dear R users,I am looking for an easy (i.e., direct) way in obtaining the F and p values from the intercept in anovas with within-subject designs. Data are from a psychophysics expeirment where I am using d' (d-prime) values obtained from 3 modalities of presentation in each subject.I would

Re: [R] Simple Graph

amna khan wrote: Hi Sir If I want 10 ticks on x-axis and 10 ticks on y-axis then I use lab=c(10,10, 20) Here first 10 is for 10 ticks on x-axis, second 10 is for 10 ticks on y-axis and 20 is the lenght of x and y labels in characters. Now my question is that how I can label all 10

[R] how to calculate weighted means by groups

Hello, I want to calculate the weighted means of a table column, aggregated by other elements of the same table. For example I want to calculate mean species numbers weighted by area and grouped by year, altitude and/or region. I think a combination of aggregate() and weighted.mean() should work,

Re: [R] How to signal the end of the table?

Have you tried read.table() and the similar functions? I think they would provide a much more simple solution for your problem. HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research

[R] Data simulation with R

Hi, I am currently on a placement here at GSK for my studies, and I'm working on Heckman Models. I have to make simulations, in order to see whether these models are efficient or not. I have to generate a dataset under the following constraints: - outcome is 0 or 1 - one control group, one

[R] R CMD BATCH: cat does not print

Dear All, I am trying to write my first R script. The code is simply cat(Hello!\n) However, when I run $ R CMD BATCH myscript.R I do not see Hello! on the console. I am using Fedora 7 (Linux) and R-2.5.1. Any ideas? Thanks in advance, Paul __

[R] bug in DEoptim package

(the same mail was sent to the author) When I called the function DEoptim with control=list(strategy=1) or control=list(strategy=2) I got the error: Error in mui[rtd + 1, i] : incorrect number of dimensions Analysis of the source code of the DEoptim reveals the following fragment if

Re: [R] R CMD BATCH: cat does not print

Paul Smith wrote: Dear All, I am trying to write my first R script. The code is simply cat(Hello!\n) However, when I run $ R CMD BATCH myscript.R I do not see Hello! on the console. I am using Fedora 7 (Linux) and R-2.5.1. Any ideas? You shouldn't see it on the console!

Re: [R] breaking the x-axis and having two different x-axis labels

Georg Ehret wrote: Dear R community, I have two questions concerning barplots that I struggle to resolve: 1) How can I break (interrupt) the x-axis (e.g.: have it display values from -100 to -90 and 90 to 100 only)? I think you mean a horizontal barplot with a gap between -90 and 90.

Re: [R] R CMD BATCH: cat does not print

Use rscript Rscript myscript.R or Rscript -e 'cat(Hello!\n)' will show Hello! on the console. R CMD BATCH writes its output to the file myscript.Rout Paul Smith wrote: Dear All, I am trying to write my first R script. The code is simply cat(Hello!\n) However, when I run $ R

Re: [R] R CMD BATCH: cat does not print

On 8/30/07, Barry Rowlingson [EMAIL PROTECTED] wrote: I am trying to write my first R script. The code is simply cat(Hello!\n) However, when I run $ R CMD BATCH myscript.R I do not see Hello! on the console. I am using Fedora 7 (Linux) and R-2.5.1. Any ideas? You

Re: [R] R CMD BATCH: cat does not print

On 8/30/07, Vladimir Eremeev [EMAIL PROTECTED] wrote: Use rscript Rscript myscript.R or Rscript -e 'cat(Hello!\n)' will show Hello! on the console. R CMD BATCH writes its output to the file myscript.Rout Thanks, Vladimir. Rscript is exactly what I was looking for! Paul Paul Smith

Re: [R] Simple Graph

amna khan wrote: Hi Sir We use the function mpg() for ticks on axes. If we define more number of ticks then each tick is not labeled. How to label each tick on axes? Hi Amina, The staxlab function in the plotrix package will display tick labels even if they would overlap on the standard

Re: [R] piecewise linear approximation

If you want to minimize absolute error for this, then you can try the rqss fitting in the quantreg package and tune lambda to get one break in the fitted function. url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox:

[R] Assigning line colors in xyplot

Hi, I have a dataframe containing data from individuals 1, ..., 12 (grouping variable g in the data frame below), which belong either to A or B (grouping variable f): set.seed(1) tmp - data.frame(

Re: [R] Converting into time series object

You can use the 'ts' function. Example df - ts(your_data_frame, start=c(1990,1), #The initial date of your data end=c(2000,12), #The final date of your data or frequency=12)#The frequency of your data, so you don't need of the 'end' argument See ?ts and ?as.ts

Re: [R] Q: Mean, median and confidence intervals with functions summary boxplot.stats

If you look at ?boxplot.stats, you will find that the confidence interval it reports is centred on the median and : The notches (if requested) extend to '+/-1.58 IQR/sqrt(n)'. If you have skewed data it is very possible (as you have found) that the mean is outside median+/-1.58 IQR/sqrt(n).

Re: [R] Excel (off-topic, sort of)

Try prolog. You can do this sort of programming there. Albicelli, Nicholas (Exchange) wrote: Except for the ability to perform circular recalculation, I believe that the closest programming analogy to a spreadsheet is a functional programming language. Check out Haskell (or LISP or

[R] boxplot will remember the factor levels

R-help, I'm trying to do a simple box-and-whisker plot to some data. The data are a subset of a large data frame but when running the boxplot function on the subset data all the factors are still present in the graph leaving a huge empty space until the actuals factors are shown. This produces a

Re: [R] R CMD BATCH: cat does not print

Paul Smith wrote: Thanks, Barry. Indeed, the file myscript.Rout exists and contains the output of cat. I was expecting a behavior similar to the bash scripts. And by the way, cannot a R script write only on the console and just what one tells it to write, likewise bash scripts? Not easily,

[R] How to mask or escape =

-- View this message in context: http://www.nabble.com/How-to-mask-or-escape-%22%3D%22-tf4354174.html#a12406926 Sent from the R help mailing list archive at Nabble.com. __ R-help@stat.math.ethz.ch mailing list

[R] How to mask or escape = in Windows command prompt?

I have defined a function with several arguments and have it stored in the .RData file. The 'function head' is defined as follows EstimALIConc -function(sdname,SZ,W,farea,watri,biomodel,start.part=1,nparts=20,method=c(optim,DEoptim)) { [ blah-blah-blah ] (function body doesn't matter) }

Re: [R] boxplot will remember the factor levels

Try: boxplot(mLength ~ puntar[drop=T],data=test) -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 30/08/2007, Luis Ridao Cruz [EMAIL PROTECTED] wrote: R-help, I'm trying to do a simple box-and-whisker plot to some data. The data are a subset of a large data

Re: [R] xeon processor and ATLAS

hui xie [EMAIL PROTECTED] writes: The reason why I would like to use ATLAS is that R FAQ said : The savings can be appreciable: on a 2.6GHz P4 and a 1000 x 1000 matrix svd took 16.2 sec with the standard BLAS and 7.8 sec with ATLAS. Because ATLAS is tuned to a particular chip we can't

[R] Another issue with the Matrix package.

Hi all, I am encountering a strange issue with the Matrix package. I have just built R-devel from source on my macbook pro, and I wonder if others can reproduce this problem. I will give example code to go along: Starting a fresh R session: R version 2.6.0 Under development (unstable)

[R] categorical variable coefficients in QSAR

Dear list: I am interested in the following sort of problem, as is found frequently in the field of QSAR. I have biological activity as a function of chemical structure, with structure defined in a categorical manner in that the SUBSTITUENT is the levels of the POSITION factor. For example, data

[R] categorical variable coefficients in QSAR

Dear list: I am interested in the following sort of problem, as is found frequently in the field of QSAR. I have biological activity as a function of chemical structure, with structure defined in a categorical manner in that the SUBSTITUENT is the levels of the POSITION factor. For example, data

[R] How to signal the end of the table?

Well, i am surprise you have problems to read a table that is small enough to be opened entirely in Excel. I work with csv tables with hundreds of thousands of rows, and sometimes even millions with no problems except that sometimes i have to wait up to 1 or 2 minutes for R to read the

Re: [R] Single plot multiple levels in x?

Uwe, I have looked into lattice and can't seem to make this work. I can easily make multiple panels but this isn't what I am looking to do. Any suggestions on which functions to use? the axis function seems a natural place to start but I still can't seem to make it happen. HELP! Richard

Re: [R] Month end calculations

Shubha, I apologize if this is a bit late - consequence of digest summary preference. If I understand what you need, it is to calculate the value at month end given a data object. For zoo objects the following should do what you need. Actually is part of my new package on CRAN quantmod -

Re: [R] xeon processor and ATLAS

On Thu, Aug 30, 2007 at 09:33:32AM -0400, Jeffrey J. Hallman wrote: saving 8 seconds of computation time is worth going through the ATLAS hassle? For some of us, the 'hassle' is typing $ sudo apt-get install atlas3-base which gives you the low-hanging fruit of basic Atlas. Locally

Re: [R] retrieve p-value from a cox.obj

Several responders have given input on how to find and read the code for summary.coxph and print.coxph, or using str() to examine the output of a coxph object. An even better starting place would be the documentation (yes there is some). help(coxph) gives a page that contains: Value:

Re: [R] Another issue with the Matrix package *under R-devel*

I suspect you have not using a re-installed Matrix after re-building R. I can reproduce the problem using a version of Matrix I installed under 2.5.1, but not with one installed under R-devel this week. Since R-devel is 'Under development' you may need to reinstall packages when it changes.

Re: [R] How to signal the end of the table?

For example, fp-file(abc.csv,r) c.row-scan(file=fp,sep=,,nlines=1) # what argument is omitted for bevity as it doesn't matter rows-c.row while(length(c.row)0) { c.row-scan(file=fp,sep=;,nlines=1); rows-rbind(rows,c.row) } close(fp) If you want to read a file by parts, then you do

[R] Behaviour of very large numbers

Dear all, I am struggling to understand this. What happens when you raise a negative value to a power and the result is a very large number? B [1] 47.73092 -51^B [1] -3.190824e+81 # seems fine # now this: x - seq(-51,-49,length=100) x^B [1] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN

[R] R and Web Applications

Hello, I'm curious to know how people are calling R from web applications (I've been looking for Perl but I'm open to other languages). After doing a search, I came across the R package RSPerl, but I'm having difficulties getting it installed (on Mac OSX). I believe the problem probably has to

Re: [R] correlation structure in lmer

On 8/29/07, Fränzi Korner [EMAIL PROTECTED] wrote: how can I specify a correlation structure in the lmer-function as it is possible in lme(formula, ..., corr=corAR1(form=...))? The short answer is you can't. __ R-help@stat.math.ethz.ch mailing list

Re: [R] Assigning line colors in xyplot

Hi Christof, You can do this in ggplot, with one exception: install.packages(ggplot2) library(ggplot2) qplot(x, y, data=tmp, facets = . ~ g, geom=line, colour=f) Unfortunately I don't yet have an implementation of facetting that works like lattice, wrapping the line of plots in to 2d

Re: [R] Excel (off-topic, sort of)

On Wed, 29 Aug 2007, Alberto Monteiro wrote: Do you know what's in my wish list? I wish spreadsheets and computer languages had gone one step further. I mean, it's nice to define Cell X to be equal to Cell Y + 10, and then when we change Cell Y, magically we see Cell X change. But why

[R] mtext, picking subscript text from colnames(df)

Hi! Could someone give advise on how to plot Greek letters where subscript text is picked from the names of a data.frame. I want what appears using this code mtext(expression(mu[Mountain]),side=2,at=max(y)-7,las=1,cex=1.5) but I would like to do it using something like this:

Re: [R] Month end calculations

Shubha et al, I forgot that the period methods may not work without the full package. So either just download the current CRAN version 0.1-0 or version 0.1-3 from quantmod.com, a new more complete version will be up with a week or two. For completeness here is my months method for zoo objects:

Re: [R] xeon processor and ATLAS

On 8/29/07, hui xie [EMAIL PROTECTED] wrote: Thanks very much for all your advice. To be clear, my OS is window XP. I bought this server last year. It's Dell Precision PWS690. THe processor is Xeon(TM) CPU 3GHZ, 2G RAM. I am not sure how to check more details of processor on my computer.

[R] t-test within tapply

Hello, My data are as following: Data - data.frame(Ind=rep(1:3,c(10,10,10)), Replicate=rep(c(rep(a,5),rep(b,5)),3), EggSize=rep(rnorm(5,mean=10),6) ) attach(Data) Using a t-test, I want to check if the mean egg sizes are

Re: [R] Single plot multiple levels in x?

Richard Yanicky wrote: Uwe, I have looked into lattice and can't seem to make this work. I can easily make multiple panels but this isn't what I am looking to do. Any suggestions on which functions to use? the axis function seems a natural place to start but I still can't seem to make

Re: [R] Behaviour of very large numbers

Hello, it seems to be an R bug. It gives strange errors for non-integer exponents: version platform i686-redhat-linux-gnu version.string R version 2.4.1 (2006-12-18) x^47.0 [1] -1.802180e+80 -1.768932e+80 -1.736284e+80 -1.704227e+80 -1.672748e+80 [...] x^47.10 [1] NaN NaN NaN NaN

[R] How to obtain intercept statistics in anova with within-subject factors?

Dear R users, I am looking for an easy (i.e., direct) way of obtaining the F and p values from the intercept in anovas with within-subject designs. My data are from a psychophysics experiment where I am using d' (d-prime) values obtained from 3 modalities of presentation in each subject. I would

Re: [R] Excel

On 8/28/2007 3:16 AM, J Dougherty wrote: On Monday 27 August 2007 22:21, David Scott wrote: On Tue, 28 Aug 2007, Robert A LaBudde wrote: If you format the column as Text, you won't have this problem. By leaving the cells as General, you leave it up to Excel to guess at the correct

Re: [R] Incomplete output with `sn' library package

You have my sincere apologies for the incompleteness of my message. I have given the details below, including my dataset and my code. I'm using R, version 2.5.0. My OS is a Mac, (version Tiger). The sn package is Version 0.4-1 My code was as follows: mydata -

Re: [R] R and Web Applications

The R packages and projects for the web and R are listed here: http://www.lmbe.seu.edu.cn/CRAN/doc/FAQ/R-FAQ.html#R-Web-Interfaces On 8/30/07, Chris Parkin [EMAIL PROTECTED] wrote: Hello, I'm curious to know how people are calling R from web applications (I've been looking for Perl but I'm

Re: [R] Behaviour of very large numbers

On 8/30/2007 11:08 AM, willem vervoort wrote: Dear all, I am struggling to understand this. What happens when you raise a negative value to a power and the result is a very large number? B [1] 47.73092 -51^B [1] -3.190824e+81 You should be using parentheses. You evaluated -(51^B),

Re: [R] customizing the color and point shape for each line drawn using lattice's xyplot

Thank you for your suggestion. Unfortunately it does not appear to matter if the data is subsetted before the xyplot command or within it. The problem remains. I was able to remove method 4 entirely from the key by using your suggestion to subset the dataset “common,” and then creating a new

Re: [R] R and Web Applications

Hi Chris -- RWebServices provides a way to expose R functionality as (SOAP-based) web services. http://bioconductor.org/packages/2.0/bioc/html/RWebServices.html This might be more than you are looking for, and has not been tested on MacOS (which seems to be your platform). Martin Chris Parkin

Re: [R] Q: how to interrupt long calculation?

Paul Smith said the following at 08/29/2007 04:32 PM : The instance of R running will be immediately killed and then you can start R again. But then I would lose all the work. There must be some way to merely interrupt the current calculation. Mustn't there?

Re: [R] correlation structure in lmer

On 8/30/07, Douglas Bates [EMAIL PROTECTED] wrote: On 8/29/07, Fränzi Korner [EMAIL PROTECTED] wrote: how can I specify a correlation structure in the lmer-function as it is possible in lme(formula, ..., corr=corAR1(form=...))? The short answer is you can't. My response was incorrect. I

Re: [R] Behaviour of very large numbers

On 8/30/07, Gustaf Rydevik [EMAIL PROTECTED] wrote: On 8/30/07, willem vervoort [EMAIL PROTECTED] wrote: Dear all, I am struggling to understand this. What happens when you raise a negative value to a power and the result is a very large number? B [1] 47.73092 -51^B [1]

Re: [R] Excel (off-topic, sort of)

On 8/30/2007 10:43 AM, Thomas Lumley wrote: On Wed, 29 Aug 2007, Alberto Monteiro wrote: Do you know what's in my wish list? I wish spreadsheets and computer languages had gone one step further. I mean, it's nice to define Cell X to be equal to Cell Y + 10, and then when we change Cell Y,

Re: [R] Behaviour of very large numbers

-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of willem vervoort Sent: Thursday, August 30, 2007 8:09 AM To: r-help@stat.math.ethz.ch Subject: [R] Behaviour of very large numbers Dear all, I am struggling to understand this. What happens when

Re: [R] Nested functions.

nalluri pratap schreef: Hi All, I have two variables X, Y. The question is if the value of X is equal to one, then the values in Y have to be reversed other wise it should not perfom any action. I think this should be done using lapply function? Example Y values : 1 2

Re: [R] Behaviour of very large numbers

willem vervoort wrote: Dear all, I am struggling to understand this. What happens when you raise a negative value to a power and the result is a very large number? B [1] 47.73092 -51^B [1] -3.190824e+81 # seems fine Well, this seems not to be what you intended to do,

[R] xyplot() groups scope issue

Hi all, Tinkering with a wrapper for xyplot that will help me plot a bunch of plots in a data analysis I'm doing and I ran into an odd error that I'm guessing is a scope issue. Here's a very simple version of the code: ### #load lattice library(lattice) #create the wrapper

Re: [R] R and Web Applications

Dear Chris, I use Python (http://www.python.org) in combination with Rpy (http://rpy.sourceforge.net/). Rpy enables you to use R commands inside Python, not the other way around. Works quite well, also for different R versions (I currently run R 2.5.1 under Linux). cheers, Paul Chris Parkin

Re: [R] Behaviour of very large numbers

On Thu, 30 Aug 2007, willem vervoort wrote: Dear all, I am struggling to understand this. What happens when you raise a negative value to a power and the result is a very large number? Where are the 'very large numbers' here? R can cope with much larger numbers (over 10^300). B [1]

[R] 7 Courses: Upcoming September-October 2007 R/S+ schedule by XLSolutions Corp

Our upcoming September-October 2007 R/S+ course schedule is now available. Please check out this link for additional information and direct enquiries to Sue Turner [EMAIL PROTECTED] Phone: 206 686 1578 http://www.xlsolutions-corp.com/courselist.htm

Re: [R] mtext, picking subscript text from colnames(df)

Try: mtext(eval(parse(tex=paste(expression(mu[, colnames(PresEsts)[2],]), sep=))),side=2,at=max(y)-15,las=1) -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O [[alternative HTML version deleted]] __

Re: [R] Behaviour of very large numbers

On 8/30/2007 12:11 PM, Martin Becker wrote: willem vervoort wrote: Dear all, I am struggling to understand this. What happens when you raise a negative value to a power and the result is a very large number? B [1] 47.73092 -51^B [1] -3.190824e+81 # seems fine Well,

Re: [R] xyplot() groups scope issue

On 8/30/07, Mike Lawrence [EMAIL PROTECTED] wrote: Hi all, Tinkering with a wrapper for xyplot that will help me plot a bunch of plots in a data analysis I'm doing and I ran into an odd error that I'm guessing is a scope issue. Here's a very simple version of the code: It's indeed a scoping

Re: [R] Assigning line colors in xyplot

On 8/30/07, Christof Bigler [EMAIL PROTECTED] wrote: Hi, I have a dataframe containing data from individuals 1, ..., 12 (grouping variable g in the data frame below), which belong either to A or B (grouping variable f): set.seed(1) tmp - data.frame(

[R] percentage explained by fixed effects in random model

Hi, I realise this has come up before in various reincarnations but I couldnt find the answer... I wish to quote the percentage variance explained by each of three components in my mixed model. If I didnt have a random effect I would just use r squared. I can work out the percentage explained

Re: [R] Month end calculations

The zoo package includes the yearmon class to facilitate such manipulations. Here are a few solutions assuming you store you series in a zoo variable: # test data library(zoo) z - zoo(1001:1100, as.Date(101:200))[-(45:55)] # Solution 1. tapply produces indexes of last of month tt - time(z) z[

Re: [R] Month end calculations

The last line is wrong (see below for correction): On 8/30/07, Gabor Grothendieck [EMAIL PROTECTED] wrote: The zoo package includes the yearmon class to facilitate such manipulations. Here are a few solutions assuming you store you series in a zoo variable: # test data library(zoo) z -

Re: [R] Q: how to interrupt long calculation?

On Thu, 30 Aug 2007, D. R. Evans wrote: Paul Smith said the following at 08/29/2007 04:32 PM : The instance of R running will be immediately killed and then you can start R again. But then I would lose all the work. There must be some way to merely interrupt the current calculation.

[R] Barplot2 using for loop, how to adjust margins?

Hi R-users, I inted to make multiple plots using for loop. The question is how can I adjust the left hand side margin of the plot according to the names.arg argument in barplot2. In every plot I have different annotations in the y axis and they vary in length. Now when I have fixed margins opar

[R] Additions to xyplot (lattice)? - legend, ticks, axis label size, text

I have created an xyplot of a time series with the following code... win.graph(width = 10, height = 7) panel1 = function(x, y) { panel.loess(x, y, lwd=2.5, span=0.5, col=gray) panel.xyplot(x, y, pch=19, col=blue, cex=1.25) } xyplot(oneplusdensity ~ year, data=figdata,

Re: [R] Q: how to interrupt long calculation?

On linux and solaris, I have had some luck doing a CTRL-Z (STOP signal) and the resuming later... some libraries have issues (I don't recall which) but...usually pretty good. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of D. R. Evans Sent: Thursday,

Re: [R] Q: how to interrupt long calculation?

Prof Brian Ripley said the following at 08/30/2007 11:00 AM : On Thu, 30 Aug 2007, D. R. Evans wrote: Paul Smith said the following at 08/29/2007 04:32 PM : The instance of R running will be immediately killed and then you can start R again. But then I would lose all the work. There must

Re: [R] Additions to xyplot (lattice)? - legend, ticks, axis label size, text

On 8/30/07, Dave Hewitt [EMAIL PROTECTED] wrote: I have created an xyplot of a time series with the following code... win.graph(width = 10, height = 7) panel1 = function(x, y) { panel.loess(x, y, lwd=2.5, span=0.5, col=gray) panel.xyplot(x, y, pch=19, col=blue, cex=1.25)

Re: [R] Month end calculations

And one more yearmon solution. Here z is a zoo series as before: tt - time(z) aggregate(z, ave(tt, as.yearmon(tt), FUN = max), tail, 1) On 8/30/07, Gabor Grothendieck [EMAIL PROTECTED] wrote: The last line is wrong (see below for correction): On 8/30/07, Gabor Grothendieck [EMAIL PROTECTED]

Re: [R] Q: how to interrupt long calculation?

On Thu, 30 Aug 2007, D. R. Evans wrote: Prof Brian Ripley said the following at 08/30/2007 11:00 AM : On Thu, 30 Aug 2007, D. R. Evans wrote: Paul Smith said the following at 08/29/2007 04:32 PM : The instance of R running will be immediately killed and then you can start R again. But

  1   2   >