Re: [R] how to include bar values in a barplot?

Quoting Greg Snow <[EMAIL PROTECTED]>: > My original intent was to get the original posters out of the mode of > thinking they want to match what the spreadsheet does and into thinking > about what message they are trying to get across. To get them (and > possibly others) thinking I made the stat

Re: [R] Changing font in boxplots

ey <[EMAIL PROTECTED]> wrote: > On Wed, 8 Aug 2007, G Iossa, School Biological Sciences wrote: > >> Hi John, >> >> Thanks so much for such a quick reply. >> I have tried to set all to Times font running >> >> par(font.lab=6) (not 4, maybe this is a loca

Re: [R] Changing font in boxplots

gt; boxplot > > Example: > par(font.lab=4) > boxplot(mass ~ family, data=mydata, ylab="mass %", > xlab="family",las=1, cex.axis=1) > > --- "G Iossa, School Biological Sciences" > <[EMAIL PROTECTED]> wrote: > >> Hi all, >> >&g

[R] Changing font in boxplots

Hi all, I am very new to R and this might be a simple question but I have looked everywhere you suggest before writing to you. I am trying to change font type from san-serif to a serif (Times New Romans) on all labels and axis of my boxplot. I have used this function in other plots before, e.g

Re: [R] how to include bar values in a barplot?

On Wednesday 08 August 2007 00:40:56 Donatas G. wrote: > On Tuesday 07 August 2007 22:09:52 Donatas G. wrote: > > How do I include bar values in a barplot (or other R graphics, where this > > could be applicable)? > > > > To make sure I am clear I am attac

Re: [R] how to include bar values in a barplot?

On Tuesday 07 August 2007 22:09:52 Donatas G. wrote: > How do I include bar values in a barplot (or other R graphics, where this > could be applicable)? > > To make sure I am clear I am attaching a barplot created with > OpenOffice.org which has barplot values written on top

Re: [R] how to include bar values in a barplot?

On Tuesday 07 August 2007 22:09:52 Donatas G. wrote: > How do I include bar values in a barplot (or other R graphics, where this > could be applicable)? > > To make sure I am clear I am attaching a barplot created with > OpenOffice.org which has barplot values written on top of each

[R] how to include bar values in a barplot?

How do I include bar values in a barplot (or other R graphics, where this could be applicable)? To make sure I am clear I am attaching a barplot created with OpenOffice.org which has barplot values written on top of each barplot. -- Donatas Glodenis http://dg.lapas.info

Re: [R] KMO sampling adequacy and SPSS -- partial solution

to the output; the function and an example are copied below. You can just copy/paste to the console to try it out, provided you have the MASS package. Cheers, Jay *** G. Jay Kerns, Ph.D. Assistant Professor / Statistics Coordinator Department of

[R] using loops to create multiple images

I have a data.frame with ~100 columns and I need a barplot for each column produced and saved in some directory. I am not sure it is possible - so please help me. this is my loop that does not work... vars <- list (substitute (G01_01), substitute (G01_02), substitute (G01_03), substitute (G01_

Re: [R] Simple table with frequency variable

nrow = li, dimnames = list(i, j)) > for (r in 1:li) > for (s in 1:lj) > A[r, s] <- 10*r + s > A > > HTH > /Johan > > > > G. Draisma skrev: >> Thank you Jim, >> Sorry, that I was not clear enough. >> Each case has a frequency variable N.

Re: [R] Simple table with frequency variable

> row col > 1 1 1 > 2 2 1 > 3 3 1 > 4 4 1 > 5 5 1 > 1 1 2 > 2 2 2 > 3 3 2 > 4 4 2 > 5 5 2 >> x <- which(z==1,arr.ind=T) >> paste(rownames(z)[x[,'row']], colnames(z)[x[,'col']], sep=''

[R] Simple table with frequency variable

Hallo, Im trying to find out how to tabulate frequencies of factors when the data have a frequency variable. e,g: i<-rep(1:5,2) j<-rep(1:2,5) N<-10*i+j table(i,j) gives a table of ones as each combination occurs only once. How does one get a table with the corresponding N's?

[R] quantreg behavior changes for N>1000

Hello again R-experts and novices (like me), This seems like a bug to me - or maybe it's intentional...can anyone confirm? Up to 1000 reps, summary() of a rq object gives different output and subtly different confidence interval estimates. ThanksJeff testx=runif(1200) testy=rnorm(1200, 5)

[R] plotting a summary.rq object in using pkg quantreg

Hello, I am having problems adjusting the plot output from the quantreg package. Anyone know what I'm doing wrong? For example (borrowing from the help files): plot(summary(rq(foodexp~income,tau = 1:49/50,data=engel)), nrow=1, ncol=2,alpha = .4, ols = TRUE, xlab="test") The "alpha=" paramete

[R] how do I draw such a barplot?

;) dettach(tolerancija.data) However, this is a lot of writing. How do I make into a function to be reusable? Donatas On Jul 16, 2007, at 9:06 , Donatas G. wrote: > Hi, > > I cannot figure out how to draw a certain plot: could someone help   > me out? > > I have this data.fram

[R] how to determine/assign a numeric vector to "Y" in the cor.test function for spearman's correlations?

does not have any default value assigned... (excuse my translation) WHAT SHOULD I DO??? HOW CAN I DEFINE "Y" AS A NUMERIC VECTOR??? thanx guys, i just can't express myself as a computer would... G :) -- View this message in context: http://www.nabble.com/how-to-determine-as

[R] how do I draw such a barplot?

Hi, I cannot figure out how to draw a certain plot: could someone help me out? I have this data.frame from a survey my.data that looks like something like this: col1 col2 col3 col4 1 5 5 4 5 2 3 5 3 1 3 2 3 4 5 4 3 1 1 2 5

Re: [R] charset in graphics

On Friday 13 July 2007 17:56:45 Prof Brian Ripley wrote: > On Fri, 13 Jul 2007, Peter Dalgaard wrote: > > Prof Brian Ripley wrote: > >> On Fri, 13 Jul 2007, Peter Dalgaard wrote: > >>> The correct incantation seems to be > >>> > >>> postscript(font="URWHelvetica", encoding="ISOLatin7") > >>> plot(0

[R] charset in graphics

How do I make Lithuanian characters display correctly in R graphics? Instead of the special characters for Lithuanian language I get question marks... I use Ubuntu Feisty, the locale is utf-8 ... Do I need to specify somewhere the locale for R, or - default font for the graphics? -- Donatas

[R] elementary statistics with R (rkward?)

stead... -- Donatas G. __ R-help@stat.math.ethz.ch 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.

Re: [R] Sweave bug? when writing figures / deleting variable in chunk

> On Wed, 27 Jun 2007 21:06:04 -0400 D G Rossiter wrote: > >> I have found a quite strange (to me) behaviour in Sweave. It only >> occurs in the following situation: >> >> 1. define a variable in one chunk >> 2. use it within a subsequent figure-generating ch

[R] R function command on a list

Hello R: I am working with a self-defined function and I wish to subject a list (t) to this function. My list is a list of tables: $rs7609589 2/2 2/4 4/4 2/2 2/4 4/4 89 188 87 89 188 87 $rs3909907 1/1 1/4 4/4 94 178 92 $rs12748004 0/0 1/3 3/3 37 150 177 $rs6695928 2/2 2/4 4/4 35 12

[R] Sweave bug? when writing figures / deleting variable in chunk

letions into a later code chunk. This only has the disadvantage of making more chunks, so now that I know what's happening it's no big deal. But it's driving me crazy... am I missing something? Thanks! D. G. Rossiter Senior University Lecturer Department of Earth Systems Analys

[R] codamenu crash

lve this behaviour Sincerely Pablo -- Pablo G Goicoechea Conde Don Vela 54, 2º Izda 01009 Vitoria-Gasteiz Spain Phone: +34 945 282417 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gui

[R] general question about use of list

Hi, is this list only related to R issues or it has a broader context regarding questions and discussions about statistics. Is there any other email list or forum for that? For example, I have a question regarding variance. It is defined as: variance = sum(sq(Xi-mean)) / (N-1) and I never understo

Re: [R] barchart producing incorrect number of barcharts when columns renamed

xlab='Mean Intensity', main=paste("Mean Cluster Intensity for Errors", clipping), ylab='Group', col=rainbow(length(rownames (errs.Means))), border=rainbow(length(rownames(errs.Means)

[R] barchart producing incorrect number of barcharts when columns renamed

quot;\n"); print(barchart(stops.Means, groups=rownames(stops.Means), xlab='Mean Intensity', main=paste("Mean Cluster Intensity for Stops", clipping), ylab='Group', col=rainbow(3), border=rainbow(3))) # dev.off() } dev.off

Re: [R] help interpreting the output of functions - any sources of information

> > Sarah Thanks for the hint, that will do... On 4/24/07, Y G <[EMAIL PROTECTED]> wrote: > > Hi, > > > > I am looking for documentation, reference guides, etc. that explain the > > output of functions... For example using cor.test(, > method="

[R] help interpreting the output of functions - any sources of information

Hi, I am looking for documentation, reference guides, etc. that explain the output of functions... For example using cor.test(, method="pearson") with Pearson's corr coeff the output is: Pearson's product-moment correlation data: a and b t = 0.2878, df = 14, p-value = 0. alternativ

Re: [R] Power analysis and mixed model

-Original Message- From: Jean G. Orelien [mailto:[EMAIL PROTECTED] Sent: Wednesday, April 04, 2007 11:22 AM To: 'Julien' Subject: RE: [R] Power analysis and mixed model Julien, You could write your Linear Mixed Model as a General Linear Multivariate Model (GLMM) and use

Re: [R] Cannot allocate vector size of... ?

or Applied Social Research University of Oklahoma 2 Partners Place 3100 Monitor, Suite 100 Norman, OK 73072 -Original Message- From: Patrick Burns [mailto:[EMAIL PROTECTED] Sent: Thursday, March 15, 2007 12:31 PM To: Wingfield, Jerad G. Subject: Re: [R] Cannot allocate vector size of...

[R] Cannot allocate vector size of... ?

Hello all, I've been working with R & Fridolin Wild's lsa package a bit over the past few months, but I'm still pretty much a novice. I have a lot of files that I want to use to create a semantic space. When I begin to run the initial textmatrix( ), it runs for about 3-4 hours and eventually g

[R] Calling Optim() from C

Hello: I am sure this question was dealt with several years ago. Is the function vmmin() available from Rmath Standalone? If not is it possible to call optim() or nlm() from Rmath in C. Thank you. Mervyn __ R-help@stat.math.ethz.ch mailing list

Re: [R] Fwd: Re: [friday topic]: what exactly is statistical com

Ted: I agree with what you and Richard have to say on this topic but disagree with your time frame. I would like to believe that rudiments of what we call statistical computing today started long before the mid-seventies. Most probably it goes back to the mid-fifties. There is a book that is

[R] nlm() problem : extra parameters

Hello: Below is a toy logistic regression problem. When I wrote my own code, Newton-Raphson converged in three iterations using both the gradient and the Hessian and the starting values given below. But I can't get nlm() to work! I would much appreciate any help. > x [1] 10.2 7.7 5.1 3.8

[R] lme problem : extra parameters

Hello: Below is a toy logistic regression problem. When I wrote my own code, Newton-Raphson converged in three iterations using both the gradient and the Hessian and the starting values given below. But I can't get nlm() to work! I would much appreciate any help. > x [1] 10.2 7.7 5.1 3.8

Re: [R] reading very large files

that was only 1/4th done after 25 h! Many thanks Juli On 02/02/07, juli g. pausas <[EMAIL PROTECTED]> wrote: > Hi all, > I have a large file (1.8 GB) with 900,000 lines that I would like to read. > Each line is a string characters. Specifically I would like to randomly >

Re: [R] reading very large files

-02-03 at 19:06 +0100, juli g. pausas wrote: > > Thank so much for your help and comments. > > The approach proposed by Jim Holtman was the simplest and fastest. The > > approach by Marc Schwartz also worked (after a very small > > modification). > > > > It is clear

[R] reading very large files

trs<- trs[sample(length(trs), 3000)] And this works OK; however my computer seems not able to handle the 1.8 G file. I thought of an alternative way that not require to read the whole file: sel <- sample(1:90, 3000) for (i in 1:3000) { un <- scan("myfile", what= character

Re: [R] Calculate NAs from known data: how to?

Torleif Markussen Lunde wrote: > In a dataset I have length and age for cod. The age, however, is ony > given for 40-100% of the fish. What I need to do is to fill inn the NAs > in a correct way, so that age has a value for each length. This is to be > done for each sample seperately (there are 324

[R] new to R: don't understand errors

Hello all, I'm brand new to the use of R, and I'm trying to quickly learning the rudiments for a couple of projects here at work. I'm working with the lsa package and trying to generate various semantic spaces. I seem to do well with small collections of clean text files, but now that I am try

[R] [R-pkgs] the IPSUR package

verse team of dedicated contributors. Indeed, without all of their countless hours of effort the IPSUR package would not be possible. Kudos to them for providing excellent software to the R community. Cheers, Jay *** G. Jay Kerns, Ph.D. Dep

[R] Significant integer digits in summary.default

ms non-intuitive to me. Thank you, D. G. Rossiter Senior University Lecturer Department of Earth Systems Analysis International Institute for Geo-Information Science and Earth Observation (ITC) Enschede, The Netherlands __ R-help@stat.math.ethz.ch mai

[R] Antwort: Re: Antwort: Buying more computer for GLM

Prof Brian Ripley wrote: > I would not have expected glm to be more than say 5x slower than lm if CPU > cycles and not memory were the limiting factor. In that case more RAM > might be all you need. The ratio between glm and lm might well be about 5x, but that's still a big difference for us.

Re: [R] Antwort: Buying more computer for GLM

Prof Brian Ripley wrote > Probably not, but you have the ability to profile in R and find out. Thanks. This is certainly something I could check, and I shall do so. > > > Some more comments; > > 1) The Fortran code that underlies glm is that of lm.fit that only makes >use of level-1 BLA

[R] Antwort: Re: Antwort: Buying more computer for GLM

Peter Dalgaard wrote > Is this floating point bound? (When you say 30 factors does that mean > 30 parameters or factors representing a much larger number of groups). > If it is integer bound, I don't think you can do much better than > increase CPU speed and - note - memory bandwidth (look for larg

[R] Antwort: Buying more computer for GLM

Hello, at the moment I am doing quite a lot of regression, especially logistic regression, on 2 or more records with 30 or more factors, using the "step" function to search for the model with the smallest AIC. This takes a lot of time on this 1.8 GHZ Pentium box. Memory does not seem t

[R] question

Hi Mr Plate, I have a little question How to convert a rowname vector of numbers into a real column of the matrix, My problem is that I applied a rowsum function on a matrix. Then I get a matrix in which the names of the columns are the values of the group (numbers) Now I need to make calculati

[R] Running R on a 64 bit processor

Greetings, We recently obtained a new computer in our lab with a Pentium 4 3.86 GHz processor and 4 gb of ram running windows xp with service pack 2. After installing R on this machine, I ran a bit of code and found that the execution time was actually significantly slower than a machine runni

[R] sparse matrix, rnorm, malloc

Hi, I'm Sorry for any cross-posting. I've reviewed the archives and could not find an exact answer to my question below. I'm trying to generate very large sparse matrices (< 1% non-zero entries per row). I have a sparse matrix function below which works well until the row/col count exceeds

[R] difference of means as response?

st, I am also unsure as to whether this makes any statistical sense to attempt. I'd appreciate any advice, thank you. Regards, Douglas G. Scofield, PhDDepartment of Biology [EMAIL PROTECTED]Indiana University off: (812) 856-0115 1001 E. 3r

[R] Multiple ablines

Thanks in advance, I have searched and read the archives extensively and have been unable to find a question (solution) exactly mimicking my problem. I have created a scatterplot simulation of data through which I have drawn a linear regression in this manner: points(y~x, pch="*", col="black")

Re: [R] R2 in mixed model

mixed models where you have conditional independence. The manuscript has not gone through the peer-review process yet. Jean G. Orelien -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Spencer Graves Sent: Friday, March 17, 2006 11:14 PM To: Angelo Colombo

[R] Error on read.table

Hello When trying to read a table from excel, the generated error message is not documented in the "R data import" document: > v2 <- read.table("v2-101-405-excel.txt", header=T,row.names=0) Error in read.table("v2-101-405-excel.txt", header = T, row.names = 0) : attempt to

[R] repeated measures with random effects

like: aov(H ~ B * A %in% B * Time + Error(id) ) where id is the factor coded for the repeated subjects. But my question is how to include the random factor C ? Any help would be appreciated. Thanks Juli -- Juli G. Pausas CEAM & UA http://www.ceam.es/lass/p

[R] GLM/GAM and unobserved heterogeneity

Hello, I'm interested in correcting for and measuring unobserved heterogeneity ("missing variables") using R. In particular, I'm searching for a simple way to measure the amount of unobserved heterogeneity remaining in a series of increasingly complex models (adding additional variabl

Re: [R] sciviews installation

Thanks for the replies 1. about Rcmdr versions Sorry, I did'nt say wich versions because I tried a lot... before asking for help First I tried the version of the site www.sciviews.org (0.9-14), after I tried the Cran version (1.0-2) and at last I got the 1.1-0 version on John Fox's web. After

[R] sciviews installation

Hello 1. a few months ago, I had sciviews working fine with R (rw2001) under windows XP 2. now, upgrading to rw2011, the stuff seems fine (every package installed),but I find a conflict when launching sciviews: - it runs, apparently - but when I try to work ("import/export In: text" for instanc

[R] bs() function of the splines package with intercept=FALSE

Hello, I'm implementing a function using non uniform B-Splines. Looking at the code of the bs() function, I realized that if the intercept was set to FALSE, the behavior of the function was the following (df is the number of degrees of freedom that I believe can be interpreted as the number of con

[R] bs() function of the splines package

Laurent 14 juin 12:02 afficher les options De : "Laurent" <[EMAIL PROTECTED]> - Rechercher les messages de cet auteur Date : Tue, 14 Jun 2005 03:02:37 -0700 Local : Mar 14 juin 2005 12:02 Objet : bs() function of the splines package Répondre | Répondre à l'auteur | Transférer | Impr

RE: [R] Table Help

text we have lol. Thanks for the reply! Ben -Original Message- From: Sean Davis To: McMurtry, Benjamin G. Cc: 'r-help@stat.math.ethz.ch' Sent: 5/24/2005 4:56 PM Subject: Re: [R] Table Help see ?aggregate Sean On May 24, 2005, at 4:39 PM, McMurtry, Benjamin G. wrote:

[R] Table Help

I have a very large table that I want to add some of the certain rows. The table is as follows: Username1 2 Username1 4 Username2 6 Username2 10 Username3 12 Username3 10 Username3 16 Etc The data is sorted by Column 1 (Ie the username) then by columm 2 (The numbers). I do not know how ma

[R] Table Help

Reply to sender Reply to all Reply to folder Forward Move/Copy Delete Read previous item Read next item Get help information on the current window From: McMurtry, Benjamin G. To: '[EMAIL PROTECTED]' Cc: Subject:Table Help Sent: 5/24/2005 4:33 PM

[R] GAM/MGCV: UBRE Score

explained deviance. Is that a reliable statistic? Are there any other statistics that one would use for model selection? Jean G. Orelien Senior Biostatistician *** SciMetrika, LLC 2 Davis Drive RTP, NC 27709 Tel: (919)765

[R] Output a dataframe from R to excel

Hi, I am trying to output an dataframe from R to Excel file. Can anyone tell me how to do it? Thanks a lot. Eg. R dataframe: ABC 121 342 ... [[alternative HTML version deleted]] __ R-help@stat.math.ethz

[R] GAM: Remedial measures

I fitted a GAM model with Poisson distribution to a data with about 200 observations. I noticed that the plot of the residuals versus fitted values show a trend. Residuals tend to be lower for higher fitted values. Because, I'm dealing with count data, I'm thinking that this might be due to overd

[R] animation without intermediate files?

Hello, Does anyone know how to make "movies" in R by making a sequence of plots? I'd like to animate a long trajectory for exploratory purposes only, without creating a bunch of image files and then using another program to string them together. In Splus I would do this using double.buffer()

[R] splot.screen: multiple plots

I have multiple conditions that I would like to plot in a grid. To save space, I don't need the ticks, labels etc on the plots which are not on the outside. I tried split.screen, but if it is clearer to use, it seems pretty rigid since it decides of each screen size apriori while I need more space

[R] GAM: Overfitting

I am analyzing particulate matter data (PM10) on a small data set (147 observations). I fitted a semi-parametric model and am worried about overfitting. How can one check for model fit in GAM? Jean G. Orelien __ [EMAIL PROTECTED] mailing list

[R] GAM: Getting standard errors from the parametric terms in a GAM model

o get standard errors for the parametric terms? Thanks. Jean G. Orelien __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] Starting with R

I am studying R and within the first 3 lines of demo("graphics") I get completely stuck. could you 1. give me the answer so I can move on 2. tell me how I might find the answer. On my mac OSX, the first command is equivalent to get("quartz")() now get("quartz") would be a regular function format

[R] Matrix decomposition: orthogonal complement

R? Best Regards, /Florin -- Florin G. Maican Ph.D. candidate, Department of Economics School of Economics and Commercial Law Gothenburg University, Sweden E-mail: [EMAIL PROTECTED] P.O.Box 640 SE-405 30 Gothenburg, Sweden Phone no +46 31 773 4866 Fax no +46 31

[R] THanks, random/repl. prob. SOLVED.

Thanks for the help, it worked. Mark G. Orr Postdoctoral Research Fellow Dept. of Neuroscience RM 825 Kennedy Center Albert Einstein College of Medicine Bronx, NY 10461 718-430-2610 __ [EMAIL PROTECTED] mailing list

[R] Sample without replacement

d devs (although close to the desired, not right on). Any help would be greatly apprecitated. -Mark Orr Mark G. Orr Postdoctoral Research Fellow Dept. of Neuroscience RM 825 Kennedy Center Albert Einstein College of Medicine Bronx, NY 10461 718-430-2610

RE: [R] D'agostino test

> > Hi, See below, taken pretty much verbatim from Zar (1999) Douglas G. Scofield, PhDDepartment of Biology [EMAIL PROTECTED]Indiana University off: (812) 856-0115 1001 E. 3rd St. fax: (812) 855-6705 Bloom

[R] removing observations

Hello R list. I am sure this is an easy question, please forgive my ignorance of R. I have a panel data set (25 years) and I would like to break it up by year, and run a few tests. How is this done in R? In STATA it is quite simple (keep if year =1970). If someone could give me an example

[R] creating your own glm link function

In sas and matlab, I have been using a modified link function. I use a binomial distribution with the link log(p-.5/1-p) instead of log(p/1-p) which is the usual logit link. This allows the probabilities to go from .5 to 1 instead of 0 to 1. which is typical when probabilities vary between chance

Re: [R] SJAVA error

; (B>To: Prof Brian Ripley <[EMAIL PROTECTED]> (B>CC: $BLpLn(B $B=g;R(B <[EMAIL PROTECTED]>,<[EMAIL PROTECTED]> (B>Subject: Re: [R] SJAVA error (B>Date: Wed, 09 Jun 2004 10:48:15 +0200 (B> (B>Thanks (B>I am using windows NT, R 1.9.0 and SJava 0

[R] SJAVA error

Hi (B (BI'm trying to use SJava and I have troubles. (BI try to run examples from "Calling R from Java" (Bbut,I have an error that (B"fatal error: enable to open the base package" (B (BI heard SJAVA bug, (Bso,could you send me your compiled SJava package with the modified (BREmbed.c b

[R] Expected Mean Squares

I am an applied statistician new to R, and I was wondering if R can calculate expected mean squares for unbalanced mixed models. Any help would be appreciated. Thanks, Bill Moore [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mai

[R] Question on deriv3()

Hello: Why is deriv3() functioning differently in R from that in Splus using library(MASS) ? For example deriv3(~(t1*log(t2)+lgamma(t1)+(1-t1)*log(y)+y/t2),c("t1","t2"),function(y,t1,t2)NULL) complains of lgamma. Mervyn __ [EMAIL PROTECTED] mailing lis

[R] APE: compar.gee( )

Dear all, I don't understand the following behaviour: Running compar.gee (in library ape ) with and without the option 'data', it give me different results Example: Start R load("eiber.RData") ls() [1] "gee.na" "mydata" "mytree" library(ape) # runnig with the option data= mydata co

[R] subset select within a function

Dear all, I'd like to subset a df within a function, and use select for choosing the variable. Something like (simplified example): mydf <- data.frame(a= 0:9, b= 10:19) ttt <- function(vv) { tmpdf <- subset(mydf, select= vv) mean(tmpdf$vv) } ttt(mydf$b) But this is not the correct way. Any he

[R] NA, deleting rows

Dear colleges, I do not understand the following behaviour: aa <- data.frame(a1= 1:10, a2= c(rep(NA, 5), 1:5) ) aa a1 a2 1 1 NA 2 2 NA 3 3 NA 4 4 NA 5 5 NA 6 6 1 7 7 2 8 8 3 9 9 4 10 10 5 aa[!aa$a2==1, ] # removing rows with a2==1 a1 a2 NA NA NA NA.1 NA NA NA.2 NA NA

Re: [R] barplot & plot together

type="o", xlim=c(0.5,9.5), pch=19, ylim=c(0, 140), axes=F, xlab="", ylab="") axis(4); mtext("Number of fires", side=4, line=3, col=1) legend(2, 155, "Number", lty= 1, pch= 19, bty="n", cex=0.8) legend(2, 163, "Area", f

[R] barplot & plot together

the first. Any trick for setting the 2 plots at the same scale? I have unsuccessfully tried: plot(t, data, type="b", xlim=c(0,7)) plot(t, data, type="b", xlim=c(min(t),max(t))) (R 1.8.1, for Windows) Thanks Juli -- Juli G. Pausas Centro de Estudios Ambientales del Mediterr

[R] Evaluating outer, numeric, variables in an lme object.

Hello! I´m working with a dateset from twelve fertilizer experiments (Trial) with a block structure of three replicats (Block). The treatment levels consist of application method and product but only one intensity. The factor TrCode could for example be BC(CAN) for broadcast calcium ammonium nit

[R] the name of a variable in a function

Dear collegues, How can I get the name of a variable (and not the variable) within a function ? For instance, in the following function, I'd like to create a variable in the dataframe df with the same name to the variable passed in var: prova <- function( var ) { df <- as.data.frame(matrix(nr=2

[R] 1.7.1 console unresponsive following "paste into"

Dear R Help, I've just installed 1.7.1 and have found that the console now exhibits odd behavior when I paste commands into it. I'm on a PC running Windows XP Professional. I use Word as my script editor, and then paste commands into the R console using the standard Windows "copy" within Word th

[R] filling a matrix who's entries are a function of the indices?

Dear R list, What's the best way in R to fill a matrix who's entries depend on some function of the indices? I'm currently doing: Q <- matrix(0, k, k) for (A in 1:k) { for (B in 1:k) { Q[A,B] <- my.function(A,B) } } but I wonder if there is a more terse way. Regar

[R] means comparison with seasonal time series?

Dear R list, I have a sequence of weekly observations of number of adults and larvae in various size classes from a butterfly population living in a subtropical area with pronounced wet and dry seasons. Wet and dry seasons are each defined 26 weeks long with fixed start and end dates. The data sp

[R] means comparison with seasonal time series?

Dear R list, I have a sequence of weekly observations of number of adults and larvae in various size classes from a butterfly population living in a subtropical area with pronounced wet and dry seasons. Wet and dry seasons are each defined 26 weeks long with fixed start and end dates. The data sp

[R] mca & contingency tables - error: "All variables must befactors"

Hi, I would like to do a multiple correspondence analysis with the mca function in the MASS library on data that I have as a contingency table (which I've tried converting to a data frame). For example, = > data(HairEyeColor) > hair.df <- as.data.frame(Ha

[R] formula (joint, conditional independence, etc.) - mosaicplots

Hi, Can someone set me straight as to how to write formulas in R to indicate: complete independence [A][B][C] joint independence [AB][C] conditional independence [AC][BC] nway interaction [AB][AC][BC] ? For example, if I have 4 factors: hair colour, eye colour, ag

Re: [R] plots to metafile and x/y ratio

[EMAIL PROTECTED] wrote: On Sat, 15 Mar 2003, juli g. pausas wrote: I have a problem which I'm not sure if it is due to my machine or if I'm doing something wrong. I'm ploting a little map using plot(val, type="l", asp=1) where val is the object with the latitu

[R] plots to metafile and x/y ratio

Dear all, I have a problem which I'm not sure if it is due to my machine or if I'm doing something wrong. I'm ploting a little map using plot(val, type="l", asp=1) where val is the object with the latitude/longitude data. In the screen the figure is perfect, and I can copy to the clipboard and

[R] any european maps

Dear all, I've seen that there are some maps, or at least costlines in R. The oz package is the map of Autralia and in the fields package there is a US map and a world map. This world map allows to select smaller sections such Europe: library(fields) world(xlim=c(-10,18),ylim=c(36,60)) However,

[R] variables in a correlation matrix

Hi, Given a correlation matrix, how can I know which variables have certain correlation values? aa <- matrix( rnorm(25, 1,1), 5,5) colnames(aa) <- c("v1", "v2", "v3", "v4", "v5") aacor <- cor(aa) For instance, which variables are negatively correlated? aacor[aacor<0] # provide the r values,

[R] hist() basic question

Hi, This is a very basic question, but I would like to undestand hist(). I thought that the hist( , freq=FALSE) should provide the relative frequencies (probabilities), and so they should sum 1, however: set.seed(2) ah <- hist(rnorm(100), freq=F) sum(ah$intensities) [1] 2 set.seed(2) bh <- hist(

  1   2   >