similar to: Plotting ordered nominal data

Hi I am a newbie to R but have tried a number of ways in R to do this and can't find a good solution. (I could do it out of R in perl or awk but would like to know how to do this in R). I have a large data frame 49 variables and 7000 observations however for simplicity I can express it in the following data frame Base, Image, LVEF, ES_Time A, 1, 4.32, 0.89 A, 2, 4.98, 0.67 A, 3, 3.7, 0.5

I cannot seem to get this to work? I have weighted data points that I want to fit an sd-ellipse to, but the ellipse is not right? I have attached a .tiff of what I get. I'm wondering if it's not looking at my actual data, but taking the data from the example? I ask this because I don't know where "Title" comes from in the graph? Also, I want to put calcentre=FALSE, but if I

Hi, I've seen posts (but no answers) about clearing a file lock on a Samba share. Presumably this kind of thing happens when the program abnormally terminates and doesn't close its open files. I've got a file that smbstatus reports as: 1216 DENY_WRITE 0x20089 RDONLY NONE /shares/share1/Files/december.xls Fri Jan 12 12:51:30 2007 Trying to open this in Excel

Hi I have a long data set on which I want to do Bland-Altman style plots for each rhythm type Using ggplot2, when I use geom_hline with facet_grid I get an extra set of empty panels. I can't get it to do it with the "Diamonds" data supplied with the package so here is a (much abbreviated) example: > lvexs cvd_basestudy ecd_rhythm fixed_time variable_time 1 CBP05J02

Dear Community, I have been using two types of survival programs to analyse a data set. The first one is an R function called aftreg. The second one an STATA function called streg. Both of them include the same analyisis with a weibull distribution. Yet, results are very different. Shouldn't the results be the same? Kind regards, J -- View this message in context:

  Hi r-users,   I would like to use gev and my data (annual rainfall ) is as follows:   > head(dat,20) A B C D E F G H I J 1 45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0 2 50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4 58.3 3 41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5 45.6 4 50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7 38.3 5 39.3 30.6 46.9 23.8 25.8

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 by base_study. My dataframe (800 different base_studies) looks like > afvtprelvefs basestudy quartile ef ef_std entropy CBP0908020 1 21.6 0.53 3.27

Hi - I am not R expert and I would appreciate your time if you can help me about my xyplot question. I would like to add text (p-value) in a 4 panels xyplot. I thought panel = function{} should work but I am not sure where I did it wrong. The error message from the following code is "Argument subscripts is missing with no default values" xyplot(GLG ~ PD | factor(TRT) , groups =

Hello all! I've been struggling with is for many hours today; I'm close to getting what I want, but not close enough... I have a dataframe consisting of two date-time columns followed by two numeric columns. what I need is the max value (in the first numeric column) based on the 2nd date-time column, which is essentially a factor. But, I want the result to provide both date-time values

Hi Jeff, Thank you for the suggestions -- I appreciate your help. Unfortunately, the result2 has two problems... (1) there are now 3 date columns (it looks like 2 cols are merged into 1 col) (2) the output rows should not have any of the basistime dates repeated (maybe I misstated the problem); I need the max fcst value by basistime, but also list the date value for that row; for example:

Hi With ggplot2 I can very easily create beautiful histograms but I would like to put two histograms on the same plot. The histograms may be over-lapping. When they are overlapping the bars are shown on top of each other (so that the overall height is the sum of the two). Is there any way to get them to display overlapping (with smaller value in front, larger value behind) so that the overall

Hi This is probably more of a statistics question than a specific R question, although I will be using R and need to know how to solve the problem in R. I have several sets of data (ejection fraction measurements) taken in various ways from the same set of (~400) patients (so it is paired data). For each individual measurement I can make an estimate of the percentage uncertainty in the

Dear all- I have a dataset (see a sample below - but the whole dataset is June 2005 - June 2008). The "LST" format is "YYMMDDHHmm" and I would like to get the hourly average of the "mph" for the summer months (spanning all years). I have been trying to use "aggregate" but am not having much success at all! any thoughts would be greatly appreciated.

Thanks for the dput... #### reproducible example of split-apply-combine ### dta <- structure(list(date = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L), .Label = c("2012-01-25 18:00:00",

On Thu, 26 Oct 2017, Thomas Adams wrote: > Hi Jeff, > > Thank you for the suggestions -- I appreciate your help. Unfortunately, the > result2 has two problems... > > (1) there are now 3 date columns (it looks like 2 cols are merged into 1 > col) No, there are two date columns. Result2 includes the grouping value as a row name (pulled from the names of the dta2list items

I have some data which is censored and I want to determine the median. Its actually cost data for a cohort of patients, many of whom are still on treatment and so are censored. I can do the same sort of analysis for a survival curve and get the median survival... ...but can I just use the survival curve functions to plot an X axis that is $ rather than date? If not is there some other way to

Hi all. I got blocked with the first post and and the message now has been posted multiple times! Sorry about that :( Kind Regards, Keith Begin forwarded message: > From: Mail Delivery System <MAILER-DAEMON at nhs-pd1e-esg106.ad1.nhs.net> > Subject: Undelivered Mail Returned to Sender > Date: 10 July 2015 09:30:06 BST > To: <Keith.Roberts at ecric.nhs.uk> > >

This is the same as https://github.com/networkupstools/nut/issues/254 , right? (It's on my TODO list. DDL updates take a little longer than HCL updates.) Other comments below. > On Dec 29, 2015, at 1:55 PM, Thiago Wiezbicki - NHS <thiago.wiezbicki at nhs.com.br> wrote: > > Device Manufacturer: NHS > Device Name: Laser Senoidal 5000VA > > upsc output when the UPS is

I wasn't clear. battery.runtime.low is not a counter, but a value estimated in seconds that is decreasing conforming the battery is discharging. It's mean: How much time (in seconds) will be necessary until the battery achieve the depleted status ( http://www.networkupstools.org/protocols/sec.html#Estimated_Minutes ) . If the UPS don't have a load, this value will be higher than the

Hi, Taking this to a new thread. Thoughts on getting a new list started up ? Should it be centos-sysadmin or centos-infra ? Are we going to then restrict it to admin/infra related chatter ? in which case, does the eyeball density on this list reduce for that sort of content ? How about the politics and news stuff that gets posted to this list. And how would those things be addressed with