Re: [R] Chron object in time series plot

2010-10-20 Thread Robert Baer
It seems David and Phil have given you all the pieces and you are not listening. Try: h = sprintf('%06d', seq(07,18, by=1) ) # Note that h is text myseries=sample(1:12, 12) #fake data to plot tim=strptime(h,'%H%M%S') library(zoo) test=zoo(myseries, tim) plot(test) Rob

Re: [R] Plot creates a straigth line

2010-10-20 Thread Robert Baer
od-c(10, 8, 6,4,2,1, 10.5,7.8,6.4,3.8,2.1,0.95) cyto_conc=2650 # Highest cytokine concentration user defined cyto_std_conc -c(cyto_conc) for (i in 1:5) { cyto_conc = cyto_conc /3 cyto_std_conc -c(cyto_std_conc ,cyto_conc) } cyto_std_conc-log2(rep(cyto_std_conc,2))

Re: [R] best predictive model for mixed catagorical/continuous variables

2010-10-25 Thread Robert Baer
I would expect that there will be interactions between variables eg. if the vegetation is grassland then the vegetation height variable will mediate the interaction, if the vegetation is arable then crop type will be more significant. Would it be possible to use GLM or GAM models for this type

Re: [R] JGR install problem

2010-11-12 Thread Robert Baer
You install JGR and its dependencies as you would any other package from the RGUI. Make sure you have a recent java installed and functional. You may find jselect.exe on the RForge site useful for looking at this. When I did the install, things seemed a little flaky, but eventually by

Re: [R] Searching dataset

2010-11-30 Thread Robert Baer
?RSiteSearch -- From: Stephen Liu sati...@yahoo.com Sent: Tuesday, November 23, 2010 9:09 AM To: r-help@r-project.org Subject: [R] Searching dataset Hi folks, help.search(item_name) displays a list of dataset containing the item_name. If I

Re: [R] Learning ANOVA

2010-08-16 Thread Robert Baer
modelfn(count ~ spray, data=test01) Error: could not find function modelfn What will be +var3 ? I think the general advice applied to your specific situation was meant to suggest using the construct: InsectSprays.aov - aov(count ~ spray, data=InsectSprays) rather than: InsectSprays.aov -

Re: [R] How to make flowchart in R?

2010-08-16 Thread Robert Baer
The R Journal article by Paul Murrell may be of interest: http://journal.r-project.org/2009-1/RJournal_2009-1_Murrell.pdf The documentation for the Diagram Package may also be of interest: http://cran.r-project.org/web/packages/diagram/diagram.pdf

Re: [R] Missing values

2010-08-17 Thread Robert Baer
- Original Message - From: Stephen Liu sati...@yahoo.com To: r-help@r-project.org Sent: Tuesday, August 17, 2010 8:25 AM Subject: Re: [R] Missing values - Original Message From: Michael Bedward michael.bedw...@gmail.com To: Stephen Liu sati...@yahoo.com Cc:

Re: [R] how can I plot bar plots with all the bars (negative and positive) in the same direction????

2010-09-03 Thread Robert Baer
As is often the case in R, just because you shouldn't do something, doesn't mean you can't do it. Still, I'd urge you to consider the visual honesty of what you propose. If you're still insistent: dat - (-3:4) dat1=dat-min(dat) barplot(dat1,axes=FALSE) axis(2,dat1,labels=dat) Typically, zero

Re: [R] how can I plot bar plots with all the bars (negative andpositive) in the same direction????

2010-09-03 Thread Robert Baer
It occurred to me after my initial post that you will need a little more fixing of your axis labels if you have data that is not trivial integers as in your example. Consider the following solution for some irrational random numbers: dat=rnorm(20) dat1=dat-min(dat)

Re: [R] Using McNemar's test to detect shifts in response from pre- topost-treatment

2010-09-14 Thread Robert Baer
McNemar is good for paired data. See ?mcnemar.test You need to get your data into the form of a matrix (e.g., help example), and you will need data organized by concordant and discordant pairing. This means that you will need to organize your data differently than you show us: # pairs YES

Re: [R] loop in R

2010-10-06 Thread Robert Baer
set.seed(180185) #loop: create 10 times the variables (u1,u2,u3,u4,u5) for (i in 1:10){ u1 - c(runif(200,0,1)) u2 - c(runif(200,0,1)) u3 - c(runif(200,0,1)) u4 - c(runif(200,0,1)) u5 - c(runif(200,0,1)) u - c(u1,u2,u3,u4,u5) mu - matrix(u, nrow=1000, ncol=1) } As you can see, when I

Re: [R] user input in R

2009-10-10 Thread Robert Baer
I'm just learning R (I don't know any other programming languages), and I have a question. I am trying to figure out how to ask for user input (say, a set of 3 numbers) then put those numbers into an array. I've looked around, but I haven't been able to find any answers that I understand.

Re: [R] SPSS long variable names

2009-10-13 Thread Robert Baer
- Original Message - From: Robert Baer rb...@atsu.edu To: Orvalho Augusto orvaq...@gmail.com Sent: Tuesday, October 13, 2009 9:52 AM Subject: Re: [R] SPSS long variable names I am wondering if there is a patch for the SPSS reading code on the foreign package, in order to be able

Re: [R] SPSS long variable names

2009-10-14 Thread Robert Baer
To: Robert Baer rb...@atsu.edu Cc: r-help@r-project.org Sent: Tuesday, October 13, 2009 10:39 AM Subject: Re: [R] SPSS long variable names No! That is variable labels. Caveman On Tue, Oct 13, 2009 at 4:52 PM, Robert Baer rb...@atsu.edu wrote: I am wondering if there is a patch for the SPSS

Re: [R] SPSS long variable names

2009-10-15 Thread Robert Baer
The problem is the limit of 8 characters long on variable names.[in read.spss??] I just created a file in SPSS 17 and saved it in standard format (.sav) The file had four unique variable NAMES which were much longer than 8 characters (although not 64). I entered two rows of data and saved

[R] Rjava, RImageJ, and/or S4 question.

2009-10-27 Thread Robert Baer
I am out of my league with this question. The following code starts the java imaging program ImageJ from within R, and displays an image (assuming ImageJ is installed on your computer). library(RImageJ) img - IJ$openImage( file.choose() ) #pick an available .tif file img$show()# make

[R] strsplit() and Windows file paths

2009-10-29 Thread Robert Baer
There are two ways to express file paths with the Windows environment: a=file.choose() a [1] C:\\Documents and Settings\\rbaer\\Desktop\\_VNT_Test\\coordFocused 20k F5 0ng Ki8751 t20.txt and b= paste(getwd(),/,dir()[1],sep=) b [1] C:/Documents and Settings/rbaer/Desktop/_VNT_Test/coordFocused

Re: [R] new help pages in R 2.10.0

2009-11-05 Thread Robert Baer
I think the html based one works fine (even if I definitely preferred the old chtml). Only another comment: I chose html in my installation with Windows Vista, but still the help pages are opened in text files in the newly installed version. The nice thing about the chtml that is missing from

Re: [R] Whiskers on the default boxplot {graphics}

2010-05-12 Thread Robert Baer
- Original Message - Fantastic! It would be great if the description could be modified to include the mysterious bit about the upper and lower bound whisker positions: upper whisker = min(max(x), Q_3 + 1.5 * IQR) lower whisker = max(min(x), Q_1 - 1.5 * IQR) -- snip --

Re: [R] Whiskers on the default boxplot {graphics}

2010-05-13 Thread Robert Baer
Hi Peter, You're absolutely correct! The description for 'range' in 'boxplot' help file is a little bit confusing by using the words interquartile range. I think it should be changed to the length of the box to be exact and consistent with those in the help file for boxplot.stats. The

Re: [R] Whiskers on the default boxplot {graphics}

2010-05-13 Thread Robert Baer
wrote: On May 13, 2010, at 10:25 AM, Robert Baer wrote: Hi Peter, You're absolutely correct! The description for 'range' in 'boxplot' help file is a little bit confusing by using the words interquartile range. I think it should be changed to the length of the box to be exact and consistent

Re: [R] Axes intercept

2010-06-28 Thread Robert Baer
- Original Message - From: KENT V.T. v.t.k...@durham.ac.uk To: r-help@r-project.org Sent: Monday, June 28, 2010 8:28 AM Subject: [R] Axes intercept I have a plot where the values of the y axis go from a positive number to a negative number and I want the x axis to intercept at zero

Re: [R] Display of results

2010-01-25 Thread Robert Baer
I don't know this test, but as written, LM - Lo.Mac(y,kvec), will just make an assignment. To display the result you could write: (LM - Lo.Mac(y,kvec)) - Original Message - From: Hichem Ben Khedhiri bk.hic...@googlemail.com To: r-help@r-project.org Sent: Sunday, January 24, 2010

Re: [R] How to learn R language?

2008-08-27 Thread Robert Baer
On 27-Aug-08 14:24:31, Charles Annis, P.E. wrote: There are many books, but if I had to choose one that teaches R and teaches statistics at the same time (Yes, you already know stats, so it will be that much easier) I'd choose Peter Dalgaard's book, Introductory Statistics with R. It's

Re: [R] A question about using function plot

2008-07-14 Thread Robert Baer
How about: plot.default(factor(x), y, xaxt='n', xlab='x') axis(1, at=factor(x), labels=x) - Original Message - From: Jason Liao [EMAIL PROTECTED] To: r-help@r-project.org Sent: Monday, July 14, 2008 10:00 AM Subject: [R] A question about using function plot Hello, everyone! I

Re: [R] Dividing by 0

2008-07-24 Thread Robert Baer
I'm trying to calculate the percent change for a time-series variable. Basically the first several observations often look like this, x - c(100, 0, 0, 150, 130, 0, 0, 200, 0) and then later in the life of the variable they're are generally no more 0's. So when I try to calculate the percent

[R] SPSS to R Data file conversion

2008-04-12 Thread Robert Baer
I wrote a little routine to convert multiple spss data files (as data frames) to R data files. The code is as follows: # list=dir(pattern=.sav) library(foreign) for (i in 1:length(list)){ # The saved data frame will be dat dat=read.spss(list[i],to.data.frame=TRUE)

Re: [R] Histogram Label Font Size

2008-04-14 Thread Robert Baer
Here's an example to make the text 18 point. You can set ps=6 to make small text. x=rnorm(100) opar=par(ps=18) # Make text 18 point hist(x) opar Rob Baer - Original Message - From: Sue Lee [EMAIL PROTECTED] To: r-help@r-project.org Sent: Monday, April 14, 2008 2:35 PM Subject: [R]

Re: [R] Mantel-Haenszel for 2x2

2008-04-19 Thread Robert Baer
See if this is what you are looking for: library(epicalc) data(Oswego) use(Oswego) cc(ill, chocolate) mhor(ill, chocolate, sex) You might also wish to read the details of: mantelhaen.test HTH], Rob Baer - Original Message - From: Dr. Jeff Miller [EMAIL PROTECTED] To: [EMAIL PROTECTED]

Re: [R] Documentation General Comments

2008-04-21 Thread Robert Baer
I realize the R developers are probably overwhelmed and have little time for this, but the documentation really needs some serious reorganizaton. A good through description of basic variable types would help a lot, e.g. the difference between lists, arrays, matrices and frames. Agreed, esp

Re: [R] ggplot: colours to geom_segments

2009-08-22 Thread Robert Baer
Do you want something like: p - ggplot(xx, aes(x = 0, xend = Expense, y = Food, yend = Food)) pa - p + geom_point(aes(Expense, Food)) + geom_segment(colour=c(red,blue,purple,darkgreen)) + xlab(Food) + geom_vline(xintercept=40, colour='red') pa - Original Message - From:

Re: [R] upgraging R from 2.9.0 to 2.9.1

2009-08-23 Thread Robert Baer
FWIW, the R homepage says, R 2.9.2 Release Candidates will appear August 17-24. Final release is scheduled for 2009-08-24. (Tomorrow?) Maybe a short upgrade delay is in order? ;-) On Windows, I routinely install new versions without removing the previous version. They always seem to coexist

Re: [R] SRS Required sample size for simulated data

2009-09-09 Thread Robert Baer
My guess is that by table you are really looking for a 'dataframe'. Further, I guess that an indexed format is most useful to you. Try: # Create the requested dataframe dat=data.frame(type=c(rep(Hypermarket,10),rep(Supermarket,15),rep(Minimarket,20),rep(Cornershop,20),rep(Spazashop,35)),

Re: [R] Help on percentage of random numbers for different classes

2009-09-10 Thread Robert Baer
I am sorry for asking this stupid question, but i have been running in circles. I want to randomly generate a scaling point of between 1 and 10, for say hundred entries, where the first 10% percent is has rates between 2 and 7, the next 15% 3 and 7, 20% between 3 and 9, 20% between 3 and 10,

[R] R Journal

2009-09-15 Thread Robert Baer
Is there a good search tool, either in R itself or in a public place, to search for R Journal (R News) articles on a given package or topic? I found the Table of Contents and the BibTex file, but I was looking for something a little more flexible. These lists are getting impressively long

Re: [R] two questions for R beginners

2010-02-26 Thread Robert Baer
Honestly what I remember as the most difficult thing when I 'first' started using R was figuring out how to read in my own datasets. I eventually discovered the R import/export manual, but somehow this alluded me initially. All the R tutorials I was working from simply generated data or used the

[R] Rcmdr plugins produce error

2010-03-18 Thread Robert Baer
I recently updated all my plugins, and for the fun of it, I added ALL the Rcmdr Plugins to my collection to see what functionality might exist. I started Rcmdr and loaded ALL the available Plugins from to the Rcdmdr Tools menu. To my suprise Rcmdr produced a number of warnings and finally an

Re: [R] use read.table for a partial reading

2010-04-10 Thread Robert Baer
- Original Message - From: Covelli Paolo pcove...@tele2.it To: r-help@r-project.org Sent: Thursday, April 08, 2010 9:05 AM Subject: [R] use read.table for a partial reading Hi everyone, I've got a matrix data with 20 variables (V1, V2, V3, ...) and 215 rows (observations). I'm

Re: [R] Bar chart in ascending order for each level of X

2011-07-16 Thread Robert Baer
The question is how to plot a bar chart in which bars are sorted in ascending order for each level of X. I would appreciate receiving your advice and help. This toy example might help: barlab - sample(LETTERS,10) values - sample(1:100, 10) op - par(mfrow= c(4,1)) barplot(values, names.arg =

Re: [R] barplot question

2011-07-19 Thread Robert Baer
As Sarah requested, could you at least read the posting guide and provide us with some sample data? -- Robert W. Baer, Ph.D. Professor of Physiology Kirksville College of Osteopathic Medicine A. T. Still University of Health Sciences 800 W. Jefferson

[R] library(foreign) read.spss warning

2011-03-25 Thread Robert Baer
I got the following: library(foreign) swal = read.spss(swallowing.sav, to.data.frame =TRUE) Warning message: In read.spss(swallowing.sav, to.data.frame = TRUE) : swallowing.sav: Unrecognized record type 7, subtype 21 encountered in system file The bulk of the data seems to read in a

Re: [R] save ordinary numerical calculations as pdf

2011-03-30 Thread Robert Baer
I'm not sure about .pdf, but look at ?sink to create a text file. You could then print this as a .pdf if you so desired. Rob -- Robert W. Baer, Ph.D. Professor of Physiology Kirksville College of Osteopathic Medicine A. T. Still University of Health

[R] corSpatial and nlme

2011-04-06 Thread Robert Baer
I noticed that ?corClasses in package nlme does not list corSpatial among the standard classes. This might either be intentional because corSpatial is not standard , or it might be simply an oversight that needs correcting. -- Robert W. Baer, Ph.D.

Re: [R] Using Java methods in R

2011-04-29 Thread Robert Baer
-- snip -- It clogs up my email, takes a long time to delete, and is hard to be selective enough to not delete some of my other important email. -- snip -- If you don't care about contributing to the R listserve community, it's hard to imagine why that community should care about you. Some

Re: [R] Watts Strogatz game

2011-05-04 Thread Robert Baer
I have a erdos-renyi game with 6000 nodes and probability 0.003. g1 = erdos.renyi.game(6000, 0.003) How to create a Watts Strogatz game with the same probability. g1 = watts.strogatz.game(1, 6000, ?, ?) What should be the third and fourth parameter to this argument. According to

Re: [R] [R ]Plot a colored rectangle under the points of a graph

2011-10-08 Thread Robert Baer
This is what you tried? What doesn't work? x=1:10 y=1:10 plot(x,y, type='n') rect(4,0,6,11, col=5) points(x,y) ___ Hello everyone, if I have: x=1:10 y=1:10 plot(x,y) and I plot a rectangle rect(4,0,6,11, col=5) it covers the points of the graph. Is there a way to draw the

Re: [R] How to run Rcmdr with OS 10.4?

2011-10-12 Thread Robert Baer
As the posting guide tells you, the first thing to try is to install the latest version of R. -Original Message- From: roche...@free.fr Sent: Tuesday, October 11, 2011 5:44 AM To: r-help@r-project.org Subject: [R] How to run Rcmdr with OS 10.4? I've installed Rcmdr package and it

Re: [R] Cleaning up messy Excel data

2012-02-28 Thread Robert Baer
-Original Message- From: Noah Silverman Sent: Tuesday, February 28, 2012 3:27 PM To: r-help Subject: [R] Cleaning up messy Excel data Unfortunately, some data I need to work with was delivered in a rather messy Excel file. I want to import into R and clean up some things so that I can

Re: [R] ggplot2 dot chart-start at zero

2012-02-29 Thread Robert Baer
Not sure I'm interpreting what you want to do correctly, but how about: p2 - pa + scale_x_continuous(limits=c(0,80)) p2 -Original Message- From: John Kane Sent: Wednesday, February 29, 2012 2:32 PM To: r-help@r-project.org Subject: [R] ggplot2 dot chart-start at zero I am trying to

[R] Windows - **Paste commands only** issue

2012-03-06 Thread Robert Baer
In following a thread on this mailing list, I encounter an apparent issue with **Edit | Paste commands only** in the Windows R-GUI. Reproducible steps: Go to (using IE?): http://www.r-bloggers.com/select-operations-on-r-data-frames/ Find the section entitled 'Duplicate row names' Copy the

Re: [R] sort dates

2012-03-09 Thread Robert Baer
How is it possible to sort dates in R? Try this: a = sample(as.Date(1:100, origin = '2012-01-01'),15) a [1] 2012-01-31 2012-01-22 2012-03-18 2012-03-05 2012-03-17 [6] 2012-03-08 2012-01-08 2012-01-20 2012-03-01 2012-03-21 [11] 2012-02-17 2012-01-17 2012-02-12 2012-02-28 2012-04-01 sort(a)

Re: [R] Identifying a change in events between bins

2012-03-19 Thread Robert Baer
Your description was too general for me to know exactly what you want but perhaps this will help you solve your own problem set.seed(123) evtlist = sample(c('fwd','rev'),100,replace=TRUE) evtlist [1] fwd rev fwd rev rev fwd rev rev rev fwd rev fwd rev [14] rev fwd rev fwd fwd fwd rev rev rev

Re: [R] How to get the input of a function right?

2012-03-25 Thread Robert Baer
Your definition of x1, x2 and x3 requires the c() function. Try leaving some space around operators for increased readability. Remember that the that the result of operations on sets containing NA is often set to NA as well. My guess is that if you print out x it contains NA and that your

Re: [R] Problem loading package 'JGR' using R-2.15.0 (Win32).

2012-04-02 Thread Robert Baer
On a 32-bit windows 7 machine running R 2.15.0 standard binary install, JGR loaded fine as seen below: library(JGR) Loading required package: rJava Loading required package: JavaGD Loading required package: iplots Please type JGR() to download/launch console. Launchers can also be obtained at

Re: [R] Histogram from a table in R

2012-04-03 Thread Robert Baer
intv = c('0-19','10-19','20-29','30-39') cnts = c(0, 3117, 4500, 2330) barplot(cnts, space=0, names = intv, xlab='Age Range', ylab = 'Counts', main='My Histogram') -Original Message- From: gina_alessa Sent: Tuesday, April 03, 2012 4:08 PM To: r-help@r-project.org Subject: [R]

Re: [R] legend not appearing in Word document

2010-12-16 Thread Robert Baer
On 12/12/2010 8:59 AM, Tim Clark wrote: I need help with using graphics in Word 2007 that will later be converted into a pdf document. I have tried several formats and found that I get the best quality of graphics using .wmf, .eps format, but when I convert it to .pdf I get a bunch of lines

Re: [R] Assumptions for ANOVA: the right way to check the normality

2011-01-05 Thread Robert Baer
Someone suggested me that I don´t have to check the normality of the data, but the normality of the residuals I get after the fitting of the linear model. I really ask you to help me to understand this point as I don´t find enough material online where to solve it. Try the following: #

Re: [R] how to get old packages to work on R 2.12.1

2011-01-19 Thread Robert Baer
Did you do: update.packages(checkBuilt=TRUE) Rob -- From: Joshua Wiley jwiley.ps...@gmail.com Sent: Wednesday, January 19, 2011 11:43 AM To: Joseph Boyer joseph.g.bo...@gsk.com Cc: r-help@r-project.org Subject: Re: [R] how to get old packages to

Re: [R] ANOVA plotting

2011-01-20 Thread Robert Baer
I recently started using R and I have a simple question. I am running R (v. 2.12.1) and Rcmdr (v.1-6.3) on Mac (Snow Leopard). I am using a data set I used before for practicing ANOVA with R, so I know what the results should look like. I can get ANOVA table using both Rcmdr and GUI. However,

Re: [R] new to R coding.

2011-10-25 Thread Robert Baer
Not quite sure what you were trying to do, but try the slight modification below and see if it works for you: x = c (2434.1, 2463.7, 2451.6, 2444.5, 2431.3, 2436.3, 2412.6, 2417.9, 2380.1, 2366.2, 2349.1, 2373.9, 2336.9, 2335.0, 2297.1, 2291.1, 2278.6, 2289.9, 2314.5, 2328.8, 2322.0,

Re: [R] read.table question

2011-12-08 Thread Robert Baer
Hello All, This works, results - read.table(plink.txt,T) while this doesn't. results - read.table(plink.txt) The T is the value for the second parameter which you show from the help file printout you looked at to be header. Thus, you are writing in short cut: results -

Re: [R] a quick question about rbinom

2012-01-04 Thread Robert Baer
-Original Message- From: lynn.tsai Sent: Wednesday, January 04, 2012 3:38 PM To: r-help@r-project.org Subject: [R] a quick question about rbinom Hello, I have the following code using rbinom, but I don't understand what *+1* means in the code. Could someone help? Thanks so much,

Re: [R] help to slip a file name using strsplit function

2012-01-25 Thread Robert Baer
How about: unlist(strsplit(as.character(Myfile_MyArea1_sample1.txt), split = '[[:punct:]]')) -Original Message- From: gianni lavaredo Sent: Wednesday, January 25, 2012 9:26 AM To: r-help@r-project.org Subject: [R] help to slip a file name using strsplit function Dear Researchers, I

Re: [R] x-axis label for boxplot

2012-02-08 Thread Robert Baer
xlab = needs to be names = see ?boxplot boxplot(jdata[,x],jdata[,y],jdata[,z],names=c(x values,yvalues,zvalues)) -Original Message- From: John Sorkin Sent: Wednesday, February 08, 2012 6:48 PM Cc: r-help Subject: [R] x-axis label for boxplot I am trying to produce three

[R] Creating Observation ID

2011-05-09 Thread Robert Baer
If I have a data frame something like: Value=rnorm(30) Group = sample(c('A','B','C'), 30, replace=TRUE) df = data.frame(Value, Group) It seems like it should be simple to create an 'ObsID' column which indicates the observation order of each Value within each of the 3 groups. Somehow, I can't

Re: [R] Creating Observation ID

2011-05-10 Thread Robert Baer
Kirksville College of Osteopathic Medicine A. T. Still University of Health Sciences 800 W. Jefferson St. Kirksville, MO 63501 660-626-2322 FAX 660-626-2965 -- From: William Dunlap wdun...@tibco.com Sent: Monday, May 09, 2011 5:17 PM To: Robert Baer rb

Re: [R] Excel to R

2011-05-13 Thread Robert Baer
-- From: Asan Ramzan asanram...@yahoo.com Sent: Friday, May 13, 2011 6:02 AM To: r-help@r-project.org Subject: [R] Excel to R Hello R-help I am trying to copy and paste large column of data from windows Excel into R base speed sheet using:

Re: [R] recursive function

2011-05-19 Thread Robert Baer
Perhaps this is useful: x=c(-2,0,2) sign(x)*abs(x) [1] -2 0 2 -- Robert W. Baer, Ph.D. Professor of Physiology Kirksville College of Osteopathic Medicine A. T. Still University of Health Sciences 800 W. Jefferson St. Kirksville, MO 63501 660-626-2322

Re: [R] [r] regression coefficient for different factors

2011-05-20 Thread Robert Baer
?summary produces r^2 in 2nd to last line, as in, set.seed(12); a=rnorm(100); b = runif(100); c = factor(rep(c('No', 'Yes'),50)); df = data.frame(a,b,c) head(df) a b c 1 -1.4805676 0.9729927 No 2 1.5771695 0.2172974 Yes 3 -0.9567445 0.5205087 No 4 -0.9200052 0.8279428

Re: [R] Group close numbers in a vector

2011-05-21 Thread Robert Baer
Hi everyone, i am trying to group close numbers in a vector. For example i have a vector x = [1 2 4 7 9 10 15]. I want the code to pick 1 2 4 (max difference between successive numbers is 2) and assign them to variable a, then pick 7 9 10 and assign them to b and 15 to c. But since i do not

Re: [R] Data values on graphics

2011-05-24 Thread Robert Baer
zcatav wrote: How can i add data values or proportion values on bar or pie charts? You can add any text or numbers to any place of a graph with this command: text(45,20,some text) where: 45 is position on the X-axis and 20 - on the Y. Note that numbers denote the MIDDLE of the text,not its

Re: [R] Normality test

2011-05-27 Thread Robert Baer
I am writing to inquire about normality test given in nortest package. I have a random data set consisting of 300 samples. I am curious about which normality test in R would give me precise measurement, whether data sample is following normal distribution. As p value in each test is different in

Re: [R] Line Graphs

2011-06-07 Thread Robert Baer
I want to plot 6 line graphs. I have 10 points 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0. At each point say 0.1, I have 6 variables A, B, C, D, E and F. The variables all have values between 0 and 1 (and including 0 and 1). I also want to label the x axis from 0.1 to 1.0 and the y

Re: [R] Line Graphs

2011-06-07 Thread Robert Baer
Thanks Rob, but the legend is not appearing in the plot. I think the best place for it is on the top left. Is there anyway I can also get it broken down in tenths instead of fifths? The legend is easy; just specify where you want it. The first 2 parameters specify the x, y of the top left

Re: [R] Running R from windows command prompt

2011-06-28 Thread Robert Baer
Subject: Re: [R] Running R from windows command prompt Hi Siddharth, many experts already answered your query, however I would like to share how I run R in command prompt: 1. open command prompt 2. change working directory: cd C:\\R-2.13.0\bin\i386 (put the entire path here,

Re: [R] rJava works with 32-bit but not 64

2013-02-04 Thread Robert Baer
Is it feasible to have both installed in a way that allows the each version of R to select its own version of Java? A comment on stackoverflow suggests that may not be easy (http://stackoverflow.com/questions/5272216/is-it-possible-to-install-both-32bit-and-64bit-java-on-windows-7).

Re: [R] rJava works with 32-bit but not 64

2013-02-06 Thread Robert Baer
For what it is worth Spencer, I can start rJava in both 32-bit R and 64-bit R for Windows 7. [And could even before Simon fixed the error message). And yes, I have both 32-bit Java and 64-bit Java 1.7.0_13 installed. They should be separate entries under your control panel. Rob --

Re: [R] Generating 1-bit and 8-bit BMP files using R

2013-03-05 Thread Robert Baer
On 3/5/2013 12:51 AM, Ingo Reinhold wrote: Hi, I'm trying to use the data which I generate within R to make images in .bmp format to be lateron printed by a printer. My first thought was the RImageJ package, but this seems to be discontinued. What I am currently doing is generating a matrix

Re: [R] how to check given number seq. is time series or not?

2012-05-29 Thread Robert Baer
Perhaps ?str is the command you seek? -Original Message- From: sagarnikam123 Sent: Saturday, May 26, 2012 10:24 PM To: r-help@r-project.org Subject: Re: [R] how to check given number seq. is time series or not? Yes,sir index of above/below numbers is time (both dataset are same)

[R] ggplot2, grid graphics, x11(), windows(), and device fonts

2012-06-05 Thread Robert Baer
I am trying to use ggplot() to produce a graph and am getting warnings that I don't understand. This is happening from within RStudio, but also happens if I start Windows R GUI. Can anyone help me understand the warning and how to make sure the right fonts are designated? The relevant code

Re: [R] R2wd error in wdGet

2012-06-10 Thread Robert Baer
wdGet() Error in if (wdapp[[Documents]][[Count]] == 0) wdapp[[Documents]]$Add() : argument is of length zero Not sure what your error means without more context, you should see: Loading required package: rcom Loading required package: rscproxy Do you have RDCOMClient installed and working

Re: [R] How to plot linear, cubic and quadratic fitting curve in a figure?

2012-06-13 Thread Robert Baer
dput(test) structure(list(sp = c(4L, 5L, 9L, 12L, 14L), env = c(12L, 18L, 20L, 17L, 15L)), .Names = c(sp, env), class = data.frame, row.names = c(NA, -5L)) plot(test$sp~test$env, main = S vs. temp, xlim=c(0,20), ylim=c(0,14), ylab=S,xlab=env) linear-lm(test$sp~test$env)

Re: [R] binary tree

2012-07-02 Thread Robert Baer
On 6/27/2012 3:20 AM, Peppino wrote: Hi I am new with R I Have to build a binary tree with R. I'm very confused was wondering if anyone had any R sample code they would share. Any bady can help me? You might want to look at the R Task view for phylogenetics:

[R] RMySQL install on windows

2012-10-09 Thread Robert Baer
I have been trying to install RMySQL on Windows 7 following the procedure at: http://biostat.mc.vanderbilt.edu/wiki/Main/RMySQL I think I have properly installed RTools and created a proper Renviron.site file saying: MYSQL_HOME=C:/Program Files/MySQL/MySQL Server 5.5 When I try to install

Re: [R] RMySQL install on windows

2012-10-12 Thread Robert Baer
On 10/10/2012 12:58 PM, Greg Snow wrote: I finally was able to compile/load it under windows 7. I had similar problems to what you show below. I set the MYSQL_HOME environmental variable through windows (start button control panel System and Security system Advanced System Settings

[R] Java, rJava, and Windows x64

2012-10-29 Thread Robert Baer
When running [1] R version 2.15.1 (2012-06-22) x86_64-pc-mingw32, rJava fails. I have installed both the 32-bit and 64-bit versions of Java 7 update 9. library(rJava) Error : .onLoad failed in loadNamespace() for 'rJava', details: call: stop(No CurrentVersion entry in ', key, '! Try

[R] rgl package and animation

2012-11-02 Thread Robert Baer
I am trying to figure out how to use rgl package for animation. It appears that this is done using the play3d() function. Below I have some sample code that plots a 3D path and puts a sphere at the point farthest from the origin (which in this case also appears to be at the end of the path).

Re: [R] rgl package and animation

2012-11-03 Thread Robert Baer
On 11/3/2012 6:47 AM, Duncan Murdoch wrote: On 12-11-02 7:47 PM, Robert Baer wrote: I am trying to figure out how to use rgl package for animation. It appears that this is done using the play3d() function. Below I have some sample code that plots a 3D path and puts a sphere at the point

Re: [R] Saving R Graph to a file

2012-11-04 Thread Robert Baer
Some hints: For pdf(), height and width are in inches, not pixels. dev.off() is necessary after drawing the image for pdf(). The name for the file argument (file=c:/figure.xxx) is file not filename hist(CO2[,5]) is more interesting And yes, ?pdf ?postscript ?ping On 11/3/2012 11:16 PM,

Re: [R] Saving R Graph to a file

2012-11-04 Thread Robert Baer
On 11/4/2012 4:32 AM, Robert Baer wrote: Some hints: For pdf(), height and width are in inches, not pixels. dev.off() is necessary after drawing the image for pdf(). The name for the file argument (file=c:/figure.xxx) is file not filename hist(CO2[,5]) is more interesting And yes, ?pdf

Re: [R] rgl package and animation

2012-11-05 Thread Robert Baer
-- snip -- On 11/4/2012 7:45 AM, Duncan Murdoch wrote: First, draw the new sphere at the first point and save the object id: sphereid - sphere3d(dat[1,c(X, Y, Z)], col=red, radius=1) # Also save the spinner that you like: spin - spin3d( ) #maybe with different parms # Now, the animation

Re: [R] exporting 3D dynamic graph

2012-11-05 Thread Robert Baer
On 11/5/2012 8:23 AM, Christophe Genolini wrote: Hi the list, Using misc3d, we can export 3d dynamic graph in pdf format. Is it also possible to export these graph into a format that we can publish on the web? Christophe You don't provide enough information to know exactly what your needs

Re: [R] Rank Amateur Question

2013-05-29 Thread Robert Baer
You probably want from windows GUI: source(test.R) To read help to learn about this command type: ?source at the command prompt. (Similar pattern get help on other commands too - its an important R skill/habit). If your are going to do a lot of script development, I would highly recommend

Re: [R] how to install R 3.0.1

2013-05-31 Thread Robert Baer
On 5/31/2013 11:59 AM, Ranjan Maitra wrote: Hi John, I suspect you may be missing a c()? On Fri, 31 May 2013 06:05:45 -0800 John Kane jrkrid...@inbox.com wrote: paks - install.packages( Hmisc, plyr) paks - install.packages( c(Hmisc, plyr)) install.packages(paks) You can use the

Re: [R] Workspace Question.

2012-07-23 Thread Robert Baer
The objects from an R session are saved in the RWorking directory of the session. The answer to your question will depend on whether you started the different versions of R using shortcuts located in different folders or the same folder. The objects should not be automatically deleted so I

Re: [R] R2wd package wdGet() error

2012-07-23 Thread Robert Baer
On 7/23/2012 4:25 PM, Jean V Adams wrote: I am having trouble using the R2wd package. The last time I used it successfully, I was running an earlier version of R and an earlier version of Word with an earlier Windows OS. I'm not sure which if any of these changes might be contributing to the

Re: [R] Turning off continuation prompt?

2012-07-30 Thread Robert Baer
On Mon, 30 Jul 2012 09:58:09 +0100 (BST) (Ted Harding) ted.hard...@wlandres.net wrote: Greetings All. My apologies for a question whose answer is probably readily available somewhere (for some interpetation of somewhere) ... Say I have just typed (from a sheet of paper) several lines into the R

Re: [R] ANOVA repeated measures and post-hoc

2012-08-16 Thread Robert Baer
Check out: http://rtutorialseries.blogspot.com/2011/02/r-tutorial-series-two-way-repeated.html On 8/15/2012 11:32 AM, Diego Bucci wrote: Hi, I performed an ANOVA repeated measures but I still can't find any good news regarding the possibility to perform multiple comparisons. Can anyone help

Re: [R] return first index for each unique value in a vector

2012-08-28 Thread Robert Baer
On 8/28/2012 5:52 PM, Bert Gunter wrote: Sheesh! I would have thought that someone would have noticed that on the ?unique Help page there is a link to ?duplicated, which gives a _logical_ vector of the duplicates. From this, everything else can be quickly derived -- and packaged in a simple

  1   2   >