Re: [R] Query about wilcox.test() P-value

. As to your example: there is no sense at all in doing a test on such data (other than to satisfy some hypothetical fanatical journal editor). -Peter Ehlers My version for R is 2.11.1 (2010-05-31) running on x86_64 GNU/Linux (RHEL). Thanks in advance for any help with this. Govind

Re: [R] levene.test

your groups on. SPSS uses the mean, car uses the median. -Peter Ehlers __ R-help@r-project.org 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

Re: [R] Boxplot: Scale outliers

to the boxes scaled with 0.1. Am I missing something? Maybe you're looking for the 'outcex' parameter? See ?bxp. -Peter Ehlers Thanks in advance! Robert __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Calculating weibull's distribution mean, standard deviation, and variance

Try fitdistr() in pkg MASS. -Peter Ehlers On 2010-07-12 11:17, Oscar Rodriguez wrote: Dear R community: Sorry if this question has a simple answer, but I am a new user of R. Do you know a command, or package that can estimate the weibull distribution's mean, standard deviation

Re: [R] KLdiv produces NA. Why?

of the following may be the case: 1. your R and/or flexmix is outdated; 2. dependencies were not installed (modeltools?) 3. there's something weird about your x[,] -Peter Ehlers Best, Ralf __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] strange floor rounding

On 2010-07-09 8:27, Trafim Vanishek wrote: Thanks everybody for referring me to FAQ 7.31 but I don't see how to solve it. I am giving a concrete number and I need to get 58 not 57. Seems, there is no way? Sure there is: round(100*0.58). -Peter Ehlers On Fri, Jul 9, 2010 at 4:05 PM, David

Re: [R] nls error regarding numerics vs logicals

be a bit simpler. -Peter Ehlers __ R-help@r-project.org 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

Re: [R] Not nice behaviour of nlminb (windows 32 bit, version 2.11.1)

Actually, it looks like any value other than 1.0 (and in (lower, upper)) for start will work. -Peter Ehlers On 2010-07-09 14:45, Ravi Varadhan wrote: Setting abs.tol = 0 works! This turns-off the absolute function convergence criterion. nlminb( objective=function(x) x, start=1, lower=-2

Re: [R] Boxplots over a Scatterplot

always most satisfying. And don't let my inbred curmudgeonliness prevent you from asking again when the need arises. -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] Inverting a scale(X)

On 2010-07-07 19:19, Godfrey van der Linden wrote: On 2010-07-03, at 17:43 , Peter Ehlers wrote: On 2010-07-03 0:05, Godfrey van der Linden wrote: G'day, All. I have been trying to trackdown a problem in my R analysis script. I perform a scale() operation on a matrix then do further work

Re: [R] Column header strategy

~ FTIR or y ~ X + Z. Note the word _matrices_; you may not have set up your data correctly. Compare the 'yarn' dataset str(yarn) with your data str(micetitletest) And, as David says, don't use D for the name of your predictor matrix (although it will probably work). -Peter Ehlers

Re: [R] how to define a function in R

Chapter 10 of 'An Introduction to R' has the suggestive title Writing your own functions. Isn't that the *first* place you would look? -Peter Ehlers On 2010-07-06 20:57, jd6688 wrote: Thanks Joshua for the example which has been great help me to start. Wish you the best Jason Ding On Tue

Re: [R] quantiles on rows of a matrix

. (And you would probably want trim=0.025 for your application.) -Peter Ehlers Stumped. Any help appreciated. Thanks. Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

Re: [R] Boxplots over a Scatterplot

Would it really be so hard to provide commented, minimal, self-contained, _reproducible_ code ? -Peter Ehlers On 2010-07-07 15:16, fsch wrote: Hello- I'm new to R, coding and stats. (Oh no.) Anyway, I have about 12000 data points in a data.frame (dealing with dimensions and geological

Re: [R] How to Plot With Different Marker ( ‘x ’ and ‘o’) Based on Condition in R

) -Peter Ehlers __ R-help@r-project.org 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] finding subroutines

elements with getNativeSymbolInfo(setup_starma) getNativeSymbolInfo(free_starma) or with stats:::R_setup_starma -Peter Ehlers Thanks for any help. Have a great weekend. Sincerely, Erin __ R-help@r-project.org mailing list https

Re: [R] Conditional Splitting a Vectors into Two Vectors

extraction: x[count 0] x[!(count 0)] ##or equivalent -Peter Ehlers __ R-help@r-project.org 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

Re: [R] Assign Formulas to Arrays or Matrices?

thing as a 3-element vector. Why not use a list? ps. for (i in seq(from=1, to=31, by=1)) is equivalent to for(i in seq_len(31)) -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Could not find createData function

() function and read.madata used to have arguments 'probeID' and 'pmt', but those are now 'probeid' and, I think, 'intensity'. Try the demo in your (*unstated*) version of maanova. -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Error message using mi() in mi package

see why I can sometimes impute a variable labeled as ordered-categorical and sometimes cannot. Thanks! Andrew Miles This looks suspiciously like a syntax problem. I would get my text editor to search for 'c14ordered' in the code. You might have missed some punctuation. -Peter Ehlers

Re: [R] plotting data with ellipse confidence intervals

, then you can't plot CIs and ellipses. If you have the original paired data, then have a look at the car::ellipse function and friends. -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Issue with write.table and read.table : I'm not getting out what I put in

== in situations such as this. Try instead: all.equal(test2, allpos, check.attributes = FALSE) (why the check.attributes argument may be needed is left as an exercise) -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] help with predict.lda

problems arise, *first* use str() on your objects to see what you've got. 5. finally, do provide *reproducible* code; here, you could have used the example on the help page. -Peter Ehlers Thanks so much! __ R-help@r-project.org mailing list

Re: [R] Separating out data values

quadrants. Getting R to count the points is left as an exercise (which should be easy after working through An Introduction to R). -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Inverting a scale(X)

by your further work, then I'm afraid you're out of luck. -Peter Ehlers __ R-help@r-project.org 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

Re: [R] How best to set library search path so user libraries come first

as the second element of .libPaths(). (This is on Windows.) -Peter Ehlers __ R-help@r-project.org 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

Re: [R] rms::ols I(.) in formulas

/2010-June/241587.html -Peter Ehlers On 2010-07-02 4:51, Otto Kässi wrote: Hi, Lexi! I am aware that lm() is the standard way to do ols regression in R. The reason why I opted for rms::ols() is that later on in my work I need some rms functions which are not available for lm(). In retrospect

Re: [R] rms::ols I(.) in formulas

On 2010-07-02 7:36, Frank E Harrell Jr wrote: ols has always been fussy about this. And don't use poly; use pol with rms/Design. Frank Quite right; I did mean pol(), not poly(). Darn fingers sometimes just won't obey. -Peter Ehlers On 07/02/2010 07:28 AM, Peter Ehlers wrote: Otto

Re: [R] xyplot: key inside the plot region / lme: confidence bands for predicted

=0.8, corner=c(0,1))) See description of 'key' in ?xyplot. -Peter Ehlers __ R-help@r-project.org 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

Re: [R] merging plot labels in a lattice plot

as to how I could acheive this would be appreciated Thanks, I like the useOuterStrips() function in the latticeExtra package for this: useOuterStrips(histogram(~d | f1 + f2, x)) -Peter Ehlers __ R-help@r-project.org mailing list https

Re: [R] s.e. of arima

is the estimated variance matrix of the coefficients. If you're talking about some other arima() function, then I have no idea. -Peter Ehlers Thank you Sabrina __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] plot focus

You could just open two devices. See ?dev.cur for an example. -Peter Ehlers On 2010-06-30 4:26, pdb wrote: I am doing calculations in a loop and then plotting the results by adding a point to each of 2 charts at the end of the loop. Its very informative as you can see the progression

Re: [R] plot focus

of screen(), why not use lattice with trellis.focus(). Have a look at the examples in ?trellis.focus. Also, there are some examples in ?xyplot.zoo in the zoo package. -Peter Ehlers Contact Details:--- Contact me: tal.gal

Re: [R] What is wrong with this code?

want to get a FALSE value. Is there a way to do this that's less messy than what I came up with? Thanks again. Str1 - abc Str2 - abcd Str3 - efgh sapply(list(Str2, Str3), grepl, pattern = Str1) [1] TRUE FALSE I think I'd just go with grepl(Str1, c(Str2, Str3)) -Peter Ehlers

Re: [R] string question

cat() is probably what you want, but note that print() has a 'quote=' argument that you could set to FALSE: print(qr2, quote = FALSE) See ?print.default -Peter Ehlers On 2010-06-30 13:16, Phil Spector wrote: Paul - When you print a string, it escapes the quotes with a backslash. That's

Re: [R] Embed function strips out date index

series.d (a matrix) to have the same rownames: rownames(series.d) - rownames(series) But I'm not sure how much sense that makes, especially for columns 3 and 4. -Peter Ehlers Thank you, Manussawee __ R-help@r-project.org mailing list https

Re: [R] What is wrong with this code?

certainly should upgrade. But in the meantime, try wrapping the grep call in length(): if(length(grep(Str1, Str2)) 0) { stuff } else { other stuff } -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] mhplot error with test example: ylim not found

the function (it's short and the fix is easy). -Peter Ehlers On 2010-06-23 8:47, vaneet wrote: Hello all, I am trying to make a genome association plot for p-values related to SNPs and was fortunate to find that R contains a package that produces Manhattan plots which is what's preferred for my

Re: [R] become a member of R user community

the homepage 'Foundation' link or go directly to http://www.r-project.org/foundation/membership.html Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] Remove squares from scatter3D

not tried to use Rcmdr's menus. -Peter Ehlers Best, //M __ R-help@r-project.org 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

Re: [R] About normality tests...

you have good reason to suspect that the cdf lies above or below the Normal cdf. -Peter Ehlers Thank you, Ralf __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] installing packages from a directory on the computer

) : cannot open the connection Please help Just use file.choose() for the pkgs argument (and perhaps set repos=NULL). -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] applying ifelse to dataframe

(11,12,13,16,19,50,27,24,22,26) mask-apply(DF[,3:5],2, %in% ,targets) is.na(DF[3:5]) - !mask -Peter Ehlers mask-data.frame(a=TRUE,b=TRUE,!mask) DF[mask]-NA Error in FUN(X[[1L]], ...) : only defined on a data frame with all numeric variables DF2-data.frame(DF[,3:5]) mask-apply(sample,2

Re: [R] installing packages from a directory on the computer

have to copy the file PACKAGES into your download folder. 3. Specify lib= if you don't want to install into the default directory. 4. You may get warnings about unavailable packages. -Peter Ehlers On 2010-06-22 7:38, David Winsemius wrote: On Jun 22, 2010, at 5:48 AM, Subramanian S wrote

Re: [R] Is there a similar R function like stpower in STATA

On 2010-06-20 23:33, zhu yao wrote: Dear R users: In stat, there is a stpower function for power analysis and sample-size determination in survival models. Is there a counterpart in R? library(sos) findFn(power) -Peter Ehlers Thanks Yao Zhu Department of Urology Fudan University

Re: [R] list() assigning the same value to two items

to define an additional variable? Something like z- list(a=1, b=a) or even z- list(a=b=1)? Both commands don't work of course, but I'm just curious whether this is possible in principle (a rather academic question I've to admit) f - function(x){ list(a=x, b=x) } z - f(1) etc. -Peter

Re: [R] mean duration (persistence) of daily data

I'm not sure what exactly you want, but maybe ?rle can help you. -Peter Ehlers On 2010-06-20 15:52, Schmidt Martin wrote: Dear R users A quite simple question: how can I get the mean duration (persistence) of daily data, when the data set looks as follows: [1] 7 7 7 9 9 5 7 5 5 5 9 5 5 6

Re: [R] How to read data from this csv?

. Try count.fields() on your file and see if all lines do have the same number of fields. -Peter Ehlers Thanks __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] Find the 50 highest values in a matrix

m - matrix(round(rnorm(4000 * 2000), 4), nr = 4000) is.na(m) - sample(8e6, 1e6) system.time( idx - which( matrix(m %in% head(sort(m, TRUE), 50), nr = nrow(m)), arr.ind = TRUE)) # user system elapsed # 3.120.193.18 -Peter Ehlers On 2010-06-18 5:13, Dennis

Re: [R] NPMC

are in a table? Is this a data.frame (as required by npmc)? I think that npmc can only handle one numerical and one categorical variable and they must be named var and class. -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] ERROR need finite 'ylim' values

, data = OrchardSprays, subset = treatment != A) I'll be submitting a bug report (and I think the fix is easy). -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] generating ordered, random decimal fractions

I must be missing something; can't you just use, e.g., one of: sort(round(runif(10), 3)) sort(sample(1000, 10, TRUE)/1000) ? -Peter Ehlers On 2010-06-12 7:49, kurt_h...@nps.gov wrote: Greetings How do I do this in R? Checking the Cran site produces a bewildering array of packages

Re: [R] Calculation of r squared from a linear regression

summary.lm reports is clearly described on the help page. See r.squared in the Value section. -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] Misplacement of Greek letter

lines. -Peter Ehlers On 2010-06-10 23:36, beloitstudent wrote: Hello. I am trying to get my axis label to read as follows (The symbol) Delta AUC blah blah... then below it...(some other text) The problem is the Delta symbol shows up beside the (some other text) rather than the AUC. Does any

Re: [R] ff package when reading .csv files

for read.table.ffdf() with FUN set to read.csv. You should not use any FUN= argument in read.csv.ffdf. Either call read.table.ffdf with FUN=read.csv or call read.csv with no FUN specification. Nor is there any need to specify all the other arguments to take their default values. -Peter Ehlers

Re: [R] Cforest and Random Forest memory use

You say you're on a 64-bit box, but are you running 64-bit R? -Peter Ehlers On 2010-06-10 4:36, Matthew OKane wrote: Hi all, I'm having great trouble working with the Cforest (from the party package) and Random forest functions. Large data set seem to create very large model objects which

Re: [R] r code to broaden the boarder of the bars of a histogram

Use truehist() in pkg:MASS. -Peter Ehlers On 2010-06-11 7:09, Andreas Baranowski wrote: To whom it may concern, I have a problem concerning the design of a histogram. How do I change the border widths of the bars of a histogram. The initial command is: hist (punkte,breaks=30, xlab=Punkte

Re: [R] lattice: how to remove ticks from splom()?

It gets my vote. -Peter Ehlers On 2010-06-10 9:30, Bert Gunter wrote: Is the quote below from Deepayan a possible Fortune? It put a sardonic smile on my lips. Bert Gunter Genentech Nonclinical Biostatistics ... You should really get into the habit of reading documentation; it can often

Re: [R] nls model fitting errors

((ldose - B)/C)) / (1 + exp((ldose - B)/C)) with starting values: A = 0.2, B = 3, C = 1 -Peter Ehlers Gregory A. Graves, Lead Scientist Everglades REstoration COoordination and VERification (RECOVER) Restoration Sciences Department South Florida Water Management District Phones: DESK: 561

Re: [R] bug? in stats::cor for use=complete.obs with NAs

before the casewise deletion (else why the na.last=keep?). If anyone knows the rationale and/or can give a reference, I'd be glad to receive such. -Peter Ehlers On 2010-06-09 11:36, hugh.ge...@thomsonreuters.com wrote: Ar, I think I've found a bug in the behavior of the stats::cor

Re: [R] How to add a new plot in the same graph using add=T at the command plot?

that Petr knows better.) -Peter Ehlers On 2010-06-09 9:08, Larissa Lucena wrote: Thanks so much!!! I'm using R for the first time, and so, I have many stupid doubts! Sorry and thanks again! Regards! 2010/6/9 Petr PIKALpetr.pi...@precheza.cz Hi where did you find parameter add=T. You can use

Re: [R] Rglpk

modified function is myfun(). After you source myfun, set its environment with environment(myfun) - environment(Rglpk_solve_LP) -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Draw text with a box surround in plot.

Your transformation assumes that the x- and y-axes are on the same scale. Add 'asp = 1' to your plot() call to set the appropriate aspect ratio. -Peter Ehlers On 2010-06-09 10:13, g...@ucalgary.ca wrote: Rectangle R centered at (x,y) with width 2w and height 2h is given by x1=x-w y1=y-h x2

Re: [R] Draw text with a box surround in plot.

: plot(0:30, 0:30, type='n', asp=1, axes=FALSE, xlab=, ylab=) abline(v=c(0,30), h=c(0,30)) abline(a=0, b=1, col=2) abline(a=30, b=-1, col=4) polygon(c(x1,x2,x3,x4), c(y1,y2,y3,y4), border=blue, lwd=2) polygon(x.rotated, y.rotated, border=red, lwd=2) Works for me. -Peter Ehlers plot(1

Re: [R] Draw text with a box surround in plot.

Maybe this can be display-dependent. I just measured the two diagonals on the plot and they're equal to within less than one mm (possibly closer but that was the accuracy of my measuring device (pencil and strip of paper)). -Peter Ehlers On 2010-06-09 16:57, g...@ucalgary.ca wrote: For me

Re: [R] import text file into R

also that read.delim is convenient for tab-delimited files. -Peter Ehlers ca_los.txt is my tab delimited large text file which contains about 16lakhs observations. thank you On Fri, Jun 4, 2010 at 9:44 AM, Erik Iversoner...@ccbr.umn.edu wrote: Dhanasekaran wrote: Sorry guys

Re: [R] Boxplot: what is shown by default?

. Note that the box limits are the 'hinges' given by fivenum(x) and their difference need not equal IQR(x) which uses quantile(..., type=7). [For your data above, they do happen to coincide.] -Peter Ehlers /thanks in advance for any help! cheers, tom human genetics, bern

Re: [R] lattice, xyplot, using panel.segments by just addressing one panel

reading the posting guide. It may be that trellis.focus() solves your problem. -Peter Ehlers On 2010-06-02 8:46, Doris wrote: Hi R experts, I'm using the xyplot function in lattice to draw a multipanel plot consiting of 5x6 scatterplots. Now I need to link single points in each of those

Re: [R] pdf function, resize xyplot plot automatically

, 45, 60, 90, 150, 210, 270), cex = 0.7))) dev.off() -Peter Ehlers On 2010-06-02 12:34, Afshartous, David wrote: All, When saving plots to a pdf file via the pdf function, I would like to be able to automatically expand the graphics device to achieve the same result as when one

Re: [R] Draw text with a box surround in plot.

On 2010-06-02 16:57, g...@ucalgary.ca wrote: text() can draw text on a plot. Do we have a way/function to draw text with a box surround it? Thanks, There may be a function in some package, but I would just use rect(). ?rect -Peter Ehlers __ R

Re: [R] Faster matrix operation?

all.equal(vx, wx) [1] TRUE Or, for a very slight further reduction in time in the case of larger matrices/vectors: as.vector(tcrossprod(V, xyzs)) I mention this merely to remind new users of the excellent speed of [t]crossprod(). -Peter Ehlers ? -Original Message- From: r-help-boun

Re: [R] data frame manipulation with zero rows

(x) 1) out - x[, c(1,3,2)] # or whatever else out - ddply(x, c(DESCRIPTION,SETTLEMENT), summarise, POSITION=sum(QUANTITY))[,c(1,3,2)] out } f(futures) -Peter Ehlers -Original Message- From: Prof Brian Ripley [mailto:rip...@stats.ox.ac.uk] Sent: Tuesday, June

Re: [R] New behavior of 'write.csv' append argument?

) # It's always a good idea to check the NEWS file. Here is a quote: write.csv[2] no longer allow 'append' to be changed: as ever, direct calls to write.table() give more flexibility as well as more room for error. So use write.table. -Peter Ehlers

Re: [R] Help barplots

() to place horizontal marks. -Peter Ehlers Thank you Jeet __ R-help@r-project.org 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

Re: [R] Issue with assigning text to matrix

Jessica, Two further comments: 1. When you read your data, your character vectors are stored as factors; is that what you want? 2. You may be better off using a list. -Peter Ehlers On 2010-06-01 7:47, Sarah Goslee wrote: On Tue, Jun 1, 2010 at 5:04 AM, Jessica Queree j.j.que

Re: [R] How to make R automatic?

[, mtime]) For (4) see help(Sys.sleep) -Peter Ehlers Any suggestion/better solution? Thank you very much!! Best regards, Ted __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] geepack installation problem?

Forum articles. Has the package been attached to the search tree? See Chapter 13 of 'An Introduction to R. ?library library(geepack) -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] read in data file into R

to first make two files, for each object one? I would go the two-files route. -Peter Ehlers Below I have copied a subset of my data files for ilustration. This seems a very silly question but I just didn't manage to to it. Thanks a lot for the help cheers beni This is a subset of my data file

Re: [R] writing function : can't find an object

looking for 'PosB'. -Peter Ehlers __ R-help@r-project.org 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

Re: [R] Custom axis function in lattice:::xyplot

for you.) -Peter Ehlers __ R-help@r-project.org 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] website address for the pseuso-XLS files

Maura, These are html files. Rename the downloaded file(s) to *.html and open with your favourite browser. Follow links from there. -Peter Ehlers On 2010-05-25 6:24, mau...@alice.it wrote: http://gigamail.rossoalice.alice.it/messages/readMessageFrameset.aspx?DeliveryID=ba40cf18-29db-4404

Re: [R] Cohen's Kappa for beginners

cohen.kappa(cbind(x,y)) works for me. -Peter Ehlers On 2010-05-25 7:42, Scot W. McNary wrote: Hi, It doesn't seem happy with non-numeric data in the data frame version. Maybe a recode would work? x1 - c(1, 2, 3, 1) y1 - c(1, 3, 3, 1) cohen.kappa(data.frame(x = x1, y = y1)) Call

Re: [R] Lattice: relation = 'free' in scales

= ylim, at = list(NA,NULL,NULL,NA,NULL,NULL,NA,NULL,NULL))), # If you want to move panels closer together, you can add between = list(x = -1.5), # to your pplot code. -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] Cohen's Kappa for beginners

Jason, It seems that your version of psych may be outdated. CRAN currently has version 1.0-88, in which wkappa is listed as deprecated. Might be time to upgrade. -Peter Ehlers On 2010-05-25 10:55, Jason Priem wrote: Thanks for you quick responses, all! Bill, it looks like that's just what

Re: [R] doubt about samr siggenes.table$genes.up

of the data you're saving (probably rownames from a larger dataframe or matrix). -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] retrieve path analysis coefficients (package agricolae)

)) invisible(list(Coeff = Coeff)) } -Peter Ehlers Thanks in advance, Zack [[alternative HTML version deleted]] __ __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] import data from a csv file

getting a pairs plot?) 3. Have you read 'An Introduction to R', especially Chapter 12? -Peter Ehlers Hannah __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] lattice 'scales' option help

where I'm going wrong? OR someway to achieve this? bwplot(treatment ~ decrease | factor(rowpos), data=dta, pch=|, scales=list(y=list(alternating=1)) ) You've switched x and y; try x=list(alternating). -Peter Ehlers Any help would be highly appreciated, thank you in advance... Regards

Re: [R] adding one line to a plot

plot(x,dt(x,df=18)) will provide a clue. (What's x?) You probably want x - tstat[order(tstat)] -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] need help in understanding R code, and maybe some math

, Q) -Peter Ehlers Thanks, Barisdad. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] Lattice Wireframe Plot into LaTex

Could be FAQ#7.22. -Peter Ehlers On 2010-05-22 13:00, Lars Bishop wrote: Dear R/LaTex user, I'm simply trying to include a plot created with the Lattice wireframe function into LaTex. I have no problems including other R plots into LaTex by exporting as a Postcript and then including

Re: [R] Levelplot

for a lattice function, look at ?xyplot where the common arguments are described. Here you will find descriptions of 'xlab', 'ylab'. -Peter Ehlers Regards, Benedikt __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Levelplot

on the help page which has no x- or y-labels. It would seem that the OP didn't bother to read the help page. Tsk, tsk. -Peter Ehlers Cheers, Andy __ Andy Jaworski 518-1-01 Process Laboratory 3M Corporate Research Laboratory - E-mail: apjawor...@mmm.com Tel: (651

Re: [R] lattice: How to specify strip text in an xyplot?

(var.name,...) strip.default(var.name=c(plot 1,plot 2),...)) xyplot(y[,1]+y[,2]~x,outer=TRUE, strip = strip.custom(factor.levels = c(plot 1,plot 2))) -Peter Ehlers __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] save in for loop

(filename, temp) save(list=c(filename), file=paste(filename, .rda, sep=)) } -Peter Ehlers - Original Message From: Ivan Calandraivan.calan...@uni-hamburg.de To: r-help@r-project.org Sent: Wed, May 19, 2010 7:56:44 AM Subject: [R] save in for loop Dear users, My problem concerns save

Re: [R] colored venn diagram

Tao, You might also have a look at the venneuler package. -Peter Ehlers On 2010-05-19 21:37, David Winsemius wrote: On May 19, 2010, at 7:15 PM, Shi, Tao wrote: Hi list, This is probably too much to ask, but I'm wondering if there is a ready-to-use function somewhere that allows me

Re: [R] (no subject)

Arantzazu, Your problem is that the data were probably imported from Excel where you had at least one cell containing #N/A. You need to replace those cases in your dataframe with NA. Then you should be able to do as.numeric(as.character(arcilla)). -Peter Ehlers On 2010-05-18 10:07, Arantzazu

Re: [R] avoiding reinstall already installed *package*

* reside. And it's a good idea to have Martin remind us now and again that books themselves are not libraries. But I must confess that I'm no longer much bothered by the misuse. If it ever leads someone astray in their code, then, well, they have only themselves to blame. Cheers, Peter Ehlers

Re: [R] Function that is giving me a headache- any help appreciated (automatic read )

- with(precip.1, tapply(precipitation, gauge_name, cumsum)) b - unlist(b) You'll also find that you need to wrap your qplot() calls in print(). And your if/else logic (at the end of USGS) looks faulty. -Peter Ehlers have a look at ?is.na --- On Tue, 5/18/10, stephen sefickssef...@gmail.com wrote

Re: [R] Adding a row at top of dataframe

of your variables are not numeric or character, you'll get a warning. -Peter Ehlers __ R-help@r-project.org 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

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