[R] ggplot2 - expand range?

Dear UserRs, I am trying to use systematically ggplot2 for most of my plots, but I am fighting some lack of documentation, which I try to overcome. I want to build a scatterplot where the axes cross exactly at (0,0). I tried using scale_y_continuous(limits=c(0,10)), but I always get an extra

[R] Start plot really at baselines x=0, y=0

Dear R Plots without par arguments do start not at zero (means, the box around is somwhere outside the specified plotrange). How to start really from zero, pe. basline y=0? every standart par works like this: x-seq(1,10,0.1) y-sin(x) plot(sin(x)) Thanks for help Marc -- Ist Ihr Browser

[R] Some problem in opening connection with .dat extention file in matrix(scan) function of R 2.5

Dear helpers please provide me some helpful answer to my problem while I m trying to run a program .I m attaching both the program and the data to which I have to obtain my estimation results. Motives.dat is the data file, and OBTfile4.3 is the complete code of program. by Running this //

Re: [R] Updating packages

On Thu, 1 Nov 2007, John Field wrote: I've just updated from R-2.5.1 to R-2.6.0, under Windows XP. Is there a simple way to reload into R-2.6.0 all the packages I had loaded under R-2.5.1? Do you mean 'installed'? If so, this is covered in the rw-FAQ. -- Brian D. Ripley,

Re: [R] F distribution from lme()?

Dear Andreas, try: anova(incub.lme2, type=marginal) Read more about marginal vs. sequential tests in the help page ?anova.lme. Christian __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] computing the mean of a few variables

hello, I have a simple question: I want to compute the average of 6 variables, my data looks like this: t0t1 t2t3 t4 t5 1 32 34 36 40 38 40 2 35 37 38 36 33 39 . . . I want the mean of t0-t5 for each row (observation). thank you, sigalit. [[alternative

Re: [R] computing the mean of a few variables

On Thu, 2007-11-01 at 10:34 +0200, sigalit mangut-leiba wrote: hello, I have a simple question: I want to compute the average of 6 variables, my data looks like this: t0t1 t2t3 t4 t5 1 32 34 36 40 38 40 2 35 37 38 36 33 39 . . . I want the mean of

[R] Matrix package problem: dsyMatrix %*% vector gives error

Dear UseRs, here is an example scenario presenting my problem: Multiplying a dsyMatrix with a numeric vector results in an error (unfortunately in German due to my locale): (M1 - Matrix( c( 1, 2, 2, 2, 1, 2, 2, 2, 1), nrow = 3)) 3 x 3 Matrix of class dsyMatrix [,1] [,2] [,3] [1,]1

Re: [R] computing the mean of a few variables

sigalit mangut-leiba wrote: I want to compute the average of 6 variables, my data looks like this: t0t1 t2t3 t4 t5 1 32 34 36 40 38 40 2 35 37 38 36 33 39 . . . I want the mean of t0-t5 for each row (observation). If your data are available in a

Re: [R] Find A, given B where B=A'A

On Wednesday 31 October 2007 22:45, Peter Dalgaard wrote: Michael Gormley wrote: Thanks for your help, all those who submitted responses. I do not need a specific matrix A, any solution will do. With this said, is it possible to specify the dimensions of the A matrix in the decompostion?

Re: [R] Matrix package problem: dsyMatrix %*% vector gives error

As the posting guide says: If you are using an old version of R and think it does not work properly. upgrade to the latest version and try that, before posting. This works in current R and current Matrix (0.999375-3) It also says If the question relates to a contributed package , e.g.,

Re: [R] computing the mean of a few variables

Thank you! Sigalit. On 11/1/07, Tobias Verbeke [EMAIL PROTECTED] wrote: sigalit mangut-leiba wrote: I want to compute the average of 6 variables, my data looks like this: t0t1 t2t3 t4 t5 1 32 34 36 40 38 40 2 35 37 38 36 33 39 . . . I

Re: [R] R GUI for Linux

This is a reported bug in Suse10.3. It is described in https://bugzilla.novell.com/show_bug.cgi?id=330991, where it is listed as a major, but unresolved, bug. It has been like that for at least a week . . . :( Ted. Dr EA Catchpole Visiting Fellow [EMAIL PROTECTED] and Univ of Kent

Re: [R] drawing direction fields

Jong-Hoon: If a direction field is the equivalent of a phase plane diagram, then you can use a function written by Daniel Kaplan, Macalester College. You can find the function in his online notes for an applied calculus course: http://www.macalester.edu/~kaplan/math135/index.html Jong-Hoon

[R] Problem with compiling 64bit R(2.5.1) under HP-UX(ia64)

Hi there, We are trying to compile a 64bit version of R (2.5.1) on HP-UX (B.11.23 U ia64), but are running into some problems. This is our configure step: ../configure --prefix=/rnd/homes/lfan/R251 --enable-R-shlib CC=cc CFLAGS=+z +DD64 CXX=aCC CXXFLAGS=-b -lxnet +z +DD64 FC=f90 FCFLAGS=+DD64

[R] Zelig and the blogit model

Hi Folks, According to the PDF file blogit.pdf in the Zelig documentation: Use the bivariate logistic regression model [blogit] if you have two binary dependent variables (Y1,Y2), and and wish to model them jointly as a function of some explanatory variables. Each pair of dependent

Re: [R] Some problem in opening connection with .dat extention file in matrix(scan) function of R 2.5

I would assume that the file that you are trying to read is not in the current directory your program is in. Do 'list.files()' to see what is in the current directory. Also use 'setwd' to get you to the directory you should be in. On 11/1/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Dear

Re: [R] F distribution from lme()?

Using anova with the default setting generates a sequential analysis of variance table. You can see this by noting that if you change the order of terms in the model, it gives you a different result: incub.lme2 - lme(egg.temp ~ kjday + treat, random = ~1|ind, data = incub.df) incub.lme3 -

[R] structure data

Hello, I have a quite simple question, I guess. I have a data frame and I would like to process the data in several ways. The processing is dependent on grouping (factors) and the parameter itself. Meaning for parameter1 in the df, I would like to calculate the mean and sd when grouping by

[R] ploting a comparison of two scores, including the labels in the plot

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*

[R] A question about lme object

I have a question about the lme function in R. My question is: After I got the object from function lme, why the numIter value of the object is always NULL? Following is my code: jjww-lme(y~x*zz,data=simul,random=~x|group, control=lmeControl(returnObject=TRUE)) attributes(jjww) jjww$numIter

Re: [R] RWeka and naiveBayes

Sancar, I'm trying to use RWeka to use a NaiveBayes Classifier(the Weka version). However it crashes whenever there is a NA in the class Gender It is difficult to tell what the problem is without the data (or a reproducible toy example), but it seems like you are trying to remove NA values

Re: [R] seg fault with randomForest ( ... , xtest )

Hi Clayton, Using formula like that, IMHO, is a genuinely bad idea. You subset Y on the LHS of the formula, but what about the RHS? Also, I've never tried to use xtest via the formula interface, and that is never the intended usage. I'll look into more checking in the code to prevent such

[R] Question about lme object

Hi, I have a question about the lme function in R. My question is: After I got the object from function lme, why the numIter value of the object is always NULL? Following is my code: jjww-lme(y~x*zz,data=simul,random=~x|group, control=lmeControl(returnObject=TRUE)) attributes(jjww)

Re: [R] Don't understand write.csv default: Why is column name forcolumn of row names blank?

Gabor Grothendieck [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Would this be good enough: # example using builtin BOD data frame write.csv(cbind(Z = row.names(BOD), BOD), row.names = FALSE) Thanks, Gabor, I like your solution better than mine. But, I still don't understand

Re: [R] thicker axis in levelplot

Thank you for your reply - I'm afriad that just doesn't work, any other suggestions or is there something I need to do WITH this command to make it work? I really am stuck on this one and would appreciate any help offered, Jenny From: jim holtman [EMAIL PROTECTED] To: Jenny Barnes [EMAIL

[R] Calculate aggregate differences

Hi everyone I am trying to summarize a table with yield estimates of a forest plantation. For that I have four blocks and four treatments measured over a period of 10 years (every year). In each plot trees are measured (diameters and heights). With aggregate function I can calculate the average

[R] cluster analysis

Hi Sir How can we select the optimum number of clusters? Best Regards -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] ctest package

Hi Sir I did find ctest package from CRAN. Where it can be found? Regards -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Calculate aggregate differences

Thanks Andrew, but what you gave me can actually be done simpler like my_agg$vol.diff - diff(my_agg$Vol) or my_agg$vol.diff -c(NA, diff(my_agg$Vol)) # for a list with the same #length as the aggregated mean list However, what I need is to be able to

[R] problem with log axis

Hello, if I do: x - c(0.5,1,3,6,10,20,40) y - 10-log(x)+rnorm(7,0,0.05) r1 - lm(y ~ log(x)) plot(log(x),y) abline(r1) # I get a nice plot with the regression line almost over the points. but: plot(x,y,log=x) abline(r1) gives me exactly the same plot for the points but the regression

[R] nlme installation/ loading problem

After installation of the latest version of R (2.6.0) I have a problem with loading the nlme library. When I use RGui under Windows XP, and try to install the nlme package, I get a warning that installation was unsuccessful due to a previously installed version of nlme. But the nlme library does

Re: [R] problem with log axis

R Heberto Ghezzo, Dr wrote: Hello, if I do: x - c(0.5,1,3,6,10,20,40) y - 10-log(x)+rnorm(7,0,0.05) r1 - lm(y ~ log(x)) plot(log(x),y) abline(r1) # I get a nice plot with the regression line almost over the points. but: plot(x,y,log=x) abline(r1) gives me exactly

Re: [R] textplot() in gplots causes problems (0x9)

Dear Gregory, How can I avoid using tab character when all I want to do is to print a model summary on my pdf device using textplot()? How do I set the font size? If you mean using cex inside textplot, then it does not work. Whether cex is 1 or 0.2 I get the same result, exemplified here: Call:

Re: [R] problem with log axis

Well, here are two attempts that I would have bet on to work, but don't: #Doesn't seems to show up any line at all: abline(a=as.numeric(r1$coefficients[(Intercept)]), b=as.numeric(r1$coefficients[log(x)])) #Line doesn't match points: abline(r1, untf=TRUE) So much for furthering knowledge and this

Re: [R] structure data

Have a look at ?aggregate or the doBy library. They may be what you want. --- Antje [EMAIL PROTECTED] wrote: Hello, I have a quite simple question, I guess. I have a data frame and I would like to process the data in several ways. The processing is dependent on grouping (factors) and

Re: [R] textplot() in gplots causes problems (0x9)

Dear Gregory and members of R-Help I found what was causing the problem, it was the fact that my R Console window is minimized by me horizontally, and therefore summary output does not fit and gets split, and function capture.output() copies that into PDF! I maximized R Console window so that

[R] [R-pkgs] SASxport v. 1.2.0

SASxport Version 1.2.0 is now available --- The SASxport package provides R with full support for reading and writing SAS xport format files. Version 1.2.0 corrects a critical issues with storage of negative numbers, as well as adding additional improvements

Re: [R] ploting a comparison of two scores, including the labels in the plot

I gave it a try with conventional plot and it does not look easy to get a good result. x - YRKE SAMHM INDM Ambassadör 8.32 7.2771 Läkare 8.15 8.1029 Domare 8.14 7.5965 Professor 8.13 7.5618 Advokat7.95

Re: [R] thicker axis in levelplot

On 11/1/07, Jenny Barnes [EMAIL PROTECTED] wrote: Thank you for your reply - I'm afriad that just doesn't work, any other suggestions or is there something I need to do WITH this command to make it work? I really am stuck on this one and would appreciate any help offered, Does this work for

[R] Rcmdr Plugin and menus.txt

Hi all, I'm trying to write a plugin for R Commander, following the model of the TeachingDemos plugin. I am struggling trying to even add items by editing the menus.txt file. I would welcome any help from anyone who has messed with it. Essentially the problem I am having is: I

Re: [R] pdf() device uses fonts to represent points - data alteration?

Hello, I had the same problem. When opening PDFs with a recent developer version of inkscape all circles were replaced by the letter q, see a screenshoot of the imported figure: http://www.simecol.de/figs/R_pdf_inkscape.png I spent at least two hours trying different development versions of

Re: [R] loops sampling

Hi Garth, Your code is really confusing! You should start by reading the help file on the for() function and understanding what it does: ?for Your line for(i in 1:nboot){ } is simply starting a loop around the variable 'i', which will change values following the sequence 1:nboot. It seems

Re: [R] pdf() device uses fonts to represent points - data alteration?

On 2007-November-01 , at 20:18 , Thomas Petzoldt wrote: I had the same problem. When opening PDFs with a recent developer version of inkscape all circles were replaced by the letter q, see a screenshoot of the imported figure: http://www.simecol.de/figs/R_pdf_inkscape.png I spent at

[R] help

My name is enrique diaz of chile and not if it is the correct place, but I would like to know if implemented this one the estimation of the parameter in the transformation of aranda-ordaz, thank you [[alternative HTML version deleted]] __

Re: [R] pdf() device uses fonts to represent points - data alteration?

On Thu, 2007-11-01 at 20:18 +0100, Thomas Petzoldt wrote: Hello, I had the same problem. When opening PDFs with a recent developer version of inkscape all circles were replaced by the letter q, see a screenshoot of the imported figure: http://www.simecol.de/figs/R_pdf_inkscape.png For

Re: [R] Lattice: par.ylab.text alpha ?

wohoo! A bug. I know I wasn't the first to find it, but I had a hunch this time that it was going to be a bug and not just my poor understanding. I wasn't really thinking as just ylab= works fine too. I think I was hoping to squeeze it (the ylab) closer to the axis w/o axis labels. Now

[R] Reading R-help Digests with Mozilla Thunderbird

This is somewhat off-topic but I think that an answer may help other users of R-help. Thunderbird trys to help in the display of messages by greying out quoted text. However when reading r-help in digest form it gets thoroughly confused and usually ends up greying out the fresh text in a

Re: [R] pdf() device uses fonts to represent points - data alteration?

Hi, this was my former suggestion but as the RSvgDevice doc says: This driver currently does not have any font metric information, so the use of plotmath is not supported. Maybe, the future (or already recent) alternative will be: library(Cairo) that can among others also create svg and pdf.

Re: [R] Some problem in opening connection with .dat extention file in matrix(scan) function of R 2.5

The error message is telling you that R cannot find your file. Is your 'motives_pc.dat' file in your R working directory? If not, you have to give a complete path to the scan() function. Julian [EMAIL PROTECTED] wrote: Dear helpers please provide me some helpful answer to my problem while

Re: [R] subsetting problem with multiple criteria: Works in some but not all cases.

I think the problem is with your use of == instead of %in%; try matching - subset(mydata[,c(j+1,j+7)], mydata[,j+1] %in% lone.word) On 11/1/07, John Kane [EMAIL PROTECTED] wrote: I am trying to compare some word lists which have an associate set of numbers. I want to compare word list aa

[R] Need help on an error, when looding library(uroot)

Hi When I'm looding the library uroot something happens with the basic R code that handles numeric classes. Its looks like a bug but I cant isolate the problem. So hopfully one of you smart guys out there can give me some help. This is how the problem looks. Take for example the vector TEST -

Re: [R] Need help on an error, when looding library(uroot)

On 2/11/2007, at 9:39 AM, Tom.O wrote: Hi When I'm looding the library uroot something happens with the basic R code that handles numeric classes. Its looks like a bug but I cant isolate the problem. So hopfully one of you smart guys out there can give me some help. This is how

Re: [R] Need help on an error, when looding library(uroot)

ok, I should have figured that out by myself that you might need some info about my platform. But here it is. Im on my PC running Vista and I have latest version of R. version _ platform i386-pc-mingw32 arch i386

Re: [R] Need help on an error, when looding library(uroot)

I noticed, upon further digging around, that the help for the uroot package lists a function ``as.numeric'' in its index. However in my version of the package no such function appears to exist. If I do an ls() of the appropriate search list item, the only ``as.'' item that shows up is

Re: [R] Calculate aggregate differences

Hi Cristian, yes, indeed. What I'm not sure about is how doing this by groups will give you a different result, if the object is already sorted according to the group structure. As you are working with differences, each group 'loses' an observation, which will be the observation that is

Re: [R] ploting a comparison of two scores, including the labels in the plot

On Thu, Nov 01, 2007 at 02:52:08PM -0400, John Kane wrote: I gave it a try with conventional plot and it does not look easy to get a good result. Thanks alot John Kane! While I see what you mean, I think your solution does a good job and provides a basis for me to work on. If someone would

[R] simple averaging question?

Hi all, Suppose I have a column vector of 600 measurements taken in 1s intervals. What I want is a new vector with the averages for each min (so there would be 10 entries). Is there an efficient way to do this? I’ve been doing it with a ‘for’ loop but something tells me there is a

[R] missing data

Hi There, Will somebody know if there is a function in R which can compute the proximity from an incomplete data matrix? Or any other software which can do this? Thank you. Van __ R-help@r-project.org mailing list

Re: [R] Zelig and the blogit model

It seems no-one was interested in my Zelig query (part reproduced below)! Anyone out there using Zelig? I'm most obliged to Thomas Yee for his response to my query about vglm() in VGAM, confirming what I had begun to suspect from my experiments that vgam() makes separate independent fits to the

Re: [R] simple averaging question?

Jeff Miller wrote: Hi all, Suppose I have a column vector of 600 measurements taken in 1s intervals. What I want is a new vector with the averages for each min (so there would be 10 entries). Is there an efficient way to do this? I’ve been doing it with a ‘for’ loop but something tells

[R] unable to install package ff

Hi all, I've had one of my most miserable R weeks in memory. I'm trying to deal with huge datasets (1GB each) but am running up against those pesky memory limits. The libraries filehash and g.data are not very suitable for what I need. I haven't gotten into the sql thing yet. Most recently I've

Re: [R] simple averaging question?

Moshe, Ted, Chuck, All THREE approaches worked fine for me. Thank you! Sincerely, Jeff Miller Internal Virus Database is out-of-date. Checked by AVG Free Edition. 3:09 PM ---BeginMessage--- Jeff Miller wrote: Hi all, Suppose I have a column vector of 600 measurements

[R] Mplus and R

Has anyone successfully linked Mplus and R and have code on how to do it? It would be great to use R to simulate datasets and then ship to Mplus for SEM analysis. I know one researcher who does this quite regularly but I can’t get ahold of him right now. Sincerely, Jeff Miller

[R] R and Debian Linux

Hi Folks, Lately I've been testing out recent versions of various Linux distributions, preparatory to choosing which one I want to upgrade to. One of these is debian 4.0 (Etch). I can install R on it as a debian package (it comes through as R version 2.4.0 Patched (2006-11-25 r39997). Many of

[R] Shading area under PDF of t -distribution

Dear R experts, I have plot the PDF of t distribution with df = 74. curve(dt(x,df=74),from=-4, to=4) how can I shade the area under curve (for example, col=red) from t=+- 1.996? Thank you. Regards, CH -- CH Chan Research Assistant - KWH http://www.macgrass.com

[R] density calculation

Hello, I have another, more statistic related question. I have several different measurements as an output of an image analysis process (e.g. size, intensity, ... of objects of interests). Now, I would like to quantify the difference for each parameter measured in a special condition comparing to

Re: [R] R and Debian Linux

Ted, On 2 November 2007 at 00:21, (Ted Harding) wrote: | Hi Folks, | Lately I've been testing out recent versions of various | Linux distributions, preparatory to choosing which one | I want to upgrade to. | | One of these is debian 4.0 (Etch). | | I can install R on it as a debian package (it

Re: [R] Strange behavior with time-series x-axis

Allen, I can reproduce the problems you described, but I can't find where the problem is (in axis.POSIXct). Anyway, here is an alternative I wrote a while ago: # for base graphics timeAxis - function(side, at=NULL, labels=TRUE, ..., tz=GMT) { range - par(usr)[if (side%%2) 1:2 else 3:4]

Re: [R] unable to install package ff

Hi Benilton all, I'm fairly clueless when it comes to compiling etc... never done it. Thanks for your response. Here's the info you requested: ~/Documents/Teaching/R matthewkeller$ make make: *** No targets specified and no makefile found. Stop. ~/Documents/Teaching/R matthewkeller$ gcc

[R] Die Chance des Jahres 2007 für [EMAIL PROTECTED] t.math.ethz.ch - Jetzt zuschlagen und endlich nebenb ei Geld verdienen

Hallo [EMAIL PROTECTED], sehen Sie sich jetzt Ihre Chance des Jahres 2007 an. Sie wurden ausgewaehlt fuer diese einmalige Gelegenheit! Setzen sie jetzt zum Schlussspurt an und kassieren Sie für Ihre Arbeit endlich das Geld, was Sie verdient haben. Verdoppeln und verdreifachen, ach was,

Re: [R] Rcmdr Plugin and menus.txt

Thanks John, no hurry. I tried to debug it a little bit, without much success. Two things I discovered: 1) The exact same lines work just fine if I place them in the Rcmdr- menus.txt file in the Rcmdr package (and add the testOne function also.) 2) While trying to do debug(Commander), the