similar to: Sorting problem

I'm writing an publication in two column format and need to shrink some plots. After increasing the axis labels it does not look nice at all. The y-axis label and tick labels almost touch each other and the x-axis tick labels expand into the plot instead of away from it. Is there a better way than "cex" to control the: 1) font size of axis and tick labels 2) font thickness 3)

Hi all! I've got a datstructure like this: subject week value 1 1 4 2 1 8 3 1 3 4 1 5 1 2 5 2 2 6 3 2 2 4 2 6 1 3 3 2 3 7 3 3 3 4 3 7 I'd like to plot the mean of 'value' against week. Is there a direct way of doing this or do I have to make a new structure with the calculated values and then plot it? All the best! -- ########################################### Kare

I'd like to avoid looping through an array in order to change values in the array as it takes too long. I red from an earlier post it can be done by "do.call" but never got it to work. The Idea is to change the value of "y" according to values in "x". Wherever "x" holds the value 3, the corresponding value in "y" should be set to 1. So I

Dear R-gurus! Is it possible within boxplot to break the y-scale into two (or anything)? I'd like to have a normal linear y-range from 0-10 and the next tick mark starting at, say 60 and continue to 90. The reason is for better visualising the outliers. All the best, Kare [[alternative HTML version deleted]]

Hi r-helpers... Why do I get this strange huge jump of 36524 days when changing "origin" from 1969-01-01 to 1968-12-31. It should still be close to zero! This really messes up my calculations of follow-up times in my analyses. > julian(strptime("010169", format = "%d%m%y"),origin = as.Date("1969-01-01")) > Time difference of -0.04166667 days >

Is there a function providing more descriptive statistics than "summary()"? I'm working with a coxph analyses and would like to have more info on certain numbers. If my call is something like: Call: coxph(formula = Surv(followup, CasesCancer) ~ age + BMI + parity + HRT) I'd like to know: * How many CasesCancer was excluded (not only the total number of excluded due to

Hello r-help! I have data with two kind of ratings on status of 100 occupations. The first kind of rating is on the percieved "objective" status that these occupations have in society at large, and the second kind or rating is on the status that the respondents think that these occuption *should* have. The ratings were originally integer values in the rage 1-9, but in the current data,

Hi all. If I run the simple regression when x is a categorical variable ( x <- factor(x) ): > MyFit <-coxph( Surv(start, stop, event) ~ x ) How can I get the overall p-value on x other than for each dummy variable? > anova(MyFit) does NOT provide that information as previously suggested on the list. All the best, Kare [[alternative HTML version deleted]]

Dear R-helpers, I'm for sure not familiar with R, but it seem like a nice sofware tool, so I've decided to try using it. Here is my problem I just can't figure out: I'd like to do least square fit of a straight horizontal (a = 0) line y = ax + b through some data points x = (3,4,5,6,7,8) y = (0.62, 0.99, 0.83, 0.69, 0.76, 0.82) How would i find b????? All the best, Ked

Hello, I have the following data.frame structure(list(Study = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L,

Hi Gang, I've got the latest SVN branch of 1.4 downloaded onto SUSE 11.0. Everything is happy EXCEPT, I can't get fax to be recognized by make menuselect. I tried copying app_rxfax.c and app_txfax.c to the apps directory and starting again from ./configure, but no joy. Any suggestions? Danny Nicholas -------------- next part -------------- An HTML attachment was

I'm trying to add 'gplots' from install.packages("gplots", lib = "/usr/lib/R/library", dependencies = TRUE) on a linux host, but R does not seem to figure out there's a new package installed (yes, I've restarted the R-session). I've tried the default, and different library folders, without success. Installation does not report any warnings or bad exit

I'm trying to devide x into tertiles, but ends up with integer limits even x holds one decimal. The analysis is extremely sensitive to the limits and I like to keep them right. How can that be done? quartiles <- quantcut( x[x >= 0], q=seq(0,1, by=(1/3)) > table(quartiles) quartiles [180,344] (344,448] (448,644] 16467 16476 16452 [[alternative HTML version deleted]]

Hello all, I am not a Samba expert by any means, so please bear with me if I make a fool of myself. We run a small network in this office. There is a Windows domain running off of a Debian Linux server, running Linux 2.4.25. There are mostly Windows 2000 and Windows XP machines, but there are still a few users clinging on to Windows 98 (this is because they run specialized testing software

Hi! First of all: I'm a newbie to both statistics and R, so please be patient with me... I do however, like R because I've been programming (pascal, IDL, perl, C etc) and designing models since -92, but never related to statistics. Ok, here we go: I've got a set of 15 people, all of them observed over 10 weeks (10 analysed blood samples) with - let us kall it the A-value - influenced

Hi, I want to run t.test() for several variables among two groups, and I would like to skip the tedious process of collecting information to assemble a table, but I am not sure if the function I want already exists. Any suggestion would be appreciated. I have a working example, as required by the posting guide: my_swiss = swiss[-1,] my_swiss$facto = rep(1:2,nrow(my_swiss)/2)

Full_Name: K?re Jonsson Version: 2.10 OS: WinXP Submission from: (NULL) (85.225.178.169) I made a brand new installation of R 2.10 in two machines and got the same problem. Machine 1 is Win XP on metal installation machine 2 is virtual in VMWare running a hacky OS named Tiny XP in which any Microsoft windows XP siftware always successfully is installed. Operation 1. Add package RMySQL from

Hi All, I have a small dataset named das (43 cases) in which I am trying to create a binary outcome (1/0) based on the following code: if (das$age>65 && das$bmi>30) {das$danger<-1} else das$danger<-0 I am setting a flag called 'danger' to 1 of the subject is over 65 and has a BMI > 30. I find that my statement evaluates the first record in the data.frame and

Hi How would I analyse time series with - different lengths (i.e. one has 9 entries and the other has 14 entries) - different frequency (i.e. dates are random - no repeated length) - multiple values for the same time entry (e.g. 2009-10-23 below) i.e. my data takes the form: 1st time series 2009-10-07 0.009378 2009-10-19 0.014790 2009-10-23 -0.005946 2009-10-23 0.009096 2009-11-08 0.004189

Hello, I am drawing some data with ggplot2 and would like to reverse the order of the y axis and set a custom range for it. I can do the former but when I add the key `limits` to `scale_y_reverse` I get an error as shows below and, worse, no data shown in the plot. How can I properly set a reverse y scale with a defined range in ggplot2? Thank you ``` A <- rep(c(0, 3, 6, 12, 24), 3) B <-