[R] Find the prediction or the fitted values for an lm model

2013-11-27 Thread jpm miao
Hi, I would like to fit my data with a 4th order polynomial. Now I have only 5 data point, I should have a polynomial that exactly pass the five point Then I would like to compute the "fitted" or "predict" value with a relatively large x dataset. How can I do it? BTW, I thought the mode

[R] importing many csv files into separate matrices

2013-11-27 Thread yetik serbest
Hi Everyone,   I am trying to import many CSV files to their own matrices. Example, alaska_93.csv to alaska. When I execute the following, for each csv.file separately it is successful.   singleCSVFile2Matrix <- function(x,path) {  assign(gsub(pattern=".csv",x,replacement=""),read.csv(paste(path,

[R] heat map bin locations

2013-11-27 Thread Heath Blackmon
I am trying to construct a het map type of plot with rectangular bins colored to indicate the number of records for a given value within different categories. I am getting close to what I want with ggplot2 and stat_bin2d. However if you run the code you will se below that the placement of the bin

[R] error message

2013-11-27 Thread Mosiuoa Bereng
Dear Sir/Madam I have a problem with r. I was still using it on my laptop and it then later refused to open. When I open a script, being a new script or the previous one that I have saved, it displays a message saying "fatal error: unexpected error reading console input". I then uninstalled th

[R] if, apply, ifelse

2013-11-27 Thread Andrea Lamont
Hello: This seems like an obvious question, but I am having trouble answering it. I am new to R, so I apologize if its too simple to be posting. I have searched for solutions to no avail. I have data that I am trying to set up for further analysis ("training data"). What I need is 12 groups based

[R] Requirement for Ruby on rails [REQ:10217055]

2013-11-27 Thread Revanth kumar Gudimalla
[1]Click here to unsubscribe if you no longer wish to receive our emails Dear Recruiter, Here is our Direct client requirement which can be filled immediately. Kindly respond to this requirement with your consultant resume, contact and current location info to speed u

Re: [R] Y-axis label not plotting

2013-11-27 Thread Jim Lemon
On 11/28/2013 02:10 AM, Luke M wrote: Thanks for the follow-up. Per your suggestion, I was able to plot the data for the last 30 days with no problem. However... I am still confused about the y-axis labels (values, title, etc) not plotting. If the scripting works as it should, what could cause t

[R] Tinn-R: source code on GitHub

2013-11-27 Thread Jose Claudio Faria
Dear list, The source code of Tinn-R editor is available on GitHub: https://github.com/jcfaria/Tinn-R Tinn-R is free, simple but efficient replacement for the basic code editor provided by Rgui. The project is coordinate by José Cláudio Faria/UESC/DCET. All users are welcome to make it better. L

Re: [R] Etimating time to run an analysis?

2013-11-27 Thread Bert Gunter
I would say that if the OP even contemplated this, it strongly suggests that she needs to consult a local statistician for help. Cheers, Bert On Wed, Nov 27, 2013 at 1:14 PM, Ben Bolker wrote: > Erika Barthelmess stlawu.edu> writes: > >> Hi everyone, >> >> I'm new to this list and have searched

Re: [R] xtable: custom row.names, move caption to top

2013-11-27 Thread Michael Friendly
Thanks so much, Rainer. Your detailed example has taught me a lot of what I need to use xtable more productively, in particular the options for the print() method. -Michael On 11/27/2013 12:20 PM, Rainer Schuermann wrote: UPDATE: Now including the LaTeX math formatting saxtab <- t( as.data.fra

Re: [R] Etimating time to run an analysis?

2013-11-27 Thread Ben Bolker
Erika Barthelmess stlawu.edu> writes: > Hi everyone, > > I'm new to this list and have searched R help prior for an answer > to this question, without luck. If I'm > posting in error, please forgive. > > I'm thinking about using package MuMIn to do multimodel inference > with logistic regres

Re: [R] R- package for Parametric Survival Analysis with Left-censored data

2013-11-27 Thread Göran Broström
On 11/20/2013 08:17 AM, peter dalgaard wrote: On 20 Nov 2013, at 04:15 , David Winsemius wrote: On Nov 19, 2013, at 5:30 PM, Vinod Mishra wrote: Dear All, I am new to R. Can someone please direct me to an R package using which I can estimate a Parametric Survival Analysis model with Lef

Re: [R] Space in label name

2013-11-27 Thread jim holtman
try check.names = FALSE in the read.table to see if this helps. Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it. On Wed, Nov 27, 2013 at 7:21 AM, Gitte Brinch Andersen < gitt...@hum-gen.au.dk> wrote: > Hi

Re: [R] cut2 not binning interval endpoints correctly

2013-11-27 Thread William Dunlap
You can look at the source code of Hmisc::cut2() to see what is going on -- it does a lot more than calling cut() with different default arguments. Another approach to debugging this is to use trace() to see what cut2() passes down to the default cut method: > trace(cut.default, quote(cat(" x

[R] ifelse, apply, if

2013-11-27 Thread Andrea Lamont
Hello: This seems like an obvious question, but I am having trouble answering it. I am new to R, so I apologize if its too simple to be posting. I have searched for solutions to no avail. I have data that I am trying to set up for further analysis ("training data"). What I need is 12 groups based

Re: [R] Y-axis label not plotting

2013-11-27 Thread Luke M
Thanks for the follow-up. Per your suggestion, I was able to plot the data for the last 30 days with no problem. However... I am still confused about the y-axis labels (values, title, etc) not plotting. If the scripting works as it should, what could cause the labels not to show up? A while ago I

Re: [R] xtable: custom row.names, move caption to top

2013-11-27 Thread Rainer Schuermann
UPDATE: Now including the LaTeX math formatting saxtab <- t( as.data.frame( addmargins( Saxony ) ) ) rownames( saxtab ) <- c( "Males ($k$)", "Families ($n_k$)" ) saxtab <- xtable( saxtab, digits = 0, caption = "Number of male children in 6115 Saxony families of size 12", align = "l|r

Re: [R] cut2 not binning interval endpoints correctly

2013-11-27 Thread S Ellison
> -Original Message- >jim holtman > You need to look at the full accuracy of the number representation: Um... I think I did. But I'm not sure you did print(..., digits=20) has used different numbers of digits for your two print()s, probably because print() decided it needed more di

Re: [R] Coding systems.

2013-11-27 Thread gerald . jean
Hello, as Jan pointed out the problem is with the encoding in which R saves the fucntion. If I set this encoding to "UTF-8" in source everything is fine. If I go either in my .bash_profile or my .Renviron file and set all LOCALE variables to "fr_CA.UTF8" it should do the job, and to a certain p

Re: [R] xtable: custom row.names, move caption to top

2013-11-27 Thread Rainer Schuermann
You get the cation to the top of the table with print( saxtab, caption.placement = "top" ) Formatting the table the way you want can be done like this - I did not manage to carry the LaTeX math formatting for the row names over ($k$ and $n_k$)but the rest should be very much what you want: saxt

Re: [R] Conditional error bars

2013-11-27 Thread William Dunlap
Change your two lines > arrows(xvals, mean, xvals, CI.H, angle=90,length=length) > arrows(xvals, mean, xvals, CI.L, angle=90,length=length) to the one line arrows(xvals, mean, xvals, mean + sign(mean) * se, angle=90, length=length) (I would also use a scatter plot instead of a barp

[R] Conditional error bars

2013-11-27 Thread raz
How can I condition any error bar function that use the arrows() function, such as 'CI.plot' (see example below) or 'error.bars', to draw only upper error bar (upper CI) if the bar value (mean) is positive and the lower error bar (lower CI) if bar value is negative? CI.plot <- function(mean, se,le

Re: [R] Buse's GLS R2

2013-11-27 Thread Sheri O'Connor
Thanks very much Pascal. That is a very useful package! On Wed, Nov 27, 2013 at 1:26 AM, Pascal Oettli wrote: > Hello, > > Using 'sos': > > library(sos) > findFn('buse') > > Hope this helps, > Pascal > > On 27 November 2013 14:55, Sheri O'Connor wrote: >> I was wondering if anyone knew of a pack

[R] Etimating time to run an analysis?

2013-11-27 Thread Erika Barthelmess
Hi everyone, I'm new to this list and have searched R help prior for an answer to this question, without luck. If I'm posting in error, please forgive. I'm thinking about using package MuMIn to do multimodel inference with logistic regression. I have many (25) possible predictors and am curio

[R] xtable: custom row.names, move caption to top

2013-11-27 Thread Michael Friendly
With xtable, I'm producing one-way tables from table objects in horizontal form as shown below. I'd like to change the labels used for the rows and move the caption to the top of the table, as is typically standard for tables. I can hand-edit, but would prefer to do it in code. data(Saxony, p

[R] Space in label name

2013-11-27 Thread Gitte Brinch Andersen
Hi I am doing a cluster analysis and I have some troubles with the label names. In the text file I have my data in, the first row contains the names of the samples. They are called "FFPE Tumor 2", "Fresh Frozen Tumor 2" etc. But when the cluster is made the label names contains periods instead

Re: [R] Should there be an R-beginners list?

2013-11-27 Thread Joris Meys
StackOverflow has certainly its merits, although I miss a bit the good ol' Oxford sarcasm gems you find in this list. This said : Beginner's list. Bad, bad idea. First rule in my classes is: RTFI (Read The Fucking Internetzz). Anybody using R should be able to do a basic Google search. A beginner'

Re: [R] Tinn-R user guide (latex sources) available on GitHub

2013-11-27 Thread Joris Meys
Thank you for doing this! On Wed, Nov 27, 2013 at 11:22 AM, Jose Claudio Faria < joseclaudio.fa...@gmail.com> wrote: > Dears user, > > The Tinn-R User Guide is completely written in LaTeX and the idea > behind this to be available on GitHub is that it has contributions > from multiple users. > >

Re: [R] Y-axis label not plotting

2013-11-27 Thread Jim Lemon
On 11/27/2013 05:32 PM, Luke M wrote: Dear R community I am trying to make an XY plot that shows temperature (y axis) as a function of time (x axis) but I am having some problems. When I use the code shown below: 1. my plot does not show any of the y-axis labels even though there is plenty of w

[R] Tinn-R user guide (latex sources) available on GitHub

2013-11-27 Thread Jose Claudio Faria
Dears user, The Tinn-R User Guide is completely written in LaTeX and the idea behind this to be available on GitHub is that it has contributions from multiple users. If you find something that you would like to include or impruve: please, fell free to make it better. This User Guide have been de

[R] MANOVA Question

2013-11-27 Thread Chris Linton
I'm running a MANOVA on survey data. The dependent variables are the ranking of each DV. So the respondent was asked to rank 7 different subject lines, each subject line is a DV with what the respondent ranked as datum. I had the code running yesterday, had a computer crash, then had to rewrite

[R] Y-axis label not plotting

2013-11-27 Thread Luke M
Dear R community I am trying to make an XY plot that shows temperature (y axis) as a function of time (x axis) but I am having some problems. When I use the code shown below: 1. my plot does not show any of the y-axis labels even though there is plenty of white space there (i.e.: I am getting no

[R] GADM Data Download

2013-11-27 Thread Lorenzo Isella
Dear All, Please consider the snippet at the end of the email. I often download some maps (in the R format) from http://www.gadm.org/ However, when I run (typically more than once) a variation of the script below (based on http://bit.ly/1b3W0Aa ), I often get Error in load(url(paste("http:/

Re: [R] Should there be an R-beginners list?

2013-11-27 Thread Rolf Turner
On 11/25/13 09:04, Rich Shepard wrote: On Sun, 24 Nov 2013, Yihui Xie wrote: Mailing lists are good for a smaller group of people, and especially good when more focused on discussions on development (including bug reports). The better place for questions is a web forum. I disagree. Mail lis