Re: [R] Windows binary version of DPpackage

2009-03-13 Thread Dieter Menne
Debabrata Midya commerce.nsw.gov.au> writes: > This has reference to the package “DPpackage”. The binary version is > available on Mac OS, but I am using Windows XP. > > May I request you to assist me in the followings: > How can I prepare Windows binary of DPpackage from source? > >

Re: [R] changing function argument

2009-03-13 Thread Dieter Menne
Thomas Mang fiwi.at> writes: > I wonder if the following is possible in R: > > Suppose a function takes an argument, and wants to modify this argument > so that the change is visible _at the call side_. It would be what is > for example known as pass-by-reference-to-non-const in C++. > > test

Re: [R] Ggplot2: saving a grid with multiple plots

2009-03-12 Thread Dieter Menne
ONKELINX, Thierry inbo.be> writes: > I have managed to create a figure on the screen with multiple plots in > it. Something like the example below. When I save that with ggsave(), > only the last plot gets saved (pPath in the example) instead of the > entire figure. Any suggestions how I can save

Re: [R] nlme: problem with fitting logistic function

2009-03-10 Thread Dieter Menne
Douglas Bates stat.wisc.edu> writes: > > 3) Use lmer in lme4. Your mileage may vary, I could not find a speedup > >   for my problems, but larger problem might give one. > > Did you mean nlmer in the lme4 package? If so, it may be worthwhile > trying the development branch but that is not somet

Re: [R] nlme: problem with fitting logistic function

2009-03-10 Thread Dieter Menne
Douglas Bates stat.wisc.edu> writes: > > 3) Use lmer in lme4. Your mileage may vary, I could not find a speedup > >   for my problems, but larger problem might give one. > > Did you mean nlmer in the lme4 package? If so, it may be worthwhile > trying the development branch but that is not somet

Re: [R] nlme: problem with fitting logistic function

2009-03-09 Thread Dieter Menne
Dieter Menne menne-biomed.de> writes: > 4) Use C for the core function. This is very effective, and there is at least >on example coming with nlme (was it SSlogist?). I correct: it is phenoModel, calling nlme_one_comp_first in nlme.c Dieter ___

Re: [R] nlme: problem with fitting logistic function

2009-03-09 Thread Dieter Menne
jakub kreisinger seznam.cz> writes: > I am trying to analyze growth data on mice. To do this I attempted to fit logistic curve using nlme package. > However, the dataset I use is large (in total ca 20 000 measures on ca 3 000 individuals) with relatively > complicated structure (several explanato

Re: [R] Adding text to both grid and base graphs

2009-03-09 Thread Dieter Menne
Rik Schoemaker zonnet.nl> writes: > Unfortunately that doesn't help because it requires you to know beforehand > what sort of graph you're generating. I want to be able to generate a graph > (irrespective of the type) and then use a common process to label them so I > don't have to think about wh

Re: [R] merge data frames with same column names of differe nt lengths and missing values

2009-03-08 Thread Dieter Menne
Steven Lubitz yahoo.com> writes: > Thank you - this is very helpful. However I realized that with my real data sets (not the example I have here), > I also have different numbers of columns in each data frame. rbind doesn't seem to like this. Here's a > modified example: > > x <- data.frame(item

Re: [R] xyplot() - can you control how the plots are ordered?

2009-03-08 Thread Dieter Menne
AR gmail.com> writes: > I want to control the plots in the output of the xyplot(). It is > easier to explain it through an example: > .. modified example see below > > The output of the xyplot() is sorted by the alphabetical order of > months (First Apr, then Aug so on). I need to output such t

Re: [R] NonLinear Programming in R - QUERY

2009-03-08 Thread Dieter Menne
Ravi Varadhan jhmi.edu> writes: > Try the package "Rdonlp2", which can handle general, nonlinear, equality and inequality constraints for > smooth optimization problems. But be aware of the licencing. Dieter __ R-help@r-project.org mailing list htt

Re: [R] c() coverts real numbers to integers?

2009-03-07 Thread Dieter Menne
John Poulsen ufl.edu> writes: > I know I am forgetting to do something silly. I typed coordinates in > vectors (as below) but when I call them in R they come out as integers, > and I want them to be real numbers. I have tried using as.numeric, > as.real, etc... but they are still read by R as

Re: [R] NonLinear Programming in R - QUERY

2009-03-07 Thread Dieter Menne
Lars Bishop gmail.com> writes: > I'll appreciate your help on this. Do you know of any package that can be > used to solve optimization problems subject to general *non-linear* equality > constraints. Package DEoptim Dieter __ R-help@r-project.org m

Re: [R] merge data frames with same column names of differe nt lengths and missing values

2009-03-07 Thread Dieter Menne
Steven Lubitz yahoo.com> writes: > > x <- data.frame(item1=c(NA,NA,3,4,5), item2=c(1,NA,NA,4,5), id=1:5) > y <- data.frame(item1=c(NA,2,NA,4,5,6), item2=c(NA,NA,3,4,5,NA), id=1:6) > > merge(x,y,by=c("id","item1","item2"),all.x=T,all.y=T) #my rows are duplicated and the NA values are > reta

Re: [R] Re peated ANOVA or nested ANOVA, or parallel one way ANOVA six times?

2009-03-06 Thread Dieter Menne
janey ding hotmail.com> writes: > Sorry! > > In my experiment, there are two type of soil ( soil F and soil D), each half > of them were subjected to steam sterilize (result in FS and DS soil). A > equal volume of soil from two of the four soil types (F, D, FS, DS) were > mixed as follows: F+F,

Re: [R] RV: help

2009-03-06 Thread Dieter Menne
Jani Lobo grupocomar.com> writes: > I want to estimate the survival mean of a few specific teams. I'm trying to > calculate it through a Kaplan Meier estimator. For doing so, I load the > "survival" package and run the following instructions: > > "options(survfit.print.mean=TRUE)" allows s

Re: [R] sm.options

2009-03-06 Thread Dieter Menne
Viviana Ruiz gmail.com> writes: > I am doing kernel density plots, and am trying to make the lines thicker. I > comparing three groups, in sm.density.compare. I tried changing lwd to make > the line sthicker right on the density compare call, but was not able to do > it. There is not an option

Re: [R] Bug in by() with dates as levels?

2009-03-05 Thread Dieter Menne
hadley wickham gmail.com> writes: > > You might also want to take a look at the plyr package: > > install.packages("plyr") > library(plyr) > ddply(x, .(A, D), function(df) sum(df$Z)) > dlply(x, .(A, D), function(df) sum(df$Z)) > I recently for the first time used the plyr package in an introd

Re: [R] reading scanned graphs

2009-03-04 Thread Dieter Menne
Frank E Harrell Jr vanderbilt.edu> writes: > For Linux I recommend the engauge-digitizer package. It's Java and works under Windows. Dieter __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting g

Re: [R] Diff btw percentile and quantile

2009-03-04 Thread Dieter Menne
megh yahoo.com> writes: > To calculate Percentile for a set of observations Excel has percentile() > function. R function quantile() does the same thing. Is there any > significant difference btw percentile and quantile? If you check the documentation of quantile, you will note that there are

Re: [R] Selecting one row or multiple rows per ID

2009-03-04 Thread Dieter Menne
Vedula, Satyanarayana jhsph.edu> writes: > > > I need to select one row per patient i in clinic j. The data is organized > similar to that shown below. > ... > If patient has outcome recorded at visit 2, then outcome = outcome >columns at visit 2 > If patient does not have visit 2, then ou

Re: [R] RWeka dataset running

2009-03-03 Thread Dieter Menne
ahmet kocyigit hotmail.com> writes: > > I'm a newbie to R and will work on it from now. I'm preparing my thesis based on rpart and RWeka. > I've managed to run already existing test data but couldn't import a data set from out of R. > There is a special Export/Import manual R-data.pdf in Rxxx/

Re: [R] Characters To Matrix

2009-03-03 Thread Dieter Menne
Bob Roberts yahoo.com> writes: >I have a data set that I read in using readLines, but R classes it as a character string. The data set has > multiple columns and rows, but I only want the first 2 columns in the form of a matrix. Thanks for the help. data(CO2) str(CO2) col2 = as.matrix(CO2[,1

Re: [R] [R Question]: regarding exponential axis

2009-03-02 Thread Dieter Menne
Sun, Larry amgen.com> writes: > How to plot exponential axis (i.e. equally spaced .1, 1, 10, 100..)? > I want to use it with xyplot; it would be great if there is a detailed > example. > Thanks a lot! > Sound like what is usually called a logarithmic axis. See parameter log in scales of xyplo

Re: [R] ways to put multiple graphs on single page (using ggplot2)

2009-03-02 Thread Dieter Menne
Juliet Hannah gmail.com> writes: > > Hi, Here are three plots: > (How to layout on one page) library(ggplot2) data(diamonds) randind <- sample(nrow(diamonds),1000,replace=FALSE) dsmall <- diamonds[randind,] a= qplot(carat, data=dsmall, geom="histogram",binwidth=1) b= qplot(carat, data=dsmall,

Re: [R] Bold Face in Plot

2009-03-02 Thread Dieter Menne
Rau, Roland demogr.mpg.de> writes: > I am trying to plot some text in bold face which works fine: > > plot(0:1,0:1,type="n") > text(x=0.5, y=0.5, labels=expression(bold("the-actual-string"))) > > Now when I try to do the following, the displayed text reads thestring: > > thestring <- "the-act

Re: [R] ave and grouping

2009-03-02 Thread Dieter Menne
Patrick Hausmann uni-bremen.de> writes: > # I have a DF like this: > sleep$b <- c(rep(8,10), rep(9,10)) > sleep$me <- with(sleep, ave(extra, group, FUN = mean)) > sleep > > # I would like to create a new variable > # holding the b-th value of group 1 and 2. > > # This is not what I want, it

Re: [R] smoothing a matrix (interpolate in plane)

2009-03-01 Thread Dieter Menne
Žroutík gmail.com> writes: > I'd like to smooth a matrix to dismiss spikes and to interpolate in plane > interp.loess in package tgp works nicely for me. Dieter __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] Adjusting confidence intervals for paired t-tests o f multiple endpoints

2009-02-28 Thread Dieter Menne
Erich Studerus bli.uzh.ch> writes: > > > Dear R-users, > > In a randomized placebo-controlled within-subject design, subjects recieved > a psycho-active drug and placebo. Subjects filled out a questionnaire > containing 15 scales on four different time points after drug > administration. In or

Re: [R] formula formatting/grammar for regression

2009-02-27 Thread Dieter Menne
BKMooney gmail.com> writes: > > > Also, when I tried to use nls, I get an error: > nls(ypts ~ exp(xpts)) > Error in getInitial.default(func, data, mCall = as.list(match.call(func, : > no 'getInitial' method found for "function" objects > > If someone could please point out what I am doin

Re: [R] formula formatting/grammar for regression

2009-02-27 Thread Dieter Menne
Brigid Mooney gmail.com> writes: > I am doing some basic regression analysis, and am getting a bit > confused on how to enter non-polynomial formulas to be used. .. > But am confused on what the formula should be for trying to find a fit > to y = A*exp(r*x). If this example is just a placeholder

Re: [R] Help with Latex (Hmisc)

2009-02-24 Thread Dieter Menne
Frank E Harrell Jr vanderbilt.edu> writes: ... Word and pdf > It depends on how you copy. By all means use Insert ... Picture ... > from file and directly insert pdf. Please, tell me how you got this to work. Dieter __ R-help@r-project.org maili

Re: [R] Inference for R Spam

2009-02-24 Thread Dieter Menne
> > > And, since my son asked me and I am basketball ignorant: Why are > > basketball scores mostly much too close to equality? The arguments > > (loose power when leading) manchester.ac.uk> writes: > Or: Once you are in the lead, become much more defensive against > attacking play by the

Re: [R] Inference for R Spam

2009-02-24 Thread Dieter Menne
Tony Breyal googlemail.com> writes: >Cheers for that information; I've just registered for the useR meeting >in London and then about 10 minutes later got that same bit of spam >too which made me a wee bit suspicious. Welcome in the Fooled by Randomness society. 2:0 is a bit away from significa

Re: [R] Odfweave

2009-02-24 Thread Dieter Menne
Irina Ursachi itwm.fraunhofer.de> writes: > > Dear all, > > When using odfWeave, I get the following error: > > odfWeave(file.in,file.out) > Copying Example2.odt > Setting wd to /tmp/Rtmp8ekeDC/odfWeave24145003519 > Unzipping ODF file using unzip -o Example2.odt > Archive: Example2.od

Re: [R] odfWeave problem "Error in xmlEventParse"

2009-02-24 Thread Dieter Menne
Fredrik Karlsson gmail.com> writes: > > Dear list, > > Sorry for bothering you with a pure odfSweave question, but I just ran > into a problem that I cannot find the cause of. > Anyonse seen this before? This file "used to work", but not anymore. > > Would apreciate all the help I could get. >

[R] Inference for R Spam

2009-02-24 Thread Dieter Menne
Dear List, I registered for the useR conference in Rennes today; half an hour after the confirmation I received a first "requested newsletter" from a company selling a product named "Inference for R". This coincidence might be spurious. Or not, depending on frequency. Dieter ___

Re: [R] Fourier Row and spectral analysis

2009-02-23 Thread Dieter Menne
Peterko gmail.com> writes: > Original serie have 225 observing, but program use only 224. > the most domain frequencies are f1=1/224 f2=1/122 f3=1/74,66 f4=1/56 > f5=1/24,88 > When i know these frequencies a can do, new variable ,nubmers of variavle is > 2*nubmer of frequencies: > t<-1:225 > c1<-

Re: [R] Fourier Row and spectral analysis

2009-02-22 Thread Dieter Menne
Peterko gmail.com> writes: > Hello, i need to find in time serie, k=1,2,3...(how if possible) most domain > frequencies and than smooth original serie by Fourier row y=mean + > a1*Cos(2*pi/freq*t)+b1*Sin(2*pi/freq*t)+a2.. > > numbers a1,b1 ... i will have from simple regresion, i need only

Re: [R] density estimate

2009-02-20 Thread Dieter Menne
Edna Bell gmail.com> writes: > I have a tiny data set: > > xs > [1] 0.7 2.8 0.1 1.9 0.0 1.4 0.2 2.3 0.3 0.2 > > > > and generate the density function for it. > > I would like to replicate this. Is there a straightforward way to do > this, please? See chapter 5 in MASS and the examples in lib

Re: [R] modifying legend in scatterplot matrix R

2009-02-19 Thread Dieter Menne
Julien Beguin ulaval.ca> writes: > scatterplot.matrix(~ a + b + c, >groups=treatment, ... >legend.plot=T, >) > legend("topright",c("fenced","unefenced"), fill=NULL,bty="o",cex=3) > > > In my

Re: [R] Fwd: Percentiles/Quantiles with Weighting

2009-02-17 Thread Dieter Menne
Stavros Macrakis alum.mit.edu> writes: > > Some minor improvements and corrections below > > # Simple weighted quantile > # > # v A vector of sortable observations > # w A numeric vector of positive weights > # p The quantile 0<=p<=1 > # > # Nothing fancy: no interpolation etc.; NA cases not

Re: [R] controlling number of decimals printed in anova tables?

2009-02-16 Thread Dieter Menne
Gabor Grothendieck gmail.com> writes: > > On Mon, Feb 16, 2009 at 11:08 AM, Dieter Menne > > Yes, with as.integer(round(...)) It looks like this: > > > modelFit.glm(berk.mod2) > Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 &

Re: [R] controlling number of decimals printed in anova tables?

2009-02-16 Thread Dieter Menne
Gabor Grothendieck gmail.com> writes: > > Or safer: > > df <- as.integer(round(...)) > Did you try? I believe it is a problem of printCoefmat that has quite a few options for special column, but none for df. Ask Martin Mächler. Dieter __ R-help@r-

Re: [R] extracting standard errors in glm.nb

2009-02-15 Thread Dieter Menne
jpl math.unl.edu> writes: > > How do I extract the standard errors of the coefficients in a glm.nb model? library(MASS) quine.nb1 <- glm.nb(Days ~ Sex/(Age + Eth*Lrn), data = quine) sumnb1 = summary(quine.nb1) sumnb1 str(sumnb1) # tells you about the anatomy. Check for $coefficients Dieter __

Re: [R] help change auto.key colors in xyplot

2009-02-13 Thread Dieter Menne
Dimitri Liakhovitski gmail.com> writes: > the code below works just fine to produce a dotplot. However, I am not > successful changing the color of the lines in the legend (auto.key). > If I add col=..., it only changes the color of the letters in the > legend, but not the color of the lines. I

Re: [R] odfWeave & prettyR

2009-02-13 Thread Dieter Menne
lauramorgana bluewin.ch bluewin.ch> writes: > > Hello, > I've been trying to use odfWeave and prettyR packages to create documents with both text and graphs, but so far I haven't been very lucky... ... > library(odfWeave) I get this message: > > Loading: XML > Error: package 'XML' not loaded

Re: [R] Adding abline in Lattice graph

2009-02-13 Thread Dieter Menne
Sundar Dorai-Raj gmail.com> writes: > > Sorry, that was my lack of understanding on how coplot works. > Try the following: > > coplot(lbxglu~lbxgh|eth, data = reg.dat.5, >panel= function(...) { > panel.smooth(...) > abline(h = 126, col = "red") > abline(v = 6.5, col = "blue"

Re: [R] Problem with loading rJava in R

2009-02-12 Thread Dieter Menne
charis kaskiris.com> writes: > > > Thank you for helping. > > I am running the 32-bit version of R on a 64-bit Windows XP machine. After > reinstalling the rJava package I started getting a windows pop up with the > message: > > Rgui.exe - Unable to Locate Component > "This application has fa

Re: [R] Problem with lmer and wiki example

2009-02-12 Thread Dieter Menne
James Widman mi.nmfs.gov> writes: > > I am trying to duplicate the example by Spencer Graves in the wiki, > using lmer with the Nozzle data. > http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests > However the Chisq value and the fitAB values that are calculated are > different compar

Re: [R] Using nls or nls.lm with a simulation output

2009-02-12 Thread Dieter Menne
Marc Girondot u-psud.fr> writes: > > We would like to fit parameters using a simulation with stochastic > processes as theoretical values. We generate a simple exemple with nls.lm > to see the logic and the problem: > Example removed ... > pStart <- list(a=0.1) > #non-linear fit > library(min

Re: [R] barplot() x axes are not updated after removal of categories from the dataframe

2009-02-12 Thread Dieter Menne
R User R User googlemail.com> writes: > I am plotting categories (X axis) and their mean count (Y axies) with > barplot(). > > The first call to barplot works fine. > > I remove records from the dataframe using final=[!final$varname == "some > value",] > > I echo the dataframe and the records

Re: [R] How to tell if lattice is current device?

2009-02-11 Thread Dieter Menne
Remko Duursma gmail.com> writes: > >> How can I tell if the current device was made with plot() or e.g. > >> levelplot() or another lattice function? > >> dev.cur() does not help me, it just tells me "windows 2". > > Deepayan: > > A device can have both types, although that's rare. > > > > The be

Re: [R] 2 different colors for 2 groups of lines using xyplot

2009-02-11 Thread Dieter Menne
Yu, Xuesong scharp.org> writes: > I am trying to use xyplot to plot several lines in one panel with red > color for a group of lines, say 3 lines and blue color for another > group of lines, say 4. I know this can be easily done using regular > plot function. But i could not figure out how to d

Re: [R] changing settings on a barchart (lattice)

2009-02-11 Thread Dieter Menne
Dimitri Liakhovitski gmail.com> writes: > I apologize - I never used lattice before, so my question is probably > very basic - but I just can't find the answer in the archive nor in > the documentation: > > I have a named numeric vector p of 6 numbers (of the type 6 numbers > with people's names

Re: [R] how to perform 5 level nested anova

2009-02-10 Thread Dieter Menne
Mao Jianfeng gmail.com> writes: > I am new to R. And, I want to perform a multiple nested anova on a large > datasets (with 9448 observations). > > I tried the function ("lmer()" in package "lme4"). But, I failed. Can anyone > help me? > > my dataset("SeedL.txt") was attached. Data are not sort

Re: [R] Truncated labels with mosaicplot

2009-02-10 Thread Dieter Menne
Marcin Kozak gmail.com> writes: > How can I deal with truncated labels in the mosaicplot()? Look at the example: > mosaicplot(~ gear + carb, data = mtcars, color = TRUE) > > Look at "carb" (the number of carburetors): the label "8" is > truncated. How this might be dealt with? As Michael noted

Re: [R] summary of a list

2009-02-10 Thread Dieter Menne
lauramorgana bluewin.ch bluewin.ch> writes: > > Hello, > I'm using the following for loop to find regression curves using a list of functions (formList), a list of .. long non-reproducible code removed > And I got the following error message: > Error in chol2inv(object$m$Rmat()) : l'elemento

Re: [R] Fast fourier transformation

2009-02-10 Thread Dieter Menne
botto uke.uni-hamburg.de> writes: > > here is a practical problem we would like to solve. In a pneumatic post the > acceleration of the capsule is measured and plotted over time. From the > graph achieved we would like to derive some kind of statistic value that > describes the stress the capsul

Re: [R] How to create grouping in the residual plot

2009-02-09 Thread Dieter Menne
Ram Pandit gmail.com> writes: > I am working in a country level data. After running the regression, I would > like to plot the residuals of each observation based on the group created > for a particular variable. > > For example, one of my independent variable is "Income", I would like to > plot

Re: [R] Choosing a random number between x and y

2009-02-09 Thread Dieter Menne
Vie bham.ac.uk> writes: > > Ive been trying to find a function that will allow me to pull out a number > between a minimum and maximum threshold. > > I want a random decimal number between, for example, 0 and 0.5 or 0 and 0.7. > I've been searching everywhere for a function that will allow me t

Re: [R] glmmBUGS: logistic regression on proportional data

2009-02-08 Thread Dieter Menne
David Winsemius comcast.net> writes: > >> > >> Newdat<-data.frame(Newtree=rep(1:3, each=20), Newsect=rep(c("a","b"), > >> each=10), Newdist=rep(1:5, 2), > >> y=rpois(60,2), tot=rep(c(14,12,10,8,6), 12)) > >> > >> yseed<-cbind(Newdat$y, Newdat$tot) > >> > >> mod<-glmmBUGS(yseed~N

Re: [R] Using panel.grid in barchart

2009-02-08 Thread Dieter Menne
Christopher Jones MIT.EDU> writes: > I'm trying to do an example in Deepayan Sarkar's Lattice book. It > involves making a barchart based on the Titanic dataset. I can get > the barchart to plot fine; however, when I try to edit panel.grid, I > get an error: > > > titan<-barchart(Class

Re: [R] glmmBUGS: logistic regression on proportional data

2009-02-08 Thread Dieter Menne
John Poulsen zoo.ufl.edu> writes: > I am trying to run a logistic regression with random effects on > proportional data in glmmBUGS. I am a newcomer to this package, and > wondered if anyone could help me specify the model correctly. > > I am trying to specify the response variable, /yseed/,

Re: [R] updating contents of a package

2009-02-08 Thread Dieter Menne
Markus Loecher gmail.com> writes: > I read through the "Writing R Extensions" document and am able to now create > my own packages/libraries which so far are just well documented collections > of my own R functions. I use package.skeleton() and the tools package to > build these packages. > Howev

Re: [R] lme() direction

2009-02-07 Thread Dieter Menne
Mike Lawrence thatmike.com> writes: > > And if I decided to ignore the "type" variable altogether and simply > use the continuous "valence" variable, this is what I'd use? > > summary(lme( > fixed = rt~valence*color > , data = a > ,random = ~1|id > )) > > I also have continuo

Re: [R] lme() direction

2009-02-07 Thread Dieter Menne
Mike Lawrence thatmike.com> writes: > >Would it improve things if "type" were a continuous variable rather >than categorical? I chose words at the extreme ends of a valence >rating scale but I still have the raw valence ratings for each word. > > > > > With the interaction, the extreme would be

Re: [R] lme() direction

2009-02-07 Thread Dieter Menne
Mike Lawrence thatmike.com> writes: Thanks for the excellent reproducible sample set! > I'm most interested in the interaction between color and type, but I > know that there is likely an effect of word. Yet since word is not > completely crossed with type, simply adding it to an aov() won't wor

Re: [R] Output results to a single postscript document

2009-02-07 Thread Dieter Menne
Pele yahoo.com> writes: > I have been trying to output all my results (text, plots, etc) into the same > postscript file as > > one document, but have been unable to...Can anyone help me improve my code > below so that I can > > accomplish this? Currently I have to output them separately the

Re: [R] Target Plot?

2009-02-06 Thread Dieter Menne
Jason Rupert yahoo.com> writes:   > Also, do you know how to make the radial rings "black" and a little thicker?  A co-worker requested that the be made a little darker and easier to see. >   library(plotrix) ions<-c(3.2, 5.0, 2.0) ion.names<-c("Na","Ca","Mg") radial.plot(ions, labels=ion.n

Re: [R] See source code for survplot function in Design package

2009-02-06 Thread Dieter Menne
Marc Schwartz comcast.net> writes: > getAnywhere("function.method") > Which also works without the "" getAnywhere(xtabs) getAnywhere("xtabs") Dieter __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] eval and as.name

2009-02-05 Thread Dieter Menne
Wacek Kusnierczyk idi.ntnu.no> writes: > equal? you mean equivalent? mostly, yes. briefly, this is why: > > 1. a copy-over from other programming languages; > 2. to avoid learning yet another operator; > 3. after having learned the other operator, to avoid that ugly operator; > 4. after an r

Re: [R] Sweave and backslashes

2009-02-05 Thread Dieter Menne
cls59 mac.com> writes: > I am an avid Sweave user and I am trying to pretty print floating point > numbers for latex output. For example in my document, I would like: > > 4.2\cdot 10^-{8} > > Instead of: > > 4.2e-08 > > The Hmisc package has a nice function for doing this- but Hmisc has a ton

Re: [R] lattice key inside panel

2009-02-05 Thread Dieter Menne
Deepayan Sarkar gmail.com> writes: > > One other thing that has so far eluded me is how to specify two > > different kinds of symbols on a single column. The key above shows > > > > text line text rectangle > > > > but I would really like to get > > > > text line > > text rectangle > > > > but th

Re: [R] See source code for survplot function in Design package

2009-02-05 Thread Dieter Menne
Eleni Rapsomaniki mail.cryst.bbk.ac.uk> writes: > I know one way to see the code for a hidden function, say function_x, > is using default.function_x (e.g. summary.default). But how can I see > the code for imported packages that have no namespace (in this case > Design)? Just type the nam

Re: [R] Sweave and \Sexpr{}

2009-02-05 Thread Dieter Menne
Felipe Carrillo yahoo.com> writes: > It is working now although I still have to specify the path to Sweave.sty \usepackage{C:/R/R-2.8.1/share/texmf/Sweave} > and it doesn't evaluate \Sexpr{} when I simply use \usepackage{Sweave}. I placed Sweave.sty in the same > folder where my rnw and tex file

Re: [R] split dataset randomly in prediction and validation set

2009-02-05 Thread Dieter Menne
Els Verfaillie ugent.be> writes: > For a geostatistical analysis, I would like to split my dataset randomly > into 2 parts: a prediction set (with 2/3 of my data) and a validation set > (with 1/3 of my data). Both datasets will thus contain different data. Any > suggestions? Normally, you will

Re: [R] Sweave and \Sexpr{}

2009-02-05 Thread Dieter Menne
Felipe Carrillo yahoo.com> writes: > I am trying to create a dynamic latex table using \Sexpr{} but it's > not evaluating it. I also tried the > example below without Sweave and also fails. I have also copied the > Sweave.sty to my working directory but > nothing seems to work. Do I need to hav

Re: [R] How to import HTML and SQL files

2009-02-04 Thread Dieter Menne
Arup gmail.com> writes: > I can't import any HTML or SQL files into R..:confused: Also confused. HTML and SQL are like apples and bugs. For HTML (assume you want to extract stock quotes from a site) -- If you have strict XHTML, using package XML might be the best choice, but I doubt you ge

Re: [R] reading .odf spreadsheet into R

2009-02-03 Thread Dieter Menne
stephen sefick gmail.com> writes: > I have searched the archives and I did not find the answer to my > question. Is there a way to read in a .odf spreadsheet without > modification to a .csv file. I am analyzing my classes scores on > their first exam, and would like to read the grade book in w

Re: [R] Lattice histogram with vertical lines

2009-02-03 Thread Dieter Menne
David Scott auckland.ac.nz> writes: > > > I would like to add some vertical lines to a lattice plot of histograms. > What I am after is a lattice version of abline(v = 1234). The lattice > histogram plot is just: > > histogram( ~ LTSE | approach, data = arrivals) > As a little variation on

Re: [R] R DCOMServer Manager 2.0

2009-01-30 Thread Dieter Menne
rbs.com> writes: > We use R DCOMServer Manager 2.0, this application uses RGUI 2.3.0 as > one of its components. We are planning to upgrade to the latest version > of RGUI 2.8.0 but found that R DCOMServer manager is not > compatible with the latest version. Check the messages at: http://www.

Re: [R] putting match.call to good use

2009-01-28 Thread Dieter Menne
Prof Brian Ripley stats.ox.ac.uk> writes: > The way to do this is eval(substitute()). E.g. from the new Rd2HTML > What is Rd2HTML? Dieter __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Please Help for Augmented Prediction Plot

2009-01-28 Thread Dieter Menne
Chun-Hao Tu hotmail.com> writes: > Hi R users,I have a question about augmented prediction plot (?augPred). > The covariate of my data set is c(0, 0.01, 0.1, 1, 10, 100, 1000) and I > have fitted a nonlinear mixed effects model. I use > plot(augPred(out.nlme)) to get the augmented predicti

Re: [R] need help combining two datasets

2009-01-28 Thread Dieter Menne
Adam D. I. Kramer ilovebacon.org> writes: > > You probably want the merge function. > > ?merge > > On Wed, 28 Jan 2009, Somani, Dinesh K wrote: > > > I have two CSV files, one with daily stock returns using method A {date, > > stock, returnA, some uninteresting columns}, and another with meth

Re: [R] HMISC package: wtd.table()

2009-01-26 Thread Dieter Menne
Frank E Harrell Jr vanderbilt.edu> writes: > > (wtd.table(FamTyp.kurz,weigths=HGEW,normwt=FALSE,na.rm=TRUE)) > > That is one solution. The other is to spell 'weights' correctly Have pity with us German speakers. It was such a paing to learn "th" that we cannot resist to apply it whenever po

Re: [R] can't load rJava in R 2.8.1 on Windows XP

2009-01-26 Thread Dieter Menne
Duncan Murdoch stats.uwo.ca> writes: > I don't know what's going wrong on your system. I added a browser() > call to the .onLoad function in R/windows/FirstLib.R on my system, and I > see it successfully gets JAVA_HOME from the registry. It gets a number > of other files, then adds these pat

Re: [R] can't load rJava in R 2.8.1 on Windows XP

2009-01-24 Thread Dieter Menne
Duncan Murdoch stats.uwo.ca> writes: > >> shell("set JAVA_HOME") > >> > > JAVA_HOME=C:\Program Files\Java\jre6 > > > > That doesn't last beyond the shell call, as far as I know. It starts a > process to run the shell, sets the environment variable in that process, > then the process d

Re: [R] plotting curve in xYplot -- using panel.curve

2009-01-23 Thread Dieter Menne
John Poulsen wrote: > > I am trying to plot a curve over points plotted with se's in xYplot (see > example below). I can get Figure 1 below to plot the data with error. > However, I keep getting a the error message > > "Error using packet 1 object "y" not found" > > Can anyone see what I

Re: [R] can't load rJava in R 2.8.1 on Windows XP

2009-01-23 Thread Dieter Menne
Duncan Murdoch stats.uwo.ca> writes: > > Slightly different for me, but failure: > > > > jvm.dll is on path (checked successfully with "where jvm.dll") > > From within R? Remember that the PATH is process-specific. Maybe your > shell added something that put jvm.dll on its path, but you di

Re: [R] can't load rJava in R 2.8.1 on Windows XP

2009-01-23 Thread Dieter Menne
Duncan Murdoch stats.uwo.ca> writes: > > If you look at rJava:::.onLoad now that you can load the package, > you'll see that it goes through a number of tests to try to find the > right path. Slightly different for me, but failure: jvm.dll is on path (checked successfully with "where jvm.d

Re: [R] 3d scatter plot with both error bars and a flexibly fitted surface

2009-01-23 Thread Dieter Menne
Sean Zhang gmail.com> writes: > I, an entry level R user, wonder how make a 3d scatter plot with both error > bars and a flexibly fitted surface. For regular grid data, I found interp.loess in package tpg easiest to use to compute the fine grid required for the plot. For irregular data, check

Re: [R] dimnames in pkg "ipred"

2009-01-22 Thread Dieter Menne
Häring, Tim (LWF lwf.bayern.de> writes: > > OK, the information I send to the list were rather sparse. Sorry for that! > I just tried the command with the recent Version of R and ipred. The error message is the same. > I want to create a classification model. My data consist of 5414 observations

Re: [R] "latex" in Hmisc: cell formating

2009-01-22 Thread Dieter Menne
Tao Shi hotmail.com> writes: > ## I'm using R 2.8.0 on WinXP, Hmisc_3.4-3 > > > table1 <- matrix(10, 180,7) > > cell.format <- matrix("", ncol=7, nrow=180) > > cell.format[c(seq(3,180,6),seq(4,180,6)),] <- "color{red}" > > cell.format[c(seq(5,180,6),seq(6,180,6)),] <- "color{green}" > > > > lat

Re: [R] extension of ftable

2009-01-22 Thread Dieter Menne
Santosh gmail.com> writes: > Is there way to combine grouped and consolidated summaries into a formatted > multidimensional summary table? > e.g. > > malefemale total > country state > usa > new jersey10 20 30 > new york 23

Re: [R] using R how to read a one column alone from a database table from MySQL

2009-01-15 Thread Dieter Menne
sankar82 tkk.fi> writes: > i have a created a database table in MYSQL consisting of 11 columns. > throught RMYSQL i managed to read the entire table in R. but i have few > qureries which i need solutions...here they are: > > 1. Using R how to read a one column alone from a database table from MY

Re: [R] noise in time series

2009-01-14 Thread Dieter Menne
David Riano ucdavis.edu> writes: > I have two time series. Both measure the same thing and I would like > to determine which one is noisier. > Would it be a good measure of the noise in each time series the > absolute lag difference? > Is this a good measure? Any other measure I could use? You

Re: [R] Interface to open source Reporting tools

2009-01-14 Thread Dieter Menne
srinivasa raghavan gmail.com> writes: > I am interested to generate dashboard and reports based on data from MS > Access. These reports need to be posted on a weekly basis to the web. The > reporting interface should provide facilities for "what if" scenarios. > > Is it possible to interface

Re: [R] merge table rows (\multirow)

2009-01-12 Thread Dieter Menne
Felipe Carrillo wrote: > > I am trying to merge the 'Month' column using \multirow. For example for > the column 'Week' I want July to be merged into one row(weeks 27,28,29,30) > and so on for the following weeks. > Below, I am creating a PDF using Sweave, MikTex,R-2.8.1 and windows XP to > sho

Re: [R] merge table rows (\multirow)

2009-01-12 Thread Dieter Menne
Felipe Carrillo yahoo.com> writes: > I need help merging rows. > I am trying to merge the 'Month' column using \multirow. For example for the column 'Week' I want July to be > merged into one row(weeks 27,28,29,30) and so on for the following weeks. > Below, I am creating a PDF using Sweave, MikT

<    2   3   4   5   6   7   8   9   10   11   >