Re: [R] read/write connections

From the help: Not all modes are applicable to all connections: for example URLs can only be opened for reading. You omitted to tell us your platform (and requests to include that and not send HTML mail are in the posting guide), and pipe() differs by platform. But AFAIR pipe() is

Re: [R] regexpr virtue

Many thanks to you and all others for solutions together with regular expressions insights. Petr Chuck Taylor ctay...@tibco.com napsal dne 29.07.2009 19:45:10: Petr, Here is a different approach. It, in effect, works from the end of the string backwards, rather than from the beginning of

Re: [R] Similar package like SEM

Yes thank you, as far as I understood your proposal it would lead to a similar result. Is there a good package in order to do principal component analysis? Best wishes, Luba Von: Andreas Hary [mailto:andreash...@googlemail.com] Gesendet: Mittwoch, 29. Juli

[R] Grid lines

Hi everyone. I am new to R. It will be greatly appreciated if someone can help me with the following questions. Here's the graph I have just produced. http://www.nabble.com/file/p24732839/calibration.jpeg (1) How can I put the legend on top of the grid lines? (2) How can I match the grid

Re: [R] Which CRAN mirror is the fastest ?

Tal Galili wrote: Hi dear R help group. I was wondering how can I find out which CRAN mirror would be the fastest for me (I leave in Israel and therefore don't have a CRAN mirror for our country to go to by default) Hard to lee, you have to try out, I fear. The speed you see highly

Re: [R] Similar package like SEM

See ?princomp and ?prcomp Uwe Ligges Stein, Luba (AIM SE) wrote: Yes thank you, as far as I understood your proposal it would lead to a similar result. Is there a good package in order to do principal component analysis? Best wishes, Luba Von: Andreas

Re: [R] Adding picture to graph?

Thanks a lot for all your responses - I'll check them out and pass them on. Cheers Rainer On Wed, Jul 29, 2009 at 6:32 PM, Greg Snowgreg.s...@imail.org wrote: In addition to the other responses, you may want to look at the subplot function in the TeachingDemos package for a way to place the

Re: [R] Similar package like SEM

Thanks for all the answers! I have started already using the Kalman filter but I am still not finished because it takes some time to prepare the variables. I will also try princomp and prcomp. So thanks a lot! Best wishes, Luba -Ursprüngliche Nachricht- Von: Uwe Ligges

[R] lattice shingle plot axis annotation

Hello (R-)Experts I hope someone can help with this problem concerning axis annotation of a lattice shingle plot. I want a plot with three shingles to display some laboratory value over time. In the first panel over the first few days, then in the next panel some months, and in the last panel

[R] Another SEM question

Hello, I work with SEM the last days in order to fit given time series into a model and obtain latent variables. Now I am wondering why when multiplying all given values with 10 or 100 changes the statistic so much. In fact, it should make no difference for the regression if the factors have a

[R] stepwise variable selection method wanted

Hi List, I am looking for a variable selection procedure with a forward-backward selection method. Firstly, it is meant to work with the cophenetic correlation coefficient (CPCC) and intended to find the variable combination with the highest cophenetic correlation. Secondly, it is aimed at

Re: [R] Transporting data in different R windows

Patrick Connolly wrote: On Wed, 29-Jul-2009 at 06:15AM -0700, Bogaso wrote: | | Hi, | | I am working with multiple R-windows opened and would like to | transport some variables from one to another. Is there any direct | way to do this without saving them in hard-disk? With so little

Re: [R] Grid lines

On Wed, 29-Jul-2009 at 11:32PM -0700, Chris Li wrote: | | Hi everyone. I am new to R. | | It will be greatly appreciated if someone can help me with the following | questions. | | Here's the graph I have just produced. | http://www.nabble.com/file/p24732839/calibration.jpeg | | (1) How can

Re: [R] Grid lines

Chris Li wrote: Hi everyone. I am new to R. It will be greatly appreciated if someone can help me with the following questions. Here's the graph I have just produced. http://www.nabble.com/file/p24732839/calibration.jpeg (1) How can I put the legend on top of the grid lines? It

Re: [R] Transparency and pyramid.plot

ws wrote: Hi all Another question related to my earlier -- how might I get transparency to work with pyramid.plot in library(plotrix)? I can't figure it out, and the help page says it is there but doesn't give an example. Works for me just by specifying colors with some alpha value as

Re: [R] Grid lines

Maybe this helps par(mar = c(5,4,6,1)) plot(1:10, type = n, axes = F) xticks = axis(1, at = c(1,5,9), labels = c(1-1-2009,1-5-2009,1-9-2009)) #set axis manually yticks = axis(2, at = 1:10) #set axis manually abline(v = xticks, col = gray) #verticall lines abline(h = yticks, col = gray,

[R] What is the best method to produce means by categorical factors?

I am attempting to replicate some of my experience from SAS in R and assume there are best methods for using a combination of summary(), subset, and which() to produce a subset of mean values by categorical or ordinal factors. within sas I would write proc means mean data=dataset; class factor1

[R] randomized block design analysis PROBLEM

Dear All user, Hello, I'm a student and I have some trouble with the experimental (columns-experiments) design of my project. I use a randomized block design with 4 treatments including a control. For each treatment, I use 3 replicates and 3 blocks. The treatments are: -T1 = COD (300 mg/Lit)

Re: [R] Grid lines

Maybe this helps par(mar = c(5,4,6,1)) plot(1:10, type = n, axes = F) xticks = axis(1, at = c(1,5,9), labels = c(1-1-2009,1-5-2009,1-9-2009)) #set axis manually yticks = axis(2, at = 1:10) #set axis manually abline(v = xticks, col = gray) #verticall lines abline(h = yticks, col = gray,

[R] Odp: What is the best method to produce means by categorical factors?

Hi r-help-boun...@r-project.org napsal dne 30.07.2009 10:19:21: I am attempting to replicate some of my experience from SAS in R and assume there are best methods for using a combination of summary(), subset, and which() to produce a subset of mean values by categorical or ordinal factors.

Re: [R] Grid lines

Thanks for the reply guys. Below are the script and data files. http://www.nabble.com/file/p24734049/HEAD_calibration.R HEAD_calibration.R http://www.nabble.com/file/p24734049/6527-1052 6527-1052 http://www.nabble.com/file/p24734049/6527-1054 6527-1054

Re: [R] What is the best method to produce means by categorical factors?

Dear Pat, Have a look at recast from the reshape package. library(reshape) dataset - expand.grid(factor1 = c(A, B), factor2 = c(C, D), Rep = 1:3) dataset$variable1 - rnorm(nrow(dataset)) dataset$variable2 - rnorm(nrow(dataset), mean = 10) recast(factor1 + factor2 + variable ~ ., data = dataset,

[R] Problem with RGtk2 Rattle

HI Apologies for previously trying to post this question onto the Dev forum. I have recently update my versions of R and related packages. When I try to use rattle the following message appears Error in .RGtkCall(R_setGObjectProps, obj, value, PACKAGE = RGtk2) : Invalid property

Re: [R] Which CRAN mirror is the fastest ?

2009/7/30 Uwe Ligges lig...@statistik.tu-dortmund.de: Hard to lee, you have to try out, I fear. The speed you see highly depends on the connection from your country to others, but of course, there are also some mirrors that are not the fastest themselves. I figured you could write a

[R] PDF Compression

Hi, I am generating a large number of graphs with pdf() and incorporating them in pdf document using pdflatex. According to the pdf() help: 'pdf' writes uncompressed PDF. It is primarily intended for producing PDF graphics for inclusion in other documents, and PDF-includers such as

Re: [R] check for new files in a given directory

William Dunlap wdun...@tibco.com on Tue, 28 Jul 2009 08:50:29 -0700 writes: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Barry Rowlingson Sent: Tuesday, July 28, 2009 6:41 AM To:

Re: [R] PDF Compression

Hi, You could have a look at exporting the pdfs using the Cairo package. cheers, Paul David Keegan schreef: Hi, I am generating a large number of graphs with pdf() and incorporating them in pdf document using pdflatex. According to the pdf() help: 'pdf' writes uncompressed PDF. It is

Re: [R] PDF Compression

After generating the pdf try this using the free pdftk utilty: pdftk infile.pdf output outfile.pdf compress On Thu, Jul 30, 2009 at 5:13 AM, David Keegandavid.kee...@shenick.com wrote: Hi, I am generating a large number of graphs with pdf() and incorporating them in pdf document using

Re: [R] Which CRAN mirror is the fastest ?

Barry Rowlingson b.rowling...@lancaster.ac.uk on Thu, 30 Jul 2009 09:59:47 +0100 writes: 2009/7/30 Uwe Ligges lig...@statistik.tu-dortmund.de: Hard to lee, you have to try out, I fear. The speed you see highly depends on the connection from your country to others,

Re: [R] Function match.call or saving all initial parameters

This does what you want: expand.call - function(definition = NULL, call = sys.call(sys.parent()), expand.dots = TRUE) # like match.call, but with all formal args instead of only the specified ones { ans - as.list(match.call(definition, call, expand.dots)) frmls -

Re: [R] Function match.call or saving all initial parameters

it doesn't actually- this should work: expand.call - function(call=sys.call(sys.parent(1)), expand.dots = TRUE) # similar to match.call, but with all formal args not the specified ones only { ans - as.list(call) frmls - formals(deparse(ans[[1]])) add - which(!(names(frmls) %in%

Re: [R] Which CRAN mirror is the fastest ?

Barry - thank you, that is exactly what I had in mind (but had no knowledge on how to do) Martin - I am glad to see your interest in this. Should I now start following the R-devel so to know how this is being handled ? Barry and Martin, your posts made me smile, thank you both :) Tal On

Re: [R] PDF Compression

Gabor, Thanks for the suggestion. I tried pdftk but it made very little difference. Regards, David. -- [David Keegan david.kee...@shenick.com 353 1 2710818] Gabor Grothendieck writes: After generating the pdf try this using the free pdftk utilty: pdftk infile.pdf output outfile.pdf

Re: [R] Which CRAN mirror is the fastest ?

There may also be a difference in reliability, which would not so easily be measured by an individual user. I've selected the closest geographically until it seemed to be down, then tried the second closest, etc. This could be automated centrally, but then you'd have to deal with the

Re: [R] Grid lines

Chris Li wrote: Thanks for the reply guys. Below are the script and data files. http://www.nabble.com/file/p24734049/HEAD_calibration.R HEAD_calibration.R http://www.nabble.com/file/p24734049/6527-1052 6527-1052 http://www.nabble.com/file/p24734049/6527-1054 6527-1054

[R] creating subsets within lm()

Dear All, the lm() function has the possibility to create a subset of the possible explaining variables that you have. However, in the help there is no example how to use this subset option. I tried the following:

Re: [R] PDF Compression

I just tried it with a recent pdf that was generated from R on Windows Vista with R version 2.9.1 Patched (2009-07-16 r48939). This particular one was laden with many graphs and was reduced to 25% of the original size so my experience with that one was that it made a huge difference. On Thu, Jul

[R] Moving Average on Defined Intervals

Hi all I have been looking (in the help archives) for a function which does a moving average. Nothing new, I know. But I am looking for a function which is very flexible: The user should be able to input a vector of breaks which define the bins (and not just the number of observations in a bin).

[R] temporal and spatial pseudorreplication with lmer

Dear all, I am trying to fit a generalized linear mixed model to deal with both temporal and spatial pseudorreplication. I have repeated seasonal measurements (3 seasons during 3 years and 2 season during the last year, factor named percod, with 11 levels) of a bird presence/absence

Re: [R] creating subsets within lm()

Thomas A. Groen wrote: Dear All, the lm() function has the possibility to create a subset of the possible explaining variables that you have. However, in the help there is no example how to use this subset option. I tried the following:

Re: [R] PDF Compression

Gabor, Interesting. I have had success compressing the pdf generated by R with pdftk. However when I incorporate that in a pdflatex document it ends up decompressed in the final pdf. Regards, David. -- [David Keegan david.kee...@shenick.com 353 1 2710818] Gabor Grothendieck writes: I just

[R] warnings and errors

Hi, I am using R for a while now, but I have not yet been able to convince R to produce helpful warnings or error messages. If you run a script that sources other scripts it may take a while to understand which script produced an error or a warning. For example, R tells you that warnings exist

Re: [R] Adding picture to graph?

Rainer M Krug-6 wrote: Hi while teaching R, the question came up if it would be possible to add a picture (saved on the HDD) to a graph (generated by plot()), which we could not answer. It might easily kill a clean graph, but: is there a way of doing this, even one should not do

[R] Moving Average on Defined Intervals

Hi all I have been looking (in the help archives) for a function which does a moving average. Nothing new, I know. But I am looking for a function which is very flexible: The user should be able to input a vector of breaks which define the bins (and not just the number of observations in a bin).

Re: [R] warnings and errors

Place a cat statement at the top and bottom of each script. You can force warnings to be treated like errors using options(warn = 2). options(warn = 0) turns it back again. On Thu, Jul 30, 2009 at 6:30 AM, Jan M. Wienerm...@jan-wiener.net wrote: Hi, I am using R for a while now, but I have

Re: [R] Moving Average on Defined Intervals

Hadassa Brunschwig-2 wrote: Hi all I have been looking (in the help archives) for a function which does a moving average. Nothing new, I know. But I am looking for a function which is very flexible: The user should be able to input a vector of breaks which define the bins (and not just

Re: [R] Problem with RGtk2 Rattle

Hi Wayne, I have seen reports of that error at times on MS/Windows. This beta version should have fixed it: install.packages(rattle, repos=http://rattle.togaware.com;) Let me know how you go and I'll upload the new version to CRAN. It would also be useful to know what version or Rattle and

Re: [R] Moving Average on Defined Intervals

See ?cut and ?tapply for your first problem. Those together with ?rollappy in the zoo package can likely handle your second problem. In the future please read the last line to every message on r-help before posting. On Thu, Jul 30, 2009 at 6:15 AM, Hadassa

Re: [R] Rattle: Trouble opening .csv files

Hi James, I'm assuming this might be on MS/Windows? Could you try: install.packages(rattle, repos=http://rattle.togaware.com;) This beta may have fixed that problem. Let me know how you go. Regards, Graham 2009/6/30 Thompson, James james.thomp...@penfed.org I am a novice R user and

Re: [R] PDF Compression

Thanks for this advice - and to second Gabor's experience: I just tried it on a couple of my files and achieved reductions on the order of 90% (28.6MB to 4.4MB and 1.5GB to 170MB)! This file contains lots of small plots but also many scattersmooth()- images ... So I think it does quite well

[R] Relative paths in R?

Dear list, I use the R CMD BATCH to run a file of R commands (from php for info) This file of R commands contains a line to source a configuration file source('path/conf.R') In this configuration file I define some variables, e.g. the path of some file (uploaded before this R CMD

Re: [R] Problem with RGtk2 Rattle

HI Graham Thanks for responding so promptly - unfortunately downloading and running this new version of Rattle did not alter the outcome - I am however running on Windows XP Regards Wayne Wayne Murray wrote: HI Apologies for previously trying to post this question onto the Dev forum.

[R] Selecting Bootstrap Method for Quantile Regression

The help page and vignette for summary.rq(quantreg) mention that there are three different bootstrap methods available for the se=bootstrap argument, but I can't figure out how to select a particular method. For example, if I want to use the xy-pair bootstrap how do I indicate this in summary.rq?

Re: [R] pyramid.plot: x-axis scale

ws wrote: Hi R-ers Does anyone know how if there is a way to force pyramid plot where you determine the bottom axis range? I want to compare many different pyramids, so I am hoping to force the bottom axis range from 0% to 10%. Here is the code so far, if it matters: mk.pyr =

Re: [R] R User Group listings

There are now several R geographic user groups, and a few have mailing lists on the R mailing list system. Thanks to Martin M, there's also a pointer to a page I'm maintaining to list/describe the groups. The page is at http://macnash.telfer.uottawa.ca/RUG.html Contact me if you have a

[R] package reldist version 1.5-5.1 : how to not just compare two years's wage distribution

Dear all users, I read all R programs and texts in Relative Distribution Methods website, but I encounter two problems when I used it. then I used Google to try to find any solutions about the two problems,but I got no useful hints. 1、 In package reldist y mean sample from comparison

[R] reading a file

Hello there, I have the following info in MyData: 01 0.40 02 0.40 03 0.40 04 0.35 05 0.34 06 0.33 I want to read them using the following code: ifn - MyData dat - read.table(ifn) names(dat1)-c(Code,M) names(dat2)-c(Code,M) I want to have

[R] Relative paths in R?

Dear list, I use the R CMD BATCH to run a file of R commands (from php for info) This file of R commands contains a line to source a configuration file source('path/conf.R') In this configuration file I define some variables, e.g. the path of some file (uploaded before this R CMD using

Re: [R] reading a file

Hi, On Jul 30, 2009, at 9:08 AM, Mohsen Jafarikia wrote: Hello there, I have the following info in MyData: 01 0.40 02 0.40 03 0.40 04 0.35 05 0.34 06 0.33 I want to read them using the following code: ifn - MyData dat - read.table(ifn)

[R] Loess question

Hi, suppose I have a vector of values in 'A'. Suppose I use: scatter.smooth(A) to add in a smooth curve. Is there a way to obtain the maximum value of the smoothed curve? Anyways, any help will be greatly appreciated! -- View this message in context:

Re: [R] Relative paths in R?

Hi, (Sorry, noticed I didn't CC r-help originally) On Jul 30, 2009, at 8:16 AM, julien cuisinier wrote: Dear list, I use the R CMD BATCH to run a file of R commands (from php for info) This file of R commands contains a line to source a configuration file source('path/conf.R') In this

Re: [R] Loess question

Hi, On Jul 30, 2009, at 9:35 AM, njhuang86 wrote: Hi, suppose I have a vector of values in 'A'. Suppose I use: scatter.smooth(A) to add in a smooth curve. Is there a way to obtain the maximum value of the smoothed curve? Anyways, any help will be greatly appreciated! Look at the help for

Re: [R] RODBC and Oracle

I just joined and though I did quickly read the posting guide it is quite possible not going to be a perfect posting. For one I tried to figure out how to reply to an existing message on the topic. I gave up but figured I keep at least the same subject. In a gist I get an empty data frame

Re: [R] Re lative paths in R?

Maybe if you could send that path via the arguments, you could get them back using commandArgs(). Then you can setwd() and paste() to generate paths. julien cuisinier wrote: Dear list, I use the R CMD BATCH to run a file of R commands (from php for info) This file of R

Re: [R] Relative paths in R?

It does obviously - very silly of me to ask PathUpload - paste(getwd(),/uploads,sep=) Many thanks Julien CC: r-help@r-project.org From: mailinglist.honey...@gmail.com To: j_cuisin...@hotmail.com Subject: Re: [R] Relative paths in R? Date: Thu, 30 Jul 2009 09:41:55 -0400 Hi,

Re: [R] Which CRAN mirror is the fastest ?

On Thu, Jul 30, 2009 at 10:49 AM, Martin Maechlermaech...@stat.math.ethz.ch wrote: Yes!  And visualize the corresponding  nearest neigbourhood for each CRAN mirror on a world map and make this dynamically refreshing every few minutes and put it on a webserver so people can watch the CRAN

Re: [R] PDF Compression

PDF created by R is in vector format. If you really want smaller files you can try creating PNGs instead. With a high enough resolution (e.g. dpi=600) there won't be much difference in the printed version of your document. -eduardo On Thu, Jul 30, 2009 at 7:55 AM, Benno

[R] R 2.9.1: Error building target 'front'

Hi, I am getting the following error # (R-2.9.1/src/library/methods): gmake front /bin/sh: 0403-057 Syntax error at line 1 : `;' is not expected. gmake: *** [front] Error 2 However, I am not sure what script sh is trying to execute at this point. I am building on AIX 5.3.

Re: [R] PDF Compression

Eduardo Leoni wrote: PDF created by R is in vector format. If you really want smaller files you can try creating PNGs instead. With a high enough resolution (e.g. dpi=600) there won't be much difference in the printed version of your document. -eduardo Eduardo, It is not a good idea in

[R] dynamic network analysis

Dear group, I am browsing the web to find suitable software for dynamic network analysis. I came across the sna-package under R, but this seems to do static analysis only. Is this right? Or is there a separate package for dynamic network analysis? Thanks

[R] package lattice can't be loeaded

I try to load the package lattice by library(lattice). Although I have installed it through CRAN by internet as well as from a local directory using the zip-file, it can't be loaded into the editor. Can anyone give advice, please. I am using R 2.9.1 on Widows XP. Thanks a lot, Luba

Re: [R] R User Group listings

On Thu, Jul 30, 2009 at 2:33 PM, Prof. John C Nashnas...@uottawa.ca wrote: There are now several R geographic user groups, and a few have mailing lists on the R mailing list system. Thanks to Martin M, there's also a pointer to a page I'm maintaining to list/describe the groups. The page is at

[R] (no subject)

I try to load the package lattice by library(lattice). Although I have installed it through CRAN by internet as well as from a local directory using the zip-file, it can't be loaded into the editor. Can anyone give advice, please. I am using R 2.9.1 on Widows XP. Thanks a lot, Luba

Re: [R] PDF Compression

There was a previous post about this also: http://tolstoy.newcastle.edu.au/R/e2/help/07/05/17475.html I was able to use bitmap(,type=pdfwrite) on Ubuntu Linux but had trouble getting it to work on Windows. For now, it's kind of kludgey but I have in my .Rprofle (or .Rprofile.site on Windows):

Re: [R] Superscripts and rounding

I'd need access to the whole data file. I tried exporting some data from the website for it but it got too complex for me! (oh right you need the HYDAT package - sorry about that) =) However, it seems to me you have two chunks of data: you got it!!! stn[flow] - which has daily flow data in

[R] Plotting two overlapping curves

Hi. Sorry if this question might have already been asked... I'm trying to plot two curves on the same X-Y coordinate. The par(mfrow) command allows me to plot the two curves on the same window however on different coordinate systems. Is there a way to force the two curves into the same coordinate

Re: [R] Plotting two overlapping curves

one way is the following: x - seq(-2, 2, len = 11) y1 - 2*x^2 + 3*x + 1 y2 - - 2*x^2 - 3*x + 1 plot(x, y1, type = l, col = red, ylim = range(y1, y2)) lines(x, y2, col = blue) I hope it helps. Best, Dimitris njhuang86 wrote: Hi. Sorry if this question might have already been asked... I'm

Re: [R] Object equality for S4 objects

Hi Stavros -- Stavros Macrakis macra...@alum.mit.edu writes: To test two environments for object equality (Lisp EQ), I can use 'identity': e1 - environment(local(function()x)) e2 - environment(local(function()x)) identical(e1,e2) # compares object identity [1] FALSE

[R] Looping through R objects with $ operator and tags

Hi all, Suppose I want to set the values in a column to the log of the values of another column like so: object$LogDist - log10(object$Distance) How do I loop through the objects if I have object1, object2, etc to perform this function? object1$LogDist - log10(object1$Distance) object2$LogDist

Re: [R] USGS stream flow data automatic download R

An Excel Web query will also do what you want. From an Excel Workbook create a web query and link it to the USGS website. You can schedule downloads or the query will download data up-to-date everytime the workbook is open. Felipe D. Carrillo Supervisory Fishery Biologist Department of

Re: [R] Looping through R objects with $ operator and tags

I would first put the objects in a 'list' since it is easier to work with and then use 'lapply': obj1 - obj2 - obj3 - data.frame(Distance=1:10) # create a list - earier to work with list.obj - list(obj1, obj2, obj3) list.obj - lapply(list.obj, function(x){ + x$LogDist - log(x$Distance) +

Re: [R] (no subject)

Exactly why are you trying to load it into the editor? What happens when you execute 'library(lattice)' at the command line? Is there an error message associated with it? Here is what I get: library(lattice) It loads just fine. On Thu, Jul 30, 2009 at 11:08 AM, Stein, Luba (AIM

[R] CRAN + geography = Cranography

In response to an earlier discussion about finding the fastest/best/nearest CRAN site, I geocoded (using my geonames package) the locations of the CRAN sites returned by getCRANmirrors. I've mapped them here: http://www.maths.lancs.ac.uk/~rowlings/R/Cranography/ The geocoding may have got it

Re: [R] Looping through R objects with $ operator and tags

essentially the same, object1 = object2 = object3 = data.frame(Distance = 1:10) foo = function(o){ within(get(o), LogDist - log10(Distance)) } my.objects = paste(object, 1:3,sep=) lapply(my.objects, foo) It is however advisable to group the initial objects in

Re: [R] Re lative paths in R?

See also file.path() and create.dir(). You can use them with getwd() and setwd() to specify and/or create subdirectories, relative to your current working directory. Handy because they allow you to create paths and directories with names derived within functions. For example, I have used them

[R] saving only part of a computation in integrate

Hello, I am using Beta -integrate(dnorm,mean=0,sd=1,-Inf,2.3552) Beta 0.9907436 with absolute error 4.2e-06 1-Beta results in an error How can I store into Beta only .9907436? Thank you. Sincerely, Mary A. Marion __ R-help@r-project.org mailing

Re: [R] saving only part of a computation in integrate

See ?str str(Beta) List of 5 $ value : num 0.99 $ abs.error : num 4.22e-06 $ subdivisions: int 4 $ message : chr OK $ call: language integrate(f = dnorm, lower = -Inf, upper = 2.3552, mean = 0, sd = 1) - attr(*, class)= chr integrate Beta$value [1] 0.9907436

Re: [R] saving only part of a computation in integrate

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Mary A. Marion Sent: Thursday, July 30, 2009 12:56 PM To: r-help@r-project.org Subject: [R] saving only part of a computation in integrate Hello, I am using Beta

[R] Continue to finish for loop even there is an error in one of rounds.

I am trying to fit a logistic model to my 10 year data (1999-2008) by year. Codes like below: Year - c(1999: 2008) for(y in 1:length(year)) { file.input - paste(C:\\, year[y], \\data.csv, sep=) table - read.csv(file=fileinput, header=TRUE, as.is=TRUE, na.strings=c()) initial -

Re: [R] lattice shingle plot axis annotation

Hi, Chapter 8 of the lattice book has some examples (you can see the code and figures on r-forge). Perhaps you could try something like this, d = data.frame(x=1:10,y=1:10,f=sample(letters[1:2],10,repl=T)) axis.custom = function(side, ...){ if(side == bottom)

Re: [R] R User Group listings

I've also been maintaining a list at: http://blog.revolution-computing.com/local-r-groups.html The most recent addition is the user group in Japan. # David Smith On Thu, Jul 30, 2009 at 5:33 AM, Prof. John C Nash nas...@uottawa.ca wrote: There are now several R geographic user groups, and a

[R] Question regarding sqldf

Hello, I am having a problem using sqldf. I'm trying to choose a subset of observations from a data set based on the date and maximum value of a variable by date. Here is the code I am using: test-sqldf(select distinct * from bextuse group by sdate having bext=max(bext),method=raw); The

[R] weight median by count for multiple records

Hello everyone, I have a .csv file with the following format: uniqueID SubjectID Distance_miles Tag 1 10015.5 3 2 10017 1 3 10016.5

Re: [R] weight median by count for multiple records

one approach is: sp - split(dat[-1], dat$SubjectID) t(sapply(sp, function (d) c(d$SubjectID[1], median(rep(d$Distance_miles, d$Tag) where 'dat' is the name of your data.frame. I hope this helps. Best, Dimitris Kirsten Beyer wrote: Hello everyone, I have a .csv file with the

Re: [R] Continue to finish for loop even there is an error in one of rounds.

You probably want to wrap the nls(...) in a 'try' or 'tryCatch' function. Scott Scott Sherrill-Mix Department of Microbiology University of Pennsylvania 402B Johnson Pavilion 3610 Hamilton Walk Philadelphia, PA 19104-6076 Phone: 215-573-3141 On Thu, Jul 30, 2009 at 1:14 PM, Liao,

Re: [R] Selecting Bootstrap Method for Quantile Regression

Hi Tom, For example, if I want to use the xy-pair bootstrap how do I indicate this in summary.rq? The general approach is documented under summary.rq (sub se option 5). Shorter route is boot.rq, where examples are given. ## ?boot.rq y - rnorm(50) x - matrix(rnorm(100),50) fit - rq(y~x,tau =

Re: [R] Object equality for S4 objects

On Thu, Jul 30, 2009 at 12:01 PM, Martin Morganmtmor...@fhcrc.org wrote: S4 objects do not have the semantics of environments, but of lists (or of most other R objects), so it is as meaningful to ask why identical(s1, s2) returns TRUE as it is to ask why identical(list(x=1), list(x=1))

[R] Retrieving original data frame after repetition

Dear R users, Consider the first two columns of a data frame like this: z[,1:2] x y 1 1 1 2 2 2 3 3 3 4 1 4 Imagine that y represents the times that the value x happens in a population. But z is not exactly a frequency table, because in z we have x=1 twice. So, the x=1 in the

[R] Scatter Plot

Dear Sir I want to write the numbers 1,2,3,on a scatter plot instead of points, like 1 corresponding to first point on plot, 2 corresponding second point etc. Help in this regard. Regards -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. [[alternative HTML

[R] Help in Rserve

Hello! I've been studying the Rserve for the last few days, and achieved to connect it to a Java plattaform. I also could run some simple commands in R, and plot some graphics. But I did this just using the examples files as a template. Now I'm looking for discover what does all those Java

  1   2   >