Re: [R] The KJV

On 07-Feb-10 01:06:40, Ben Bolker wrote: Jim Lemon jim at bitwrit.com.au writes: On 02/06/2010 06:57 PM, Charlotte Maia wrote: Hey all, Does anyone know if there are any R packages with a copy of the KJV? I'm guessing the answer is no... So the next question, and the more

Re: [R] embedFonts with pdf files and Windows 7

Do your systems actually have the fonts you are trying to embed? I doubt it: Helvetica is a commercial font, and most likely the Linux system is embedding a substitute. It would be better to do pdf(test.pdf, family=NimbusSan, useDingbats=FALSE) plot(matrix(rnorm(200),nc=2)) dev.off() myCall -

Re: [R] Is there an R implementation for the Barnard's exact test (a substitute for fisher.test) ?

Hello all, After almost half a year, I received a friendly e-mail from Peter Calhoun, sharing his R implementation of Barnard's exact test. With his permission, I posted his code here:

Re: [R] The KJV

On Sun, Feb 7, 2010 at 8:28 AM, Ted Harding ted.hard...@manchester.ac.uk wrote: Delightful! And fascinating in the detail too.  length(tt)  # [1] 5078 with slight changes like:  barplot(rev(tt[1:50]),horiz=TRUE,las=1,cex.names=0.6,log=x)  # ...  

Re: [R] The KJV

On 07-Feb-10 12:49:23, Barry Rowlingson wrote: On Sun, Feb 7, 2010 at 8:28 AM, Ted Harding ted.hard...@manchester.ac.uk wrote: Delightful! And fascinating in the detail too. _length(tt) _# [1] 5078 with slight changes like: _barplot(rev(tt[1:50]),horiz=TRUE,las=1,cex.names=0.6,log=x)

[R] predicting with stl() decomposition

Hi mailinglist members, I’m actually working on a time series prediction and my current approach is to decompose the series first into a trend, a seasonal component and a remainder. Therefore I’m using the stl() function. But I’m wondering how to get the single components in order to predict

[R] optimized R-selection and R-replacement inside a matrix need, strings coerced to factors

Hello, I need to modify some huge arrays (2000 individuals x 50 000 variables). To format the data, I think I should benefit from optimized R-selection and R-replacement inside a matrix and prohibite a naive use of loops. Thank you in advance for providing information about the following

Re: [R] Posting an 'S4-creating Package Problem'...

On 02/06/2010 03:39 PM, Daniel Kosztyla wrote: Hello R-Team, May you help me to post a 'S4-creating Package Problem'? Thanks already now for supporting. The problem sounds like: Hello R forum, while compiling my R-package these 'Warnings' occur: ... Warnung in

[R] convert R plots into annotated web-graphics

Dear all, I would like to make a large scatter plot created with R available as an interactive web graphic, in combination with additional text-annotations for each data point in the plot. The idea is to present the text-annotations in an HTML-table and inter-link the data points in the plot

Re: [R] R-Help

On Sat, Feb 6, 2010 at 2:46 PM, David Winsemius dwinsem...@comcast.net wrote: On Feb 6, 2010, at 3:29 PM, Ravi Ramaswamy wrote: Hi - I am not familiar with R.  Could I ask you a quick question? When I read a file like this, I get an error.  Not sure what I am doing wrong.  I use a MAC.  How

Re: [R] Non-linear regression

It appears my suspicions about this being homework were unfounded. Given the additional problems with excess zeroes, you may want to examine the extremely informative material on analysis of such problems written by Zeileis, Kleiber and Jackman: (easily found in case you have misplaced it,

Re: [R] metafor package: effect sizes are not fully independent

Dear Gang, Here are just some general thoughts. Wolfgang Viechtbauer will be a better position to answer questions related to metafor. For multivariate effect sizes, we first have to estimate the asymptotic sampling covariance matrix among the effect sizes. Formulas for some common effect sizes

Re: [R] predicting with stl() decomposition

Hi, yes that error name is indeed kind of weird. But I think it’s thrown due to the missing robustness of the estimation since every weight is one and hence the fit is likely to be influenced by outliers in the provided data which should be just an example. But do you have an idea to extract

[R] Reading hierarchical data

I would like to read the following hierarchical data set. There is a family record followed by one or more personal records. If col. 7 is 1 it is a family record. If it is 2 it is a personal record. The family record is formatted as follows: col. 1-5 family id col. 71 col. 9

Re: [R] predicting with stl() decomposition

Hi Dennis, I’ve already found that matrix, but I want to predict the time series, for example with predict.loess() on every component. But actually I’m unable to extract the plotted series’ (trend and seasonal) as a loess object. This representation is what I’m looking for. At the moment I

Re: [R] Reading hierarchical data

Will this do it for you: input - readLines(textConnection(06470 1 1 + 1 232 0 + 2 230 1 + 07470 1 0 + 1 240 1 + 08470 1 0 + 1 227 0 + 09470 1 0 + 1 213 1 + 2 222 0 + 3 224 1 + 10470 1 1 + 1 220 0 + 2 211 1 + 11470 1 0 + 1 217 0 + 2 210 1 + 3 226

[R] conditioned xyplot, many y variables

The example below creates parallel time-series plots of three different y variables conditioned by a dichotomous factor. In the graphical layout, • Each y variable inhabits its own row and is plotted on its own distinct scale. • Each level of the factor has its

[R] Why does aggregate fail?

I am trying to get hourly totals, given 15-minute bins. s = seq(0, 95, 1) s = floor(s/4) # 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3 4 . . . s [1] 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 5 5 6 [26] 6 6 6 7 7 7 7 8 8 8 8 9 9 9 9 10 10 10 10 11

Re: [R] Why does aggregate fail?

On Feb 7, 2010, at 1:08 PM, James Rome wrote: I am trying to get hourly totals, given 15-minute bins. s = seq(0, 95, 1) s = floor(s/4) # 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3 4 . . . s [1] 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 5 5 6 [26] 6 6 6

Re: [R] Why does aggregate fail?

dput(d) structure(list(`0` = 0, `1` = 1, `2` = 0, `3` = 0, `4` = 0, `5` = 0, `6` = 0, `7` = 0, `8` = 0, `9` = 0, `10` = 0, `11` = 0, `12` = 0, `13` = 0, `14` = 0, `15` = 0, `16` = 0, `17` = 2, `18` = 0, `19` = 0, `20` = 0, `21` = 0, `22` = 0, `23` = 0, `24` = 0, `25` = 0, `26` =

Re: [R] Why does aggregate fail?

You have a dataframe with 96 columns and a single row named Sunday. My guess is that was not your intent. How did d come to exist? -- David. On Feb 7, 2010, at 1:29 PM, James Rome wrote: dput(d) structure(list(`0` = 0, `1` = 1, `2` = 0, `3` = 0, `4` = 0, `5` = 0, `6` = 0, `7` = 0, `8` =

Re: [R] (Another) Bates fortune?

On Fri, 5 Feb 2010, Peter Ehlers wrote: I vote to 'fortunize' Doug Bates on Hierarchical data sets: which software to use? The widespread use of spreadsheets or SPSS data sets or SAS data sets which encourage the single table with a gargantuan number of columns, most of which are missing data

Re: [R] Why does aggregate fail?

On Feb 7, 2010, at 1:32 PM, David Winsemius wrote: You have a dataframe with 96 columns and a single row named Sunday. My guess is that was not your intent. How did d come to exist? But to answer your question: apply(d, 1, function(z) aggregate(z, by=list(s), FUN=sum) ) $Sunday

Re: [R] Why does aggregate fail?

On 2/7/2010 1:32 PM, David Winsemius wrote:You have a dataframe with 96 columns and a single row named Sunday. My guess is that was not your intent. How did d come to exist? I was trying to make a simpler example. The actual code uses a data frame maxrdf: dput(maxrdf) structure(list(`0` = c(0,

Re: [R] Why does aggregate fail?

On 2/7/2010 1:35 PM, David Winsemius wrote:But to answer your question: apply(d, 1, function(z) aggregate(z, by=list(s), FUN=sum) ) David, That works, but I do not understand why I could not use aggregate directly. And the answer comes out as a list, which thus far baffles me. How do I get the

Re: [R] conditioned xyplot, many y variables

On Sun, Feb 7, 2010 at 9:32 AM, Jacob Wegelin jacobwege...@fastmail.fm wrote: The example below creates parallel time-series plots of three different y variables conditioned by a dichotomous factor. In the graphical layout,        •       Each y variable inhabits its own row and is plotted on

Re: [R] Why does aggregate fail?

This works. But I wish I could write it without a lot of trial and error. :-( ha = matrix(nrow=7, ncol=24) colnames(ha) = as.character(c(0:23)) rownames(ha) = rownames(maxrdf) m = as.matrix(maxrdf) for(j in 1:7) { x = aggregate(m[j,], by=list(s), FUN=sum) ha[j,] =

Re: [R] Why does aggregate fail?

On Feb 7, 2010, at 1:57 PM, James Rome wrote: On 2/7/2010 1:35 PM, David Winsemius wrote:But to answer your question: apply(d, 1, function(z) aggregate(z, by=list(s), FUN=sum) ) David, That works, but I do not understand why I could not use aggregate directly. And the answer comes out

Re: [R] convert R plots into annotated web-graphics

On Sun, Feb 7, 2010 at 2:35 PM, Rainer Tischler rainer_...@yahoo.de wrote: Dear all, I would like to make a large scatter plot created with R available as an interactive web graphic, in combination with additional text-annotations for each data point in the plot. The idea is to present the

Re: [R] convert R plots into annotated web-graphics

On Sun, Feb 7, 2010 at 2:35 PM, Rainer Tischler rainer_...@yahoo.de wrote: If you have alternative ideas for interlinking tabular annotations with plotted data points, I would appreciate any recommendation/suggestion. (I work with R 2.8.1 on different 32-bit PCs with both Linux and Windows

Re: [R] conditioned xyplot, many y variables

On Sun, Feb 7, 2010 at 11:32 AM, Jacob Wegelin jacobwege...@fastmail.fm wrote: The example below creates parallel time-series plots of three different y variables conditioned by a dichotomous factor. In the graphical layout,        •       Each y variable inhabits its own row and is plotted

[R] x-axis plot problem

Hi all, I tried to have plot of many vector in one plot and i have got a nice plot but i have problem with x-axis. I want to have month and year only(Jul.07 means July 2007) in x-axis without appearing other number behaind it. I would appercit any help. The R code:

Re: [R] x-axis plot problem

On Feb 7, 2010, at 3:22 PM, abotaha wrote: Hi all, I tried to have plot of many vector in one plot and i have got a nice plot but i have problem with x-axis. I want to have month and year only(Jul.07 means July 2007) in x-axis without appearing other number behaind it. I'm going to

Re: [R] x-axis plot problem

I think you just need to set axes=FALSE in your call to matplot(). You'll then need to add the y-axis manually --- do axis(2) in addition to your call which draws the x axis. You'll also need to do box() if you want a box around your graph. cheers, Rolf Turner P. S.

Re: [R] x-axis plot problem

Hi abotaha, Modify your matplot() call as matplot(model, pch = c(1,22,17,16), type = o,lty=c(2,2,2,5), col =c(gray10, gray10,gray10,gray10),xlab=Month-Year,ylab=Zinth, xaxt = n, yaxs = i,main=Model Output) and then add axis(1, 1:6, time) HTH, Jorge On Sun, Feb 7, 2010 at 3:22 PM,

[R] Noval numbers

Dear everybody, How can I transform numbers to a positional system with the base of, e.g., nine, and do further operations with them? Thank you in advance Yours, sincerely Mag. Ferri Leberl __ R-help@r-project.org mailing list

Re: [R] Noval numbers

PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. It would be useful if you could at least provide some examples of what you want to do. There are various ways of converting numbers back and forth. Are

[R] using a variable name stored in another variable?

Hello, I'm trying to figure out how to create a data object, and then save it with a user-defined name that is input as a command line argument. I know how to create the object and assign it the new name, however, I can't figure out how to refer to the new name for a future operation such as

[R] contour persp

I have this data set that both x y are ordered vectors of length 600 700 respectively; z is a 600 by 700 matrix whose entry z[i,j] is either a missing value (indicated by 'NaN') or a real number between 0 and 1. The contour function contour(x,y,z) gives me a blank picture. I guess the

Re: [R] Noval numbers

On 07/02/2010 4:25 PM, Mag. Ferri Leberl wrote: Dear everybody, How can I transform numbers to a positional system with the base of, e.g., nine, and do further operations with them? I don't understand what you want. Decimal, noval or binary are just ways to represent numbers as strings of

[R] mboost: Interpreting coefficients from glmboost if center=TRUE

I'm running R 2.10.1 with mboost 2.0 in order to build predictive models . I am performing prediction on a binomial outcome, using a linear function (glmboost). However, I am running into some confusion regarding centering. (I am not aware of an mboost-specific mailing list, so if the main R list

Re: [R] using a variable name stored in another variable?

Chris Seidel wrote: Hello, I'm trying to figure out how to create a data object, and then save it with a user-defined name that is input as a command line argument. I know how to create the object and assign it the new name, however, I can't figure out how to refer to the new name for a

[R] Data views (Re: (Another) Bates fortune?)

Note : this post has been motivated more by the hierarchical data subject than the aside joke of Douglas Bates, but might be of interest to its respondents. Le vendredi 05 février 2010 à 21:56 +0100, Peter Dalgaard a écrit : Peter Ehlers wrote: I vote to 'fortunize' Doug Bates on

[R] specifying colors in a heatmap/image -like plot

Hi, I have searched for a solution but I failed to find an answer. I am hoping you may be able to help me. I have a data set where I have observations for a number of units (n =~40) over a period of time (t =~100) and I have a variable (Z) that codes a categorical variable for each

Re: [R] Interactively editing point labels in a graph

The built-in R graphics system was not designed for interactivity -- there is no [feasible] way to detect the data point coordinates in a base graphics plot. The playwith package tries to figure out the coordinates from the data objects given in the call: this works for simple scatterplots etc,

[R] Out-of-sample prediction with VAR

Good day, I'm using a VAR model to forecast sales with some extra variables (google trends data). I have divided my dataset into a trainingset (weekly sales + vars in 2006 and 2007) and a holdout set (2008). It is unclear to me how I should predict the out-of-sample data, because using the

Re: [R] using a variable name stored in another variable?

Hi Charlie, get() will return the contents (value) of a variable. But what I want is to save the named object. Something like save(get(myobjectname), ...) doesn't work. In the environment, is that object of interest, and a variable which holds the name of the object of interest. If you don't

Re: [R] specifying colors in a heatmap/image -like plot

On Feb 7, 2010, at 4:57 PM, kerimcan wrote: Hi, I have searched for a solution but I failed to find an answer. I am hoping you may be able to help me. I have a data set where I have observations for a number of units (n =~40) over a period of time (t =~100) and I have a variable (Z)

Re: [R] using a variable name stored in another variable?

On 07/02/2010 6:05 PM, Chris Seidel wrote: Hi Charlie, get() will return the contents (value) of a variable. But what I want is to save the named object. Something like save(get(myobjectname), ...) doesn't work. I think you want save(list=myobjectname, file= ...) assuming that the object

Re: [R] using a variable name stored in another variable?

Tena koe Chris Does the following help? dfName - 'myDf' save(dfName, file='test1') save('dfName', file='test2') save('myDf', file='test3') save(myDf, file='test4') Peter Alspach -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf

Re: [R] Why does smoothScatter clip when xlim and ylim increased?

On Sat, Feb 6, 2010 at 6:15 AM, Duncan Murdoch murd...@stats.uwo.ca wrote: On 06/02/2010 7:51 AM, Jennifer Lyon wrote: Hi: Is there a way to get smoothScatter to not clip when I increase the xlim and ylim parameters? Consider the following example: set.seed(17) x1-rnorm(100)

Re: [R] Reading hierarchical data

Try this. It uses input defined in Jim's post and defines the rectype of each row (1 or 2). It then reads the rectype 1 records into DF1 using read.fwf and the rectype 2 records into DF2 also using read.fwf. ix is defined to have one component per personal record giving the row number in DF1 of

Re: [R] mboost: Interpreting coefficients from glmboost if center=TRUE

On Feb 7, 2010, at 5:03 PM, Kyle Werner wrote: I'm running R 2.10.1 with mboost 2.0 in order to build predictive models . I am performing prediction on a binomial outcome, using a linear function (glmboost). However, I am running into some confusion regarding centering. (I am not aware of an

Re: [R] Interactively editing point labels in a graph

Create your plot and save it in wmf format, e.g. DF - as.data.frame(state.x77) plot(Income ~ log(Population), DF, pch = 20) with(DF, text(log(Population), Income, rownames(state.x77), cex = 0.5, pos = 4)) savePlot(states.wmf) Then insert it into Microsoft Word, right click the image, choose Edit

Re: [R] Reading hierarchical data

Here is a further simplification. We use the colClasses= argument with NULL for the columns we do not want so we do not have to later remove those columns. # record type (1 or 2) rectype - substr(input, 7, 7) # read in record type 1 input1 - input[rectype == 1] DF1 -

[R] dataframe question

Folks: Good day. Please see the code below. three_wk_out is a dataframe with columns wk1 through wk209. I want to change the format of the columns. I am trying the code below but it does not work. I need $week in the for loop interpreted as wk1, wk2, etc. Could you please help? Thanks. Satish

Re: [R] dataframe question

Tena koe Satish Try using three_wk_out[,week] - as.numeric(tree_wk_out[,week]) HTH Peter Alspach -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Vadlamani, Satish {FLNA} Sent: Monday, 8 February 2010 1:51 p.m. To:

Re: [R] 3D plot of following data

Hi Jim Lemon wrote: On 02/02/2010 11:01 PM, walter.dju...@chello.at wrote: Hello R-experts, I am having difficulties with 3D plotting (i.e. the evolution of various forward curves through time). I have two comma seperated files both ordered by date (in the first column) one containing

Re: [R] dataframe question

On Feb 7, 2010, at 7:51 PM, Vadlamani, Satish {FLNA} wrote: Folks: Good day. Please see the code below. three_wk_out is a dataframe with columns wk1 through wk209. I want to change the format of the columns. I am trying the code below but it does not work. I need $week in the for loop

Re: [R] metafor package: effect sizes are not fully independent

Dear Gang, It seems that it is possible to use a univariate meta-analysis to handle your multivariate effect sizes. If you want to calculate a weighted average first, Hedges and Olkin (1985) has discussed this approach. Hedges, L. V., Olkin, I. (1985). Statistical methods for meta-analysis.

Re: [R] Noval numbers

The attached file gives functions to go both directions. I have used it in class for many years. This is very useful when studying machine representations of numbers, for understanding mixed-radix number systems, for example time (days, hours, minutes, seconds) or British money (pounds,

Re: [R] dataframe question

On Feb 7, 2010, at 8:14 PM, David Winsemius wrote: On Feb 7, 2010, at 7:51 PM, Vadlamani, Satish {FLNA} wrote: Folks: Good day. Please see the code below. three_wk_out is a dataframe with columns wk1 through wk209. I want to change the format of the columns. I am trying the code below

[R] split plot with aov

I have a factor SAMPLES which is at the lowest level (within) of a split plot anova model but this factor also appears in the ANOVA table at the block level. This happens for unbalanced responses but not for balanced responses. I would be grateful for an explanation of this. The block error

[R] Help with apply()

I have a 2 column data.frame: d[1:5,] a b 180015 C 280016 B 380023 C 480062 B 580069 B I want to apply a function across each row: for(i in 1:nrow(d)) { +myFun(con, d[i,]$a, d[i,]$b) + } How do I do this using apply()? I'm unsure how to tell

Re: [R] split plot with aov

The dummy variables for the factors in balanced designs are orthogonal. The treatment dummy variables are not orthogonal to the block dummy variables for unbalanced designs. That is essentially what the term balanced means. __ R-help@r-project.org

Re: [R] mboost: Interpreting coefficients from glmboost if center=TRUE

Thanks for your reply. In fact, I do use the predict method for model assessment, and it shows that centering leads to a substantial improvement using even the bluntest of assessments of 'goodness' (i.e., binary categorization accuracy). So I agree that the package authors must have internal tools

Re: [R] Help with apply()

On Feb 7, 2010, at 8:26 PM, Nathan S. Watson-Haigh wrote: I have a 2 column data.frame: d[1:5,] a b 180015 C 280016 B 380023 C 480062 B 580069 B I want to apply a function across each row: for(i in 1:nrow(d)) { +myFun(con, d[i,]$a, d[i,]$b) + }

Re: [R] dataframe question

David: Thanks for the idea. Both the one that you suggested and the one that Bill Venables suggested are very good. Unfortunately, this statement is creating out of memory issues like below (system limitations). When I had padded white space before the number, read.csv.sql is correctly

[R] Contributed packages

Folks: If you wanted to find out about what are the contributed packages and classify them, how would you go about it? For someone new like me, I would like to know what the possibilities are. When I click on install packages on my Windows version of R, it gives me a list but it is hard to

Re: [R] dataframe question

On Feb 7, 2010, at 11:15 PM, Vadlamani, Satish {FLNA} wrote: David: Thanks for the idea. Both the one that you suggested and the one that Bill Venables suggested are very good. Unfortunately, this statement is creating out of memory issues like below (system limitations). When I had

Re: [R] Contributed packages

On Feb 7, 2010, at 11:27 PM, Vadlamani, Satish {FLNA} wrote: Folks: If you wanted to find out about what are the contributed packages and classify them, how would you go about it? For someone new like me, I would like to know what the possibilities are. When I click on install packages

[R] ggplot2 stacked line plot

Hi all, I have been hunting around for hours trying to figure out how to generate a stacked line chart using ggplot2. This type of chart can be generated in excel 2007 by selecting: Chart type Line Stacked line. I can generate a stacked area chart using the following code: p - ggplot2(~,

[R] problem with Tinn-R

Hi, I install Tinn-R 2.3.4.4 and when I want to execute the calculation, it gives me this error:   The preferred Rterm not defined.   Thank you so much for any help given. [[alternative HTML version deleted]] __ R-help@r-project.org

Re: [R] ggplot2 stacked line plot

Hi, On Sun, Feb 7, 2010 at 11:40 PM, Liam Blanckenberg liam.blanckenb...@gmail.com wrote: Hi all, I have been hunting around for hours trying to figure out how to generate a stacked line chart using ggplot2. This type of chart can be generated in excel 2007 by selecting: Chart type Line

Re: [R] specifying colors in a heatmap/image -like plot

On 02/08/2010 08:57 AM, kerimcan wrote: Hi, I have searched for a solution but I failed to find an answer. I am hoping you may be able to help me. I have a data set where I have observations for a number of units (n =~40) over a period of time (t =~100) and I have a variable (Z) that codes a

Re: [R] Help with apply()

On 02/08/2010 12:26 PM, Nathan S. Watson-Haigh wrote: I have a 2 column data.frame: d[1:5,] a b 1 80015 C 2 80016 B 3 80023 C 4 80062 B 5 80069 B I want to apply a function across each row: for(i in 1:nrow(d)) { + myFun(con, d[i,]$a, d[i,]$b) + } How do I do this using apply()? I'm

Re: [R] ggplot2 stacked line plot

On 02/08/2010 03:40 PM, Liam Blanckenberg wrote: Hi all, I have been hunting around for hours trying to figure out how to generate a stacked line chart using ggplot2. This type of chart can be generated in excel 2007 by selecting: Chart type Line Stacked line. I can generate a stacked area

Re: [R] problem with Tinn-R

Roslina Zakaria wrote: I install Tinn-R 2.3.4.4 and when I want to execute the calculation, it gives me this error:   The preferred Rterm not defined.   Set the path to Rtermn in Options/Application/R/Path Dieter -- View this message in context: