### Re: [R] boxplot issue

A) you hijacked another thread. On Oct 15, 2010, at 9:50 AM, Jonas Josefsson wrote: Hi! I am trying to produce a graph which shows overlap in latitude for a number of species. I have a dataframe which looks as follows species1,species2,species3,species4. minlat

### Re: [R] p-values from coxph?

On Oct 15, 2010, at 9:21 AM, Öhagen Patrik wrote: Dear List, I each iteration of a simulation study, I would like to save the p- value generated by coxph. I fail to see how to adress the p-value. Do I have to calculate it myself from the Wald Test statistic? No. And the most important

### Re: [R] interaction contrasts

Should you need to do it again, you may want to look at the relevel function. I suppose that would meet the definition of some versions of on the fly but once I have a model, rerunning with a different factor leveling is generally pretty painless. -- David. On Oct 15, 2010, at 9:09 AM,

### [R] [R-pkgs] solaR: version 0.20

From: Oscar Perpiñan Lamigueiro oscar.perpi...@upm.es Date: October 8, 2010 4:25:34 AM EDT To: r-packa...@r-project.org Subject: [R-pkgs] solaR: version 0.20 Hello, I'd like to announce the availability of the version 0.20 of the solaR package. It provides a set of calculation methods of solar

### Re: [R] delete data row

IRD; There is a danger in applying logical tests of equality to floating point numbers. It may be safer to use all.equal or zapsmall in the construction of your tests. all.equal( (2^(0.5))^2 , 2) [1] TRUE (2^(0.5))^2 == 2 [1] FALSE -- David. On Oct 16, 2010, at 1:30 PM, Joshua Wiley

### Re: [R] reduce the size of points in plot???

On Oct 16, 2010, at 12:12 PM, Hongwei Dong wrote: Hi, R users, Can anyone tell me how I can change the size of points in my plot? For example: x - c(1,3,6,9,12) y - c(1.5,2,7,8,15) plot(x,y,pch=20) How do I reduce the size of those points? plot(x,y,pch=20, cex=0.2) plot(x,y,pch=20,

### Re: [R] delete data row

at 12:12 PM, David Winsemius dwinsem...@comcast.net wrote: IRD; There is a danger in applying logical tests of equality to floating point numbers. It may be safer to use all.equal or zapsmall in the construction of your tests. all.equal( (2^(0.5))^2 , 2) [1] TRUE (2^(0.5))^2 == 2 [1] FALSE

### Re: [R] delete data row

On Oct 17, 2010, at 3:56 PM, William Dunlap wrote: I had been thinking of: x - c(1, (2^(0.5))^2 , 3, 5, (2^(0.5))^2 , 3, 1) y - 2 x[-which(zapsmall(x-y) == 0)] [1] 1 3 5 3 1 Using which() to convert logicals into integer subscripts is almost always unnecessary and often wrong. At one

### Re: [R] read.spss warning message (Unrecognized record type 7, subtype 18 encountered in system file)

PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

### Re: [R] How to read only ten rows from a SAS dataset (read.ssd)?

/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

### Re: [R] readLines: how to make a data.frame?

-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

### Re: [R] Error: object 'short' not found

and the the ^^^ will not be aligned with the word short. You should fix the settings of you mail client in both these areas. Viki [[alternative HTML version deleted]] David Winsemius, MD West Hartford, CT

### Re: [R] Chron object in time series plot

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing

### Re: [R] Chron object in time series plot

On Oct 19, 2010, at 12:19 PM, Manta wrote: David Winsemius wrote: You seen to be under the mistaken impression that the internal representation of DateTime classes of 08:00 would be 8. Since the internal representation of time is in seconds, the even number hours would be at integer

### Re: [R] How to read only ten rows from a SAS dataset (read.ssd)?

internal variable named something like _N_ that is a line number. Perhaps ... , if=_N_ = 10, OR: -- using SAS to output a smaller file. OR: --- adding a SAS line to the output returned within the read.ssd function. -- David. -J 2010/10/19 David Winsemius dwinsem...@comcast.net

### Re: [R] Chron object in time series plot

mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r

### Re: [R] Problems with a specific calculate.

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list

### Re: [R] scatter.smooth() fitted by loess

-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

### Re: [R] Gini Coefficient

coefficient and add them together? Once again i am sorry if this is completely the wrong place to ask such a question. Peter David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

### Re: [R] How to select not continous rows?

On Oct 20, 2010, at 9:25 AM, skan wrote: Hello How can I select several not continuous rows ? If I wanted to select rows 1 to 7 I'll write mydata[,1:7] But what if I need to select rows 1 to 5 and 10 to 15? mydata[ , c(1:5, 10:15) ] -- David Winsemius, MD West Hartford, CT

### Re: [R] number format, writing 1e-5 instead of 0.00001

= 3) 0.1 [1] 0.1 0.001 [1] 1e-07 0.01 [1] 0.01 Other options for output: ?format ?sprintf -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

### Re: [R] CI using ci.numeric

-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help

### Re: [R] (no subject)

(mustart) : Value -1717986918 out of range (0, 1) Look more carefully at your 28th case. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

### Re: [R] How to select not continous rows?

On Oct 20, 2010, at 10:01 AM, David Winsemius wrote: On Oct 20, 2010, at 9:25 AM, skan wrote: Hello How can I select several not continuous rows ? If I wanted to select rows 1 to 7 I'll write mydata[,1:7] But what if I need to select rows 1 to 5 and 10 to 15? mydata[ , c(1:5, 10:15

### Re: [R] remove black square from factor plot

is produced by each species, there is a small black square in center position of the box. Is it possible to remove this. Have tried to include bxp(par(medpch=NA)) without success. Thanks! Jonas Josefsson David Winsemius, MD West Hartford, CT __ R

### Re: [R] create a list fails

])), total.patientnums.trt1=dat2[ ,2], total.patientnums.trt23=dat2[ ,2], num.countstrt1=dat2[ ,5], num.countstrt23=dat2[ ,6] ) === Dr. Jim Maas David Winsemius, MD West Hartford, CT __ R-help

### Re: [R] Changing sign on absolute numbers 0 problems

function: nonneg - function(x) { 0 + !sign(x) } nonneg(0) [1] 1 Any help is appreciated Thank you sadz ps please email me if you need more info -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch

### Re: [R] Changing sign on absolute numbers 0 problems

On Oct 20, 2010, at 3:19 PM, David Winsemius wrote: On Oct 20, 2010, at 11:47 AM, Sadz A wrote: Hi, I am trying to do some calculations turning DMS data to decimal degrees using the formula (D+(M/60)+(S/3600)), some of the D's involve -ve numbers, the easiest way to do the calculation

### Re: [R] question about masked object

: namespace:Hmisc base::units function (x) UseMethod(units) environment: namespace:base Thank you and kind regards, gigi [[alternative HTML version deleted]] -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list

### Re: [R] loading workspace- getting annoying

also skipped loading my utility functions from .Rprofile as well as skipping the workspace .Rdata file. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

### Re: [R] rounding up (always)

, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

### Re: [R] All other variables in upper scope arg for stepAIC

21.276 3.2365 + y11 0.149400 21.126 5.0956 + y31 0.025849 21.250 5.2122 Start: AIC=1 get(paste(y, i, sep = )) ~ 1 Df Sum of SqRSSAIC none 19.027 1.0023 + y11 0.227457 18.799 2.7618 + y21 0.023117 19.004 2.9780 -- David Winsemius, MD

### Re: [R] Efficient nested loops

the max of the last slice of the variable. Thanks in advance! -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

### Re: [R] Different time between date() and Sys.date()

, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

### Re: [R] 3D-scatterplots - high quality rendering?

simply create a collection of images that I animate elsewhere? require(rgl) ?movie3d -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

### Re: [R] Efficient nested loops

/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

### Re: [R] Changing sign on absolute numbers 0 problems

will be appreciated. Thank you sadz From: David Winsemius dwinsem...@comcast.net To: David Winsemius dwinsem...@comcast.net Cc: Sadz A sadz_a1...@yahoo.co.uk; r-help@r-project.org Sent: Wed, 20 October, 2010 20:45:39 Subject: Re: [R] Changing sign on absolute numbers 0 problems On Oct 20, 2010, at 3

### Re: [R] exact pattern match in grep for column headers

-- View this message in context: http://r.789695.n4.nabble.com/grep-tp3004422p3005885.html Sent from the R help mailing list archive at Nabble.com. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch

### Re: [R] data.frame query

/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

### Re: [R] Different time between date() and Sys.date()

setting and changing it in my system does the trick: Sys.time() [1] 2010-10-21 19:48:42 CEST -- David. Do you know where can I find a list of the timezone ? Message du 21/10/10 15:07 De : David Winsemius A : omerle Copie à : r-help@r-project.org Objet : Re: [R] Different time between

### Re: [R] How to access values in s4 method

3 dim(ovrl) [1] 3 3 -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal

### Re: [R] printing a variable during a loop

(0,0.025,0.25,0.5,0.75,0.975,1)) print(i) } I'd like to have the value of i printed for each loop (step). As I could see the values of i are shown on screen only after all the work is done. Thanks in advance for any suggestion. Best regards, Antonio -- David Winsemius, MD West Hartford, CT

### Re: [R] previous business day

On Oct 21, 2010, at 9:38 PM, Li, Jing Yi wrote: How to get the previous business day in R? I saw some post about using functions in timeSeries package before but can not find it anymore. ??holiday Perhaps you should look in the tis package. Thanks! David Winsemius, MD West Hartford

### Re: [R] previous business day

}} __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT

### Re: [R] R step-by-step execution

On Oct 22, 2010, alais wrote: |Moreover can you please tell me how I can concatenate variables and strings so to print some debugging messages / The value of x is + x + and the value of y is +y. Commas, not + cat(The value of x is , x , and the value of y is , y) -- David Winsemius

### Re: [R] superscript characters in title with '+'

^2+) 'onsets'), xlab = 'sec') David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

### Re: [R] issue with Matrix package

de la classe derivedDefaultMethod any idea on how I could solve this issue ? Update, first. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

### Re: [R] lm looking for weights outside of the user-defined function

of formula. Why is it looking outside the function for the object that has just been defined inside the function? David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

### Re: [R] lm looking for weights outside of the user-defined function

the rules, I just follow them. I agree that one might guess that to be the search order, but it is not what is documented. -- David Winsemius, MD West Hartford, CT Dimitri On Fri, Oct 22, 2010 at 9:15 AM, David Winsemius dwinsem...@comcast.net wrote: On Oct 22, 2010, at 9:01 AM, Dimitri

### Re: [R] superscript characters in title with '+'

-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

### Re: [R] Confusing error statement

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org

### Re: [R] lm looking for weights outside of the user-defined function

On Oct 22, 2010, at 12:17 PM, William Dunlap wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of David Winsemius Sent: Friday, October 22, 2010 6:25 AM To: Dimitri Liakhovitski Cc: r-help Subject: Re: [R] lm looking

### Re: [R] wait for graph to finish plotting

On Oct 22, 2010, at 12:50 PM, n...@aleblanc.cotse.net wrote: Tal Galili tal.gal...@gmail.com writes: I suspect that using dev.copy2eps Is not going to help you here. Please try again using: pdf(...) # Check: ?pdf for(i in something) { plot(things) } dev.off() But give pdf() the

### Re: [R] read.table input array

On Oct 22, 2010, at 8:40 PM, Balpo wrote: Hello Jim How can I ensure this reading it from the file? The thing is, I cannot use the textConnection(), because the real file has millions of rows. The advice being offered is not that you use textConnection. That is just for illustration.

### Re: [R] contour on a simplex

If Dennis' interpretation of your request is correct then Harrell rms/ Hmisc packages provide a mechanism for restricting plotting of contours to the region within a perimeter. He has a perimeter function and the bplot function takes a perim argument. The reason I did not offer it earlier

### Re: [R] help

On Oct 23, 2010, at 4:13 AM, Liviu Andronic wrote: Hello On Sat, Oct 23, 2010 at 10:43 AM, Mastaki Kambale jkmast...@hotmail.com wrote: Dears R I am a self taught novice user of R. I begin to understand its philosophy and some basics like objects, vectors, arrays, lists etc...I still

### Re: [R] convert wind direction from degrees to basic compass directions

On Oct 23, 2010, at 1:40 PM, will phillips wrote: Hello, I have a data set that includes a predictor variable wind direction. This variable is in degrees with 0=North. I've put this predictor into a linear model and its coefficient is far from significant. I was thinking about

### Re: [R] Long model formulae

On Oct 23, 2010, at 6:56 PM, Gabor Grothendieck wrote: On Sat, Oct 23, 2010 at 9:51 PM, James Hirschorn james.hirsch...@hotmail.com wrote: What is a good way to enter a very long model formula. For example: y ~ Input.2 + Input.3 + ... + Input.1000 (assuming the corresponding dataframe has

### Re: [R] If Statement Help

On Oct 23, 2010, at 7:44 PM, Jason Kwok wrote: Thanks Jorge. It works. Is there a way to keep the actual price in the price column instead of TRUE/FALSE but filtering on when price100? Huh? When I use subset I get what you ask for: subset(x, Price 100) Price 2010-10-12 101

### Re: [R] Contour Plot on a non Rectangular Grid

On Oct 24, 2010, at 4:30 AM, Lorenzo Isella wrote: Dear All, I would like to plot a scalar (e.g. a temperature) on a non- rectangular domain (or even better: I would simply like to be able to draw a contour plot on an arbitrary 2D domain). I wonder if there is any tool to achieve that

### Re: [R] Contour Plot on a non Rectangular Grid

On Oct 24, 2010, at 6:12 AM, Lorenzo Isella wrote: As to the domain of the function, at least in case (1), that should arise from the collected data points in (x,y) if the sampling is dense enough. And that is precisely what you get from the perimeter function. The earlier Design

### Re: [R] Contour Plot on a non Rectangular Grid

On Oct 24, 2010, at 9:30 AM, Lorenzo Isella wrote: Hi, And thanks for helping. I am anyway a bit puzzled, since case (1) is not only a matter of interpolation. Probably the point I did not make clear (my fault) is that case (1) in my original email does not refer to an irregular grid

### Re: [R] Contour Plot on a non Rectangular Grid

On Oct 25, 2010, at 3:41 AM, Lorenzo Isella wrote: On 10/25/2010 01:32 AM, David Winsemius wrote: You were advised to look at rms. Why have you dismissed this suggestion? Using your data setup below and packaging into a dataframe. require(rms) ddf - datadist(xysf - as.data.frame(xys

### Re: [R] How to save R file into specific type

On Oct 25, 2010, at 1:07 AM, zhiji19 wrote: Hello everyone Can you please teach me how to save my homework as .R file? ?savehistory # at least on a Mac (If you wanted the whole console session with output, then it would be select-all, copy, paste into a text editor, save as a text

### Re: [R] Artifacts in filled.contour+pdf

On Oct 25, 2010, at 6:50 AM, Mario Valle wrote: Dear all, I'm using R 2.12.0 on Windows 7 (32bits) I created a filled contour from the attached data using the following code: load('bug.RData') pdf('bug.pdf', width=14, height=7) filled.contour(o4$x, o4$y, o4$z, color=rainbow, xlim=c(4,18),

### Re: [R] divide column in a dataframe based on a character

On Oct 25, 2010, at 8:56 PM, Daisy Englert Duursma wrote: Hello, If I have a dataframe: example(data.frame) zz- c (aa_bb ,bb_cc ,cc_dd,dd_ee,ee_ff,ff_gg,gg_hh,ii_jj,jj_kk,kk_ll) ddd - cbind(dd, group = zz) and I want to divide the column named group by the _, how would I do this?

### Re: [R] (no subject)

On Oct 25, 2010, at 11:17 PM, Arsalan Fathi wrote: hello. how are you? my name is arsalan. I'm from iran. i want to write the program that done random walk(one variable and two variable). please help me. You should learn to use a search engine that is specific to R. Here's an example:

### Re: [R] Setting constraints in the glm package

On Oct 26, 2010, at 7:50 AM, Viechtbauer Wolfgang (STAT) wrote: The constraint b1=b2 in a model such as b0 + b1 x1 + b2 x2 + b3 x3 implies that b0 + b1 (x1 + x2) + b3 x3, so just add x1 and x2 (call this x12) and fit the model b0 + b1 x12 + b3 x3 and you have imposed the constraint that

### Re: [R] Forcing results from lm into datframe

On Oct 26, 2010, at 8:08 AM, Small Sandy (NHS Greater Glasgow Clyde) wrote: Hi I need some help getting results from multiple linear models into a dataframe. Let me explain the problem. I have a dataframe with ejection fraction results measured over a number of quartiles and grouped

### Re: [R] Zoom in in a plot

On Oct 26, 2010, at 8:30 AM, Alaios wrote: in a simple plot. When i do plot is it possible to zoom in or out or this is not possible at all? Zoom? Do you mean restrict the region plotted to specific ranges? xlim and ylim arguments provide that facility. -- David. Best Regards Alex

### Re: [R] Setting constraints in the glm package

On Oct 26, 2010, at 9:27 AM, David Smith wrote: Many thanks for the help. You could express your thanks by including context the next time you present a follow-up (as requested in the Posting Guide). Only a minority of readers view this list on Nabble, so we don't see the web

### Re: [R] Forcing results from lm into datframe

thought maybe I sould see a cbind() call in there, and I suppose this section as.call(c(expression(data.frame), x ... may have that effect -- David. Why doesn't fitsdf - as.data.frame(t(fits)) work? Sandy Small From: David Winsemius [dwinsem

### Re: [R] superscript characters in title with '+'

On Oct 26, 2010, at 11:15 AM, Gavin Simpson wrote: On Fri, 2010-10-22 at 15:39 +0200, Claudia Beleites wrote: On 10/22/2010 03:15 PM, DrCJones wrote: snip / Being a chemist, it seemed natural to me to put the i after the concentration brackets into a subscript - though you didn't say

### Re: [R] Cox Proportional Models and Haplotypes

On Oct 26, 2010, at 8:57 AM, sr500 wrote: Hello, I was wondering if anyone knew of a function that fits haplotype data into a cox proportional hazard model. I have computed my Haplotype frequencies using the haplo.stats package. I have also been using the haplo.glm function but this is

### Re: [R] ggplot - unwanted sorted X values

as is typically done at time of factor creation. You probably need to work with factor levels. (All guesswork in the absence of a reproducible example.) -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch

### Re: [R] coxph linear.predictors

in this area will get expert review and correction. Your comments are appreciated. Stephen Bond -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Wednesday, October 27, 2010 1:15 PM To: Bond, Stephen Cc: r-help@r-project.org Subject: Re: [R] coxph

### Re: [R] How do I read multiple rows of different lengths?

://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

### Re: [R] (no subject)

On Oct 27, 2010, at 2:20 PM, Czerminski, Ryszard wrote: I would appreciate any advice about how to add rectangular grid to the plot generated by levelplot() similarly as grid() function adds grid to the regular plot() ?panel.grid Best regards, Ryszard -- David Winsemius, MD West

### Re: [R] How do I read multiple rows of different lengths?

and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

### Re: [R] How do I read multiple rows of different lengths?

preprocess with readLines if you want an exact match of width or maybe take a guess with ... , colClasses = rep(numeric, 250) to create a sufficiently wide set of columns to hold everything. -- David Winsemius, MD West Hartford, CT __ R

### Re: [R] Fwd: as.list

-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT

### Re: [R] how to add rectangular grid to the plot generated by levelplot()?

plot() ?panel.grid Best regards, Ryszard -- David Winsemius, MD West Hartford, CT David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

### Re: [R] bugs when using R-2.12.0 Please HELP!

not a member). A word which has less implication of software deficiencies would be problems. This is really about BioConductor and it has its own mailing list. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https

### Re: [R] runtime on ising model

and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

### Re: [R] runtime on ising model

On Oct 28, 2010, at 12:20 PM, David Winsemius wrote: On Oct 28, 2010, at 11:52 AM, Michael D wrote: Mike, I'm not sure what you mean about removing foo but I think the method is sound in diagnosing a program issue and the results speak for themselves. I did invert my if statement

### Re: [R] wait for graph to finish plotting

: Could you please paste the exact code you are using? (the one with the pdf and dev.off, outside the loop ) Contact On Sat, Oct 23, 2010 at 4:53 PM, n...@aleblanc.cotse.net wrote: David Winsemius dwinsem...@comcast.net writes: On Oct 22, 2010, at 12:50 PM, n

### Re: [R] Clustering

package. Can anyone suggest a more appropriate package to use for such a large dataset? David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

### Re: [R] Clustering

of the same event, hence the time conditon. Hope this clarifies the problem. It clarifies it to the extent that it show how much more you will need to further clarify. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list

### Re: [R] strftime vs strptime ??

. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

### Re: [R] Clustering

. Those would be character values and, if you want to do calculations, would obviously need to be coerced to numeric. Thanks, Doug -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

### Re: [R] Clustering

On Oct 29, 2010, at 12:08 PM, David Winsemius wrote: On Oct 29, 2010, at 11:37 AM, dpender wrote: Apologies for being vague, The structure of the output is as follows: Still no code? $ cluster1 : Named num [1:131] 3.05 2.71 3.26 2.91 2.88 3.11 3.21 -1 2.97 3.39 ... ..- attr

### Re: [R] Differenciate numbers from reference for rows

-contained, reproducible code. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

### Re: [R] Clustering

On Oct 30, 2010, at 7:49 AM, dpender wrote: David Winsemius wrote: On Oct 29, 2010, at 12:08 PM, David Winsemius wrote: On Oct 29, 2010, at 11:37 AM, dpender wrote: Apologies for being vague, The structure of the output is as follows: Still no code? I am using the Clusters

### Re: [R] Differenciate numbers from reference for rows

I am missing) these results are concatenated in a string, and then evaluated, a step which I do get. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

### Re: [R] solving equation consisting of factorials

David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

### Re: [R] Using names function

month: max(diff(msales)) I get the numeric value displayed when I execute the above function, but not the month. cumsum(msales) displays the results under each month, but not max. How can I get both the month and value displayed? diff(msales)[which.max(diff(msales))] Mar 450 -- David

### Re: [R] for loop

example applying split to the builtin airquality dataframe. The plyr package also provides functions on dataframes. -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

### Re: [R] Confidence interval for response variable in mixed effects models

http://finzi.psych.upenn.edu/R/Rhelp02/archive/76742.html http://finzi.psych.upenn.edu/R/Rhelp02/archive/81237.html http://finzi.psych.upenn.edu/R/Rhelp02/archive/82567.html -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing

### Re: [R] for loop

the results creating a list of plots. -- DAvid. Thanks,m -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Saturday, October 30, 2010 8:24 PM To: Matevž Pavlič Cc: r-help@r-project.org Subject: Re: [R] for loop On Oct 30, 2010, at 2:07 PM, Matevž Pavlič

### Re: [R] two group cox model

. It seems to create a solid and dashed line in the legend when I use it. ?par # 1=solid (default), 2=dashed, Many thanks in advance. Cheba David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https