similar to: Trying to pass arrays as arguments to a function

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

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]]

R-gurus... I've got a 5 column dataframe where I'd like to plot each ID's "b" against "c" with "b" in ascending order (within the same ID). How do I sort "b" so that the other variables are altered equally? ID a b c d 101 1 240 26.7 21.85 101 2 335 21.8 21.85 101 3 1387 26.6 21.85 101 4 877 24.1 21.85

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

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]]

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,

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

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hello, I'm a novice user trying to figure out how to retain NA aggregate values. For example, given a data frame with data for 3 of the 4 possible factor colors("orange" is omitted from the data frame), I want to calculate the average height by color, but I'd like to retain the knowledge that "orange" is a possible factor,

hi, say i have a function f and i'd like to to call it like this: 1) f(list(a=...,b=...)) but i can't do it, because f is defined as: 2) f<-function(a=NULL,b=NULL){...} is there a way that i can approximate 1), such as mapping list(a=,...b=...) to "list(a=,...b=...)" and then replacing "list" by "f", and then evaluating the expression? thanks.

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]]

May be you are also not familiar with statistic. the solution of min \sum_{i=1}^{n} (y_{i}-b)^{2} is the mean. So the solution is b<-mean(y) Justin BEM Elève Ingénieur Statisticien Economiste BP 294 Yaoundé. Tél (00237)9597295. ----- Message d'origine ---- De : Kåre Edvardsen <ked@nilu.no> À : R-help <r-help@stat.math.ethz.ch> Envoyé le : Vendredi, 8 Décembre 2006,

Hi all, I have roughly fifty dataframes and a dataframe with the names of the fifty dataframes. I want to perform the same set of manipulations on all fifty dataframes, but can't find a way to batch process from a list with the dataframe names using a loop. Is there a way to read the file names from the dataframe with the names and then call the referenced dataframe? This would save me a

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

Dear R-helpers, I've got two functions; callTimes() calls times(), passing it an optional argument (bar) by name (bar=harry). times() then believes it has been passed a name, rather than a value ? but I want the value, not the name. Worse, if I evaluate the name, it is evaluated in the environment times() was defined, not where it is called. How can I call times(), defining its optional

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

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