[R] print elements' placings from vector

Dear R, Here is a code: z=c(1,3,0.5,6,8,10,2,2,3,4,7,3) z[z>2] I dont want print the elements, but theirs placings in vector. -- View this message in context: http://r.789695.n4.nabble.com/print-elements-placings-from-vector-tp3496645p3496645.html Sent from the R help mailing list archive at

Re: [R] print elements' placings from vector

Yes, thank you. -- View this message in context: http://r.789695.n4.nabble.com/print-elements-placings-from-vector-tp3496645p3497080.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] plot generates graph with coordinantes written over it

Hello R users, Im having this strange problem. http://r.789695.n4.nabble.com/file/n3350024/bad2.png I have txt file of this format (X.X, Y.Y - random numbers, ...) with data points: nm A nm A nm A X.X Y.Y X.X XX X.X X.X Y.Y Y.Y Y.Y Y.Y My script looks like this: require(graphics) require(stats)

Re: [R] plot generates graph with coordinantes written over it

Solution: main=HMF -> main="HMF" -- View this message in context: http://r.789695.n4.nabble.com/plot-generates-graph-with-coordinantes-written-over-it-tp3350024p3350163.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.o

[R] how to label lines

Hello R I can not figure out, how to add text string or number to lines like on this image: http://r.789695.n4.nabble.com/file/n3350452/graph.png I have some data series in text file. My script: require(graphics) require(stats) graf=read.table("file",header =FALSE,sep = "", dec = ",",fileEncodin

Re: [R] how to label lines

Thank you for reply. I think, in script you posted, I have to manually add coordinates for label. I don't want do it manually. I would like to add labels at the points, where the lines are maximally separated. (points of maximum separation). Graphic representation of my data looks very similar to i

Re: [R] how to label lines

Problem solved: require(graphics) require(stats) require(hmisc) HBA=read.table("file",header =FALSE,sep = "", dec = ",",fileEncoding = "", encoding = "unknown", skip=60,nrows=210) attach(HBA) labcurve( list( One= list( V1,V2), Two= list( V3,V4), Three=list( V5,V6), Four= list( V7,V8), Five= l

Re: [R] how to label lines

Nice, thank you. -- View this message in context: http://r.789695.n4.nabble.com/how-to-label-lines-tp3350452p3351605.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinf

[R] pdf device - different page orientation and different plot sizes

Hello R, I would like make single pdf file in which the the pages are are of mixed format. Manual states this code (down). The problem is the second "{}" overwrites the first "{}", but I want to append them. require(graphics) require(stats) HBA=read.table("file.txt",header =FALSE,sep = "", dec =

[R] Finding coordinates for maximum of a function

Hello R, I have data in txt file of n columns like this: x y1 .yn 1 100 2 50 3 10 4 200 5 20 I need to find coordinates for each maxima e.g.: [x,y1],[x,y2]...[x,yn] a=read.table() attach(a) b=y1<-x max(b) returns only value of y but no corresponding x value -- View this message in

Re: [R] Finding coordinates for maximum of a function

It occured to me: x[max(y)] simple isn't it? -- View this message in context: http://r.789695.n4.nabble.com/Finding-coordinates-for-maximum-of-a-function-tp3355369p3356117.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-projec

Re: [R] Finding coordinates for maximum of a function

I must have mixed it up. Thank you. -- View this message in context: http://r.789695.n4.nabble.com/Finding-coordinates-for-maximum-of-a-function-tp3355369p3357460.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mai

[R] adding linear regression data to plot

Hello R, I would like to print regression data in graph. I mean the output from: k=lm(formula,data) summary(k) Or somehow extract and print only coefficients and R-squared. -- View this message in context: http://r.789695.n4.nabble.com/adding-linear-regression-data-to-plot-tp3357946p3357946.h

Re: [R] adding linear regression data to plot

I know I can add line to graph with abline(), but I would like to print R-squared, F-test value, Residuals and other statistics from lm() to a graph. I don't know how to access the values from summary(), so that I can use them in a following code or print them in a graph. -- View this message in

Re: [R] adding linear regression data to plot

Thank you that is very helpful. -- View this message in context: http://r.789695.n4.nabble.com/adding-linear-regression-data-to-plot-tp3357946p3382615.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list ht

[R] Strange R squared, possible error

k=lm(y~x) summary(k) returns R^2=0.9994 lm(y~x) is supposed to find coef. a anb b in y=a*x+b l=lm(y~x+0) summary(l) returns R^2=0.9998 lm(y~x+0) is supposed to find coef. a in y=a*x+b while setting b=0 The question is why do I get better R^2, when it should be otherwise? Im sorry to use the wor

Re: [R] Strange R squared, possible error

It states summary.lm: r.squared R^2, the ‘fraction of variance explained by the model’, R^2 = 1 - Sum(R[i]^2) / Sum((y[i]- y*)^2), where y* is the mean of y[i] if there is an intercept and zero otherwise. Why to use different formula when intercept is set to zero? I tried to compute R^2

Re: [R] Strange R squared, possible error

1) In my very humble opinion R^2 can't be negative, at least for data for which it sound to use linear model. Or the data would have to be utterly wrong to fit them with linear model. 2) I don't want to fit data with linear model of zero intercept. 3) I dont know if I understand correctly. Im 100%

Re: [R] Strange R squared, possible error

Exuse me, I don't claim R^2 can't be negative. What I say if I get R^2 negative then the data are useless. I know, that what Thomas said is true in general case. But in my special case of data, using nonzero intercept is nonsense, and to get R^2 less than 0.985 is considered poor job (standard R^2>

Re: [R] Strange R squared, possible error

Thank you for your very comprehensible answer. I a priori know that model y=a*x+0 is right and that I can't get x=constant nor y=constant. I'm comparing performance of data gathering in my data set to another data sets in which performance gathering is characterized by R-squared . The data in dat

Re: [R] Strange R squared, possible error

Yes they are. I had edited the reply, but It didn't help. Correction: 2)I meant zero slope, no zero intercept. -- View this message in context: http://r.789695.n4.nabble.com/Strange-R-squared-possible-error-tp3382818p3384648.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Strange R squared, possible error

Thats exactly what I would like to do. Any idea on good text? I've consulted severel texts, but no one defined R^2 as R^2 = 1 - Sum(R[i]^2) / Sum((y[i])^2-y*)) still less why to use different formulas for similar model or why should be R^2 closer to 1 when y=a*x+0 than in general model y=a*x+b. fr

Re: [R] Strange R squared, possible error

Thank you for our reply. It's a pity, that 2 variables defined by different formula have same name. If the variables had been named differently, I wouldn't have problem at all and it looks like it's done on purpose. Because I test a quality of data (performance of collecting data) not a model which

[R] ploting continuos function

Hello R, I would like to plot continuos function (e.g. y=x^2) without creating a sequence for x. Is there a possibility to somehow specify the range - x belongs to <-100,100>? -- View this message in context: http://r.789695.n4.nabble.com/ploting-continuos-function-tp3467533p3467533.html Sent fr

Re: [R] Odp: ploting continuos function

Nice. Thank you -- View this message in context: http://r.789695.n4.nabble.com/ploting-continuos-function-tp3467533p3467816.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

[R] filling array with functions

Hello R, I would like to find out how to generate array full of functions. I tried it like this: fv=array(,dim=c(1,10)) V=c(1:10) for (i in 1:10){ fv[i]<-function(x)(V[i]-b*a*x) # b, x are constants. } But it returns: "incompatible types (from closure to logical) in subassignment type fix" --

[R] random typing over text

Dear R, When I'm writing script I probably hit accidently some hot key and then the text in script is being overwritten by subsequent typing. Please what is the hot key or how can I disable it? -- View this message in context: http://r.789695.n4.nabble.com/random-typing-over-text-tp3470428p347

Re: [R] random typing over text

Thank you very much. It was the Insert key. It was very annoying. Actually is this owerwrite function of any use? -- View this message in context: http://r.789695.n4.nabble.com/random-typing-over-text-tp3470428p3470749.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] filling array with functions

Thats not exactly what I hoped for. But for now it has to suffice. More transparent syntax would be nicer. Exactly what I would like to do is: for (i in 1:9){ f[i]<-function(x){ a*x+b) } curve(f[i], 0, 8) sol[i]<-uniroot(f[i],c(0, 8))$root points(sol[i],0,pch=16,cex=1,col="red") } Perhaps is

Re: [R] filling array with functions

Richard, that way I will have to write functions manually and that is not possible for large number of functions. derek -- View this message in context: http://r.789695.n4.nabble.com/filling-array-with-functions-tp3468313p3472885.html Sent from the R help mailing list archive at Nabble.com

Re: [R] filling array with functions

Still I haven't had any luck yet. How about defining new function and its domain, is it somehow possible? Like this: a<-function(x) x beolongs to natural numbers <0,100> -- View this message in context: http://r.789695.n4.nabble.com/filling-array-with-functions-tp3468313p3473300.html Sent from th

[R] How to extract smoothed plot data?

Dear R, I have this script which from my data points creates smoothed plot: require(graphics) require(stats) A=read.table("some.txt",header =FALSE,sep = "", dec = ",",fileEncoding = "", encoding = "unknown", skip=18,nrows=400) attach(A) plot(A$V1, A$V2,col="white") lines(smooth.spline(A$V1, A$V2,

[R] limit windows() rescale/resize options

I am developing an application that opens an image in a new window using, at times, windows(). I don't want the user to be able to resize the window (and distort the image). The new window contains a menu item called "resize" that contains three options - "R mode", "Fit to window", and "Fixed si

[R] Evaluating an expression

'm guessing that only the last evaluation of the expression is being evaluated, which I can understand. I feel certain that what I want to accomplish can be done programmatically, but am at a loss as to just how to do that. Chances are this has been covered before. If so, apologies. If not

Re: [R] Help with shading a polygon below a segment of a curve (normaldistribution)

Bertolt, The points you send to polygon() do not fully enclose the area you desire. Try adding one more point as such xt <- c(x[(length(x)-cutpoint):length(x)],linepos) yt <- c(y[(length(y)-cutpoint):length(y)],0) polygon(xt, yt, density = 10 ) -Original Message- From: [EMAIL PROTECTE

[R] contourLines() documentation

ction? Thanks, Derek Lacoursiere -- View this message in context: http://www.nabble.com/contourLines%28%29-documentation-tp25151467p25151467.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.et

[R] Converting from "date strings" to internal POSIX format?

,1,1,0,0,0, tz ="UTC")+t32/1000 return(format(timez,"%Y-%m-%d %H:%M:%OS3",tz="CET")) } > render.t32(1245137297746) [1] "2009-06-16 09:28:17.746" Have not succeeded in my attempts to understand POSIX and I thank you for your attention. ~Der

[R] How to convert numbers to words?

I am trying to convert numbers to English words, e.g. 123 -> One hundred twenty three. After some searching, I have been unable to find any info on this. I could certainly write code to do it, but my question is this: Does anyone know of a function that does this? Thanks, Derek McCrae Nor

[R] Iteratively Reweighted Least Squares of nonlinear regression

Dear all, When doing nonlinear regression, we normally use nls if e are iid normal. i learned that if the form of the variance of e is not completely known, we can use the IRWLS (Iteratively Reweighted Least Squares ) algorithm: for example, var e*i =*g0+g1*x*1 1. Start with *w**i = *1 2.

[R] maximum/minimum value of a function

Dear All, Does anyone know how to get the maximum/minimum value of a function using R, I think there should be an instruction, but i can not find it in help files Thank you in advance for your help :) [[alternative HTML version deleted]] __ R

[R] re ference

Hello R users, I have developed code in R that generates populations of non-overlapping random polygons ("polygonal maps") for testing certain stereological estimation methods. I have tried to look for articles or papers online on the generation of random polygonal maps and have found very littl

[R] piecewise nls?

I am looking into fitting a so-called double von Bertalanffy function to fish length-at-age data. Attempting to simplify the situation, the model looks like this ... Y ~ f(X; a,b,c) if x < Z Y ~ g(X; a,d,e) if x >= Z where * f and g are non-linear functions (the "traditional" "single" von Be

Re: [R] piecewise nls?

n Ritz [mailto:r...@life.ku.dk] > Sent: Saturday, April 17, 2010 2:45 PM > To: Derek Ogle > Cc: r-help@r-project.org > Subject: Re: [R] piecewise nls? > > Hi Derek, > > have a look at the following made-up example: > > f1 <- function(x){2*x} > f2 <- function(x

Re: [R] data frame's column names not the same as in CSV

Thank you all. I would like to say that it was my first time posting on the r-help mailing list. I am very impressed and grateful that I got the answer to my problem so quickly. Back to the issue, using read.csv (instead of read.csv) and turning off the check.names flag solved my problem. Derek

[R] Saving interactive modified latttice plots?

panel.identify() # left click on selected points. Right click to exit process trellis.unfocus() Derek N. Eder Gothenburg University "The most dangerous thing in the jungle is not the snakes, the spiders, the tigers, ... The most dangerous thing in the jungle is your

[R] Citing R/Packages Question

ould be interested in any comments or opinions. Dr. Derek H. Ogle Associate Professor of Mathematical Sciences and Natural Resources Northland College 1411 Ellis Avenue Box 112 Ashland, WI 715.682.1300 www.ncfaculty.net/dogle/ [[alterna

[R] interactive file choosing in Linux?

;text input prompt". This does not seem correct. > file.choose() Enter file name: I have seen gfile() function in the gWidgets library - but isn't there anything native? Thank you! Derek Eder Linux: Ubuntu 9.04, Gnome, [R] running in terminal or ESS

[R] Error: argument is of length zero

Hi, I have the following: for(j in (y.raw+1):(rownum-1)){ valsum<-tstfframed[min.x,j]+tstfframed[min.x,j+1] if(valsum == 1){ cat("valsum loop") int.num<-int.num+1 } }

[R] gpc.poly datatype

I have a list of polygons generated by the contourLines() command (each object of the list is a list in itself with two objects: a vector of x values, and a vector of y values for each vertex). I wish to convert that list into a gpc.poly object of multiple contours. How do I do this? gpclib app

[R] Coefficient of determination

Dear all, Is there a instruction that can help me obtain the coefficient of determination R^2 after doing linear/nonlinear regression using lm/nls? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/m

[R] driver file

nctions such as file(), pipe(), open(), found manuals on creating my own packages, etc. but I become increasingly aware that I don't know exactly what I am looking for (that I am a "neophyte"). Help would be appreciated Derek -- View this message in context: http://www.nabble.com/d

[R] help with formula and data= argument

Sorry for posting this twice, but I still have not solved this problem and am hoping for some assistance. I am attempting to write a function that is flexible enough to respond to the user providing a formula (with a data= argument) or not (similar to plot(x,y) versus plot(y~x,data=data)). I h

[R] plotCI (plotrix) problem

I am attempting to create a plot with intervals "stretched" in the x-direction using plotCI() in the plotrix package. The same data provides an appropriate set of intervals when "stretched" in the y-direction but I only get a lower interval when "stretched" in the x-direction. The data are as

Re: [R] Buglet in plotCI

Jim and Dieter, Thank you for the quick reply and precise fix. My plot now works as expected. Thanks again. > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of Jim Lemon > Sent: Monday, April 20, 2009 4:21 AM > To: dieter.me..

[R] RweaveHTML (R2HTML) Help

I have found Sweave() to be great for producing PDF documents. I have been experimenting with RweaveHTML (from the R2HTML) package and have had moderate success. My main issue has been that I simply want the R output to be shown verbatim in the HTML document but RweaveHTML tends to convert most o

[R] Suppress getOptions() in Stangle

Hello, I am creating PDF handouts using Sweave for a beginner's R workshop. I want to include the R code for the handout in an appendix as described (using the listings LaTeX package) as described in ... http://n4.nabble.com/including-Sweave-tangled-code-in-Rnw-document-tp875112p875112.html W

Re: [R] Inappropriate color name

I agree with, support, and appreciate this message. Thank you Ryan for the courage to write this, and Lainey for your courage to persist with your question. -Original Message- From: R-help On Behalf Of Ryan Novosielski Sent: Tuesday, November 17, 2020 9:40 AM To: r-help mailing list Cc

[R] RGB -> CYMK, with consistent colors

All, I used the very useful colorspace package for the plots in my book (pdf available here): http://knosof.co.uk/ESEUR/ The color makes the plots look great, on screen. To get lots printed, the printer requires converting the images to use cmyk (a common requirement for larger printers, I'm to

Re: [R] RGB -> CYMK, with consistent colors

Martin, Have you tried printed a few pages in CMYK? A monitor is based on mixing light using Red-Green-Blue. So it is not possible for the monitor to show CMYK which must be printed on paper to view correctly. Yes, I have printed some 'CMYK' pages. The blue is very obviously not cyan, as com

Re: [R] RGB -> CYMK, with consistent colors

Martin, However after some testing. I totally agree that CMYK handling in R using pdf(..., colormodel = "cmyk") is not correct. I thought the issue may be OS specific, but I get similar behavior on a Mac. I have discovered the Cyan tool: http://cyan.fxarena.net/ and will try it out. As vario

Re: [R] RGB -> CYMK, with consistent colors

Martin, However after some testing. I totally agree that CMYK handling in R using pdf(..., colormodel = "cmyk") is not correct. Looking at the source code of the function PostScriptSetCol in file src/librarygrDevices/src/devPS.c the conversion to CMYK looks correct. Looking at the generated p

[R] bgroup not rendering properly

bgroup() from plotmath does not render properly for me. For example plot(0,xlim=c(0,1),ylim=c(0,1)) text(0.3,0.5,expression(bgroup('(',atop(x,y),')'))) text(0.7,0.5,expression(group('(',atop(x,y),')'))) and library(ggplot2) ggplot(mtcars, aes(wt, mpg)) + annotate("text", x=2.5, y=25, label="bg

Re: [R] bgroup not rendering properly

Jinsong Zhao wrote: > > > On 2022/12/13 18:02, Martin Maechler wrote: > >>>>>> Jinsong Zhao > >>>>>> on Tue, 13 Dec 2022 17:07:00 +0800 writes: > > > I don > > > > > On 2022/12/13 10:13, Derek Ogle wrote

[R] coplot question

=19,col="red",cex=1.5,given.values=given.depth)| The resulting plot however has has bars in the conditional variable that don't cover all of the values of x and y at x=4, 7=14, and z=134 Thanks, Derek <>__ R-help@r-projec

Re: [R] coplot question

Hi José Yes indeed that does help. Thanks, Derek On 2/7/2013 5:01 AM, Jose Iparraguirre wrote: Hi, Derek, Try and modify the value for the overlap argument. Hope this helps, José -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf

[R] Need help with installing gplots

Hi, I am new to linux so please bear with me. OS is CentOS 5.9 - This cannot be changed I am following a guide given to me to setup a server. I am told to do the following: To Install R: sudo yum install gcc sudo yum install make sudo yum install telnet sudo rpm -Uvh http://download.fedorapr

Re: [R] Need help with installing gplots

re? What version of gplots will work with version 2.15.2 from the archive at http://cran.r-project.org/src/contrib/Archive/gplots/ On 13 June 2013 03:19, Prof Brian Ripley wrote: > On 13/06/2013 05:09, Derek Serianni wrote: > >> Hi, >> >> I am new to linux so please be

[R] optim() not finding optimal values

I am trying to use optim() to minimize a sum-of-squared deviations function based upon four parameters. The basic function is defined as ... SPsse <- function(par,B,CPE,SSE.only=TRUE) { n <- length(B) # get number of years of data B0 <- par["B0"]

Re: [R] optim() not finding optimal values

11:52 PM > To: Ravi Varadhan > Cc: Derek Ogle; R (r-help@R-project.org) > Subject: Re: [R] optim() not finding optimal values > > A slightly better scaling is the following: > > par.scale <- c(1.e06, 1.e06, 1.e-05, 1) # "q" is scaled differently > > >

Re: [R] when all I have is a contingency table....

A.Z., You could recreate the data (I assumed that the values are discrete integers and not a bin for a continuous variable) ... value <- 1:6 count <- c(10,8,12,9,14,7) d <- rep(value,count) table(d) ... and then do what you want with the data in the d vector ... summary(d) hist(d) > On 07-J

[R] Repeating values in a list

I have a list that looks like this ... > have <- list(a=7,b=3,c=1) > have $a [1] 7 $b [1] 3 $c [1] 1 and I want to have a simple way to change it to the following without re-typing the values ... > desire <- list(a=c(7,7),b=c(3,3),c=c(1,1)) > desire $a [1] 7 7 $b [1] 3 3 $c [1] 1 1 In othe

Re: [R] Repeating values in a list

e,c(1,1,2)) $a [1] 7 $b [1] 3 $c [1] 1 1 From: Dennis Murphy [mailto:djmu...@gmail.com] Sent: Saturday, September 18, 2010 8:25 PM To: Derek Ogle Cc: R (r-help@R-project.org) Subject: Re: [R] Repeating values in a list Hi: How about > have <- list(a=7,b=3,c=1) > lapply(have, rep, 2

Re: [R] fitting model to resampled data

Mike, Without completely knowing your end game with these questions and this procedure it does seem like you are "re-inventing the wheel" here. If that is true and given the nls() fit that you are using I would suggest that you look at boot.case() in the alr3 package or nlsBoot() in the nlstoo

[R] Error: cannot allocate vector of size x Gb (64-bit ... yet again)

.7 606849099 4629.9 550124408 4197.2 > matrix(x, ncol=16) Error: cannot allocate vector of size 3.7 Gb I don't understand how this operation violates the limits detailed in the Memory-limit help (above). Thank you! Derek Eder -

[R] Function to Define a Function

I am trying to define a general R function that has a function as the output that depends on the user's input arguments (this may make more sense by looking at the toy example below). My real use for this type of code is to allow a user to choose from many parameterizations of the same general

Re: [R] Function to Define a Function

Gabor ... that worked perfectly. Thank you. > -Original Message- > From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com] > Sent: Monday, August 09, 2010 10:20 PM > To: Derek Ogle > Cc: R (r-help@R-project.org) > Subject: Re: [R] Function to Define a Function >

Re: [R] Linear regression on several groups

Have a look at lmList() in the nlme package. > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of JesperHybel > Sent: Friday, August 13, 2010 7:56 AM > To: r-help@r-project.org > Subject: Re: [R] Linear regression on several groups

Re: [R] Equality of Vectors

Patrick, See all(). For example, > all(c(1,2,3)==c(1,2,3)) [1] TRUE > all(c(1,2,3)==c(2,1,3)) [1] FALSE > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of Downey, Patrick > Sent: Friday, August 13, 2010 1:49 PM > To: r-help@r

[R] Help Preparing Data for Heatmap creation

in/R # rdf: # dc:title textingHeatMap.R # dc:date 2010.05.14 # dc:creator http://www.mm.com/user/djdees/knows/who#derek-dees # dc:language R # dc:rights Copyright © # dc:description Creates a heatmap showing # of texts per hour per day. # doap:SVNRepository # doap:browse # doap:homepage # doap:wiki

[R] calculating "treatment effects" (differences) in a data frame?

I am overlooking a very simple function, something along the lines of tapply(). Thank you for you attention! Derek Eder ## Examples: myData = data.frame( ID = c("a","a","b","c","c","d","d"), visit=c(1,2,1,1,2,1,2), sco

[R] Finding Imported Packages

Is there a function similar to pkgDepends() that returns the packages that a particular package imports or imports from? I believe I can get this information from the matrix returned from installed.packages() but something like pkgDepends() would be more convenient. I did not find anything use

[R] Biplot - how can I make it more readable?

including outlier removal from a few of the variables, which has made the plot better, however it still is too cluttered to read the variable names. Is there any further improvment that can be made to the attached plot and if so what steps should I take? Thanks, Derek

Re: [R] Graph digitisation / tracing

I have successfully used the digitize package for this purpose ... http://cran.r-project.org/web/packages/digitize/index.html wrote: > > I want to take some published graphs and digitise them to allow me to > run some analysis on them. > Is this possible using any of R's > plugins. I don't thi

[R] Distribution Graph Idea/Help

For pedagogical reasons (i.e., moving a way from tables of distributions) I want to plot a probability density or distribution function with a shaded area corresponding to the calculated value or user input. For example, I want a plot to visually demonstrate the result of pnorm(1.3) I have

[R] R-Forge vs. RForge

I have two packages that I use locally and have been curious about moving them to "R Forge" to take advantage of SVN, builds, etc. In looking into this I stumbled across two different sites. The apparent official "R Forge" site at http://r-forge.r-project.org/ and a seemingly "unofficial" sit

[R] Ordered Multidimensional Arrays

Any ideas? A "zoo"-like package capable of handling multidimensional arrays would be optimal, but I don't believe this exists. Thanks, Derek [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://st

[R] Table Modification

I am trying to construct a two-way table where, instead of printing the two-way frequencies in the table, I would like to print the values of a third variable that correspond to the frequencies. For example, the following is easily constructed in R > fact1 <- factor(sample(LETTERS[1:3],10,r

Re: [R] Table Modification

and accurate answer. -Original Message- From: r...@quantide.com [mailto:r...@quantide.com] Sent: Friday, January 23, 2009 9:25 AM To: Derek Ogle Cc: r-h...@stat.math.ethz.ch Subject: Re: [R] Table Modification If I understood propelly > tapply(fact3, list(fact1, fact2) , pa

Re: [R] batch process file in R

i,sep=""),read.table(paste(fpath,a[i],sep="/"))) } Another option is to put them in a list... for(i in 1:length(a)){ y[[i]] <- read.table(paste(fpath,a[i],sep="/") } or with lapply... y <- lapply(1:length(a), function(i) read.table(paste(fpath,a [i],sep="/")

[R] NLS -- multiplicative errors and group comparison

Hello, I am attempting to fit a non-linear model (Von Bertalanffy growth model) to fish length-at-age data with the purpose of determining if any of the three parameters differ between male and female fish. I believe that I can successfully accomplish this goal assuming an additive error struc

[R] seq() question

Hello all, I am attempting to use the seq() function to define breaks for the cut() function but am finding performance that I do not understand. The root of my problem appears to be contained in the following simple example. > brks <- seq(0.8,1.2,0.1) > brks [1] 0.8 0.9 1.0 1.1 1.2 > brks==1.

[R] memory.limit

I ran the memory limit function in R 2.9.0 and received the 'error' below. The memory appears to update correctly, so there's probably no implication beyond cosmetic; however, thought I would make sure since the function as written did not generate the same error in my 2.8.0 version of R. Thanks

Re: [R] memory.limit

tring R version 2.9.0 (2009-04-17) > > All: is this a bug or are we missing something? > > Regards, > > Francisco > > __ > Francisco J. Zagmutt > Vose Consulting > 2891 20th Street > Boulder, CO, 80304 > USA > franci...@voseconsulting.com > www.voseconsulting

Re: [R] RGB -> CYMK, with consistent colors

t cyan, as compared to printing the RGB version. -- Derek M. Jones Evidence-based software engineering tel: +44 (0)1252 520667 blog:shape-of-code.coding-guidelines.com __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,

[R] Multivariate binary response analysis

ure, or with any of the R packages that implement it)? Any other recommendations? (NB: just discovered the bayesm package. Don't know if this will work for me or not.) Any help would be greatly appreciated. Derek Janszen, PhD Sr Research Biostatistician Computational Biology & Bioinformatics Pac

[R] log y 'axis' of histogram

All, I have been trying to get calls to hist(...) to be plotted with the y-axis having a log scale. I have tried: par(ylog=TRUE) I have also looked at the histogram package. Suggestions welcome. -- Derek M. Jones tel: +44 (0) 1252 520 667 Knowledge Software Ltd

Re: [R] log y 'axis' of histogram

C. and Olbricht, W. (1992) Statistics of particle size data. Appl. Statist., 41, 127–146. Hadley On Sunday, August 29, 2010, Joshua Wiley wrote: Hi Derek, Here is an option using the package ggplot2: library(ggplot2) x<- sample(x = 10:50, size = 50, replace = TRUE) qplot(x = x, geom

Re: [R] log y 'axis' of histogram

;) + scale_y_log() does a good job of highlighting the peaks. It may be useful for your purposes, but that doesn't necessarily make it a meaningful graphic. Doesn't being useful for my purpose make it meaningful, at least for me and I hope my readers? -- Derek M. Jones

Re: [R] log y 'axis' of histogram

ch used grap to draw all the graphs www.lunabase.org/~faber/Vault/software/grap/ with the numbers being extracted and processed by various C programs and awk scripts. -- Derek M. Jones tel: +44 (0) 1252 520 667 Knowledge Software Ltd mailto:de...@knosof.co.