Re: [R] Evaluating a string

2010-07-02 Thread Joshua Wiley
Hello Jonathan, You can use eval() and parse(), although there may be better ways of getting to your end goal (assuming you have not inherited a series of character strings held in objects that should be run as commands). At any rate, this works: a - 3 b - 2 operator - - statement_string -

[R] Visualization of coefficients

2010-07-02 Thread Wincent
Dear all, I try to show a subset of coefficients in my presentation. It seems that a standard table is not a good way to go. I found figure 9 (page 9) in this file (

Re: [R] Good Package(s) for String and URL processing?

2010-07-02 Thread Tobias Verbeke
On 07/02/2010 05:51 AM, Erik Iverson wrote: Ralf B wrote: Are there packages that allow improved String and URL processing? E.g. extract parts of a URLs such as sub-domains, top-level domain, protocols (e.g. https, http, ftp), file type based on endings, check if a URL is valid or not, etc...

Re: [R] Kite diagrams

2010-07-02 Thread RCulloch
Hi Par, I am trying to do the exact same thing with my class, I would like to use R too, as well as get them to draw it out. I have tried to follow the suggestions but with no luck. If you did get round to sorting the code I wondered if you'd be so kind as to let me into the secret on how to do

Re: [R] Threshold Cointegration test

2010-07-02 Thread Setlhare Lekgatlhamang
Thanks David. I shall read the posting guide and resend. Sorry for being absent from the R users' discussion for a day - yesterday was a holiday in Botswana. Lexi -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Wednesday, June 30, 2010 5:24 PM To: Setlhare

Re: [R] Visualization of coefficients

2010-07-02 Thread Bernd Weiss
Am 02.07.2010 08:10, schrieb Wincent: Dear all, I try to show a subset of coefficients in my presentation. It seems that a standard table is not a good way to go. I found figure 9 (page 9) in this file (

[R] Why doesn't my change of Renviron.site work on windows 7 ?

2010-07-02 Thread Tal Galili
Hi all, I added the following line on the Renviron.site file: R_LIBS=C:\Program Files (x86)\R\library And when I start R and run: .libPaths() I don't see this path. On windows XP it worked for me. I am now using windoes 7 (64 bit) with R 32. Is there a reason this shouldn't work? Thanks,

Re: [R] Threshold Cointegration test

2010-07-02 Thread Setlhare Lekgatlhamang
Hi, I have read the posting guide and also the warning on my posting that was sent to me some minutes after I posted. I realize that my posting was too large/heavy coz of the equations in the attachment. It seems I will not be able to post it. But thanks a lot David. -Original Message-

Re: [R] Files with Missing Data

2010-07-02 Thread Setlhare Lekgatlhamang
Hi Norman, If you still have not managed, google for Econometrics in R pdf paper by Farnsworth; it as helped me a lot. You will only need to replace Farnsworth's read function with the read.csv function suggested below by Erik. Lexi -Original Message- From: r-help-boun...@r-project.org

Re: [R] problems loading the twitteR package

2010-07-02 Thread Allan Engelhardt
If you are on a Unix machine try your package manager first to see if it has the packages you need. Otherwise, try install.package(twitteR, dependencies = TRUE) (or RCurl) and see if that helps. Allan On 01/07/10 22:51, lyolya wrote: Dear all, I cannot load the twitteR package. When I

Re: [R] Files with Missing Data

2010-07-02 Thread Norman Jessup
Thanks Setlhare. everything is fine now I'm on the right track. regards Norman Jessup On 2/07/10 6:21 PM, Setlhare Lekgatlhamang wrote: Hi Norman, If you still have not managed, google for Econometrics in R pdf paper by Farnsworth; it as helped me a lot. You will only need to replace

Re: [R] Visualization of coefficients

2010-07-02 Thread Allan Engelhardt
http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=114 should get you started. On 02/07/10 07:10, Wincent wrote: Dear all, I try to show a subset of coefficients in my presentation. It seems that a standard table is not a good way to go. I found figure 9 (page 9) in this file (

[R] K-means result - variance between cluster

2010-07-02 Thread Ralph Modjesch
Hi, I like to present the results from the clustering method k-means in terms of variances: within and between Cluster. The k-means object gives only the within cluster sum of squares by cluster, so the between variance part is missing,for calculation the following table, which I try to get.

[R] possible to plot number line in R?

2010-07-02 Thread Kroepfl, Julia (julia.kroe...@uni-graz.at)
Thank you very much for your answers, but I think I did not explain thoroughly enough what I needed. I attached a demo of the plot. I need the number line between 2 and 3, both values being shown on the line, interval values should be printed next to the dashes and lines should connect the

Re: [R] Why doesn't my change of Renviron.site work on windows 7 ?

2010-07-02 Thread Joshua Wiley
Does Sys.getenv(R_LIBS) show the path? My guess is that it has to do with this little part from ?.libPaths Only directories which exist at the time will be included. Try entering the path you specified into R; when I do, I get an error that says Error: '\P' is an unrecognized escape in character

Re: [R] Why doesn't my change of Renviron.site work on windows 7 ?

2010-07-02 Thread Tal Galili
Thank you Joshua and Allan, I tried it now, and indeed it works. It remember for some reason that it didn't work for windows XP, was I wrong, or did something change ? Thanks again to the two of you. Best, Tal Contact

[R] Knowledge discovery

2010-07-02 Thread abanero
Hi, I have 10 units with 10 attributes (attr1, attr2, attr3, etc...) For instance: unit attr1 attr2 attr3 ... 1 a ww 12 2 a re 11 3 b ww 09 4 c yt 02 5 a qw 02 ... I'd like to answer to the question:

Re: [R] K-means result - variance between cluster

2010-07-02 Thread Christian Hennig
Dear Ralph, between and within clusters sum of squares (if you want variances, you need to divide them by the appropriate constant!) add up to the overall sum of squares, so you can get the beween clusters ss by computing the overall ss (one possibility to get this is to run kmeans with k=1)

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

2010-07-02 Thread Setlhare Lekgatlhamang
Try this Lm(y~X + I(X^2)), data=dd) # this runs OLS regression and it worked for me Hope it helps Lexi -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Otto Kässi Sent: Thursday, July 01, 2010 3:28 PM To: r-help@r-project.org

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

2010-07-02 Thread Peter Ehlers
On 2010-07-01 10:33, Allan Engelhardt wrote: I want my local libraries to have priority over the system installed ones, which, as far as I can make out from help(.libPaths), means they have to come first in that list (it doesn't actually_say_ so, but that seems to be the idea). We have

Re: [R] Good Package(s) for String and URL processing?

2010-07-02 Thread Gabor Grothendieck
On Thu, Jul 1, 2010 at 11:08 PM, Ralf B ralf.bie...@gmail.com wrote: Are there packages that allow improved String and URL processing? E.g. extract parts of a URLs such as sub-domains, top-level domain, protocols (e.g. https, http, ftp), file type based on endings, check if a URL is valid or

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

2010-07-02 Thread Otto Kässi
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, I should have mentioned this already in my original post. Nonetheless, thanks

[R] Apply a shift which is a function of the array.

2010-07-02 Thread Jim Hargreaves
Dear List, I have a 2,000x10,000 array of time domain data which when plotted draws a distinct pulse. The matrix is 10,000 pulses of length 2000. I would like the pulse to be shifted so that the peak (which.max of the pulse data) is consistently at point 400. I don't want to use loops as

Re: [R] Apply a shift which is a function of the array.

2010-07-02 Thread Robin Hankin
Hello Jim you can use ashift() from the same library which does (I think) what you want. HTH, Robin On 07/02/2010 12:05 PM, Jim Hargreaves wrote: Dear List, I have a 2,000x10,000 array of time domain data which when plotted draws a distinct pulse. The matrix is 10,000 pulses of

Re: [R] Apply a shift which is a function of the array.

2010-07-02 Thread Robin Hankin
Hello again Jim It seems that ashift() from the same package *doesn't* do what you want. But you can use shift() as follows: myshift - function(x){shift(x,1-which.max(x))} a - matrix(runif(30),5,6) a [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.268955362

[R] Combining several plots besides a dendrogram?

2010-07-02 Thread Tal Galili
Hello all, I would like to recreate the plot shown here (from a useR 2009 presentation): http://www.agrocampus-ouest.fr/math/useR-2009/abstracts/pdf/Hocking.pdf I downloaded the code for that image, and discovered that it relies on external web services, and also having PERL installed on the

[R] S4 classes and debugging - Is there a summary?

2010-07-02 Thread Joris Meys
Dear all, I'm getting more and more frustrated with the whole S4 thing and I'm looking for a more advanced summary on how to deal with them. Often I get error messages that don't make sense at all, or the code is not doing what I think it would do. Far too often inspecting the code requires me to

[R] metafor and meta-analysis at arm-level

2010-07-02 Thread Angelo Franchini
Hi, I have been looking for an R package which allowed to do meta-analysis (both pairwise and network/mixed-treatment) at arm-level rather than at trial-level, the latter being the common way in which meta-analysis is done. By arm-level meta-analysis I mean one that accounts for data provided at

Re: [R] Knowledge discovery

2010-07-02 Thread Carlos Ortega
Hello, For the first question, you can use table function. Read how to apply it, it is very straitghforward. Once you have the most common combination, you have to investigate further this combination to apply a statistical model to it. What model?. It depends on what are you looking for.

Re: [R] Evaluating a string

2010-07-02 Thread David Winsemius
On Jul 2, 2010, at 2:05 AM, Joshua Wiley wrote: Hello Jonathan, You can use eval() and parse(), although there may be better ways of getting to your end goal do.call(operator, list(a,b)) [1] 1 (assuming you have not inherited a series of character strings held in objects that should be

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

2010-07-02 Thread Peter Ehlers
Otto, The current version of ols() is fairly fussy about the way the predictors are used. I'm not fond of the I() construction anyway and so I would either use poly() or define a new predictor as you suggest in your original post. See also this thread:

Re: [R] Good Package(s) for String and URL processing?

2010-07-02 Thread David Winsemius
On Jul 1, 2010, at 11:08 PM, Ralf B wrote: Are there packages that allow improved String and URL processing? E.g. extract parts of a URLs such as sub-domains, top-level domain, protocols (e.g. https, http, ftp), file type based on endings, check if a URL is valid or not, etc... I am currently

Re: [R] K-means result - variance between cluster

2010-07-02 Thread Gabor Grothendieck
On Fri, Jul 2, 2010 at 4:37 AM, Ralph Modjesch ralph.modje...@pfeiffer-koberstein-immobilien.de wrote: Hi, I like to present the results from the clustering method k-means in terms of variances: within and between Cluster. The k-means object gives only the within cluster sum of squares by

Re: [R] Knowledge discovery

2010-07-02 Thread abanero
with table function you can just build a contigence table. What do you think about arules package? I thought mining associative rules is the correct approach to the problem.. Thanks Abanero -- View this message in context:

Re: [R] S4 classes and debugging - Is there a summary?

2010-07-02 Thread Joris Meys
Correction: trace off course works when specifying the signature. so eg: trace(xyValues,tracer=browser,signature=c(object=RasterLayer, xy=matrix)) allows me to browse through the function. Should have specified the signature, I overlooked that. Still, if there's a manual on S4 that anybody likes

Re: [R] Knowledge discovery

2010-07-02 Thread David Winsemius
On Jul 2, 2010, at 5:37 AM, abanero wrote: Hi, I have 10 units with 10 attributes (attr1, attr2, attr3, etc...) For instance: unit attr1 attr2 attr3 ... 1 a ww 12 2 a re 11 3 b ww 09 4 c yt 02 5 a qw

Re: [R] Modelling survival with time-dependent covariates

2010-07-02 Thread Terry Therneau
1-Would informing the algorithm coxph which samples represents the same person (through the use of an Id for example) improve the ?efficiency? of the estimated model? And if so, how should i do that? Using strata()? No, it makes no change. The reason is that the (start, stop] is just a trick.

Re: [R] Documentation for library() and how to specify missing arguments

2010-07-02 Thread Jeffrey J. Hallman
Gabor Grothendieck ggrothendi...@gmail.com writes: Try this: mx - formals(identity)$x missing(mx) [1] TRUE sin(mx) Error in sin(mx) : 'mx' is missing Neat. There's no way to look at 'mx' because calling any function with it as an argument gives that same error message, including the auto

Re: [R] anyone know why package RandomForest na.roughfix is so slow??

2010-07-02 Thread Liaw, Andy
I'll incorporate some of these ideas into the next release. Thanks! Best, Andy -Original Message- From: h.wick...@gmail.com [mailto:h.wick...@gmail.com] On Behalf Of Hadley Wickham Sent: Thursday, July 01, 2010 8:08 PM To: Mike Williamson Cc: Liaw, Andy; r-help Subject: Re: [R] anyone

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

2010-07-02 Thread Frank E Harrell Jr
ols has always been fussy about this. And don't use poly; use pol with rms/Design. Frank On 07/02/2010 07:28 AM, Peter Ehlers wrote: Otto, The current version of ols() is fairly fussy about the way the predictors are used. I'm not fond of the I() construction anyway and so I would either

Re: [R] Documentation for library() and how to specify missing arguments

2010-07-02 Thread Gabor Grothendieck
On Fri, Jul 2, 2010 at 9:27 AM, Jeffrey J. Hallman jhall...@frb.gov wrote: Gabor Grothendieck ggrothendi...@gmail.com writes: Try this: mx - formals(identity)$x missing(mx) [1] TRUE sin(mx) Error in sin(mx) : 'mx' is missing Neat. There's no way to look at 'mx' because calling any

Re: [R] Kite diagrams

2010-07-02 Thread Jim Lemon
On 07/02/2010 04:06 PM, RCulloch wrote: Hi Par, I am trying to do the exact same thing with my class, I would like to use R too, as well as get them to draw it out. I have tried to follow the suggestions but with no luck. If you did get round to sorting the code I wondered if you'd be so kind

[R] Some questions about R's modelling algebra

2010-07-02 Thread Hadley Wickham
Hi all, In preparation for teaching a class next week, I've been reviewing R's standard modelling algebra. I've used it for a long time and have a pretty good intuitive feel for how it works, but would like to understand more of the technical details. The best (online) reference I've found so far

Re: [R] S4 classes and debugging - Is there a summary?

2010-07-02 Thread Martin Morgan
On 07/02/2010 05:05 AM, Joris Meys wrote: Dear all, I'm getting more and more frustrated with the whole S4 thing and I'm looking for a more advanced summary on how to deal with them. Often I get error messages that don't make sense at all, or the code is not doing what I think it would do.

Re: [R] Some questions about R's modelling algebra

2010-07-02 Thread S Ellison
Hadley Wickham had...@rice.edu 02/07/2010 14:59:53 Where is %in% documented within R? I'm pretty sure it's a different action to ?%in%, and it's not mentioned in ?formula ?formula in R 2.9.2 says in para 2: The %in% operator indicates that the terms on its left are nested within those on

Re: [R] Some questions about R's modelling algebra

2010-07-02 Thread Hadley Wickham
?formula in R 2.9.2 says in para 2: The %in% operator indicates that the terms on its left are nested within those on the right. For example a + b %in% a expands to the formula a + a:b. Ooops, missed that. So b %in% a = a:b, and that's what's meant by different coding. Hadley --

Re: [R] Combining several plots besides a dendrogram?

2010-07-02 Thread Tal Galili
Hi Carlos, I was debating with myself if that is the strategy to go (or the one using subplots). Following some thread on the subject from lately, I am a bit hesitant (but maybe layout *is* the better strategy). http://r.789695.n4.nabble.com/plot-focus-tt2272699.html I am hoping for some more

Re: [R] S4 classes and debugging - Is there a summary?

2010-07-02 Thread Joris Meys
On Fri, Jul 2, 2010 at 4:01 PM, Martin Morgan mtmor...@fhcrc.org wrote:  selectMethod(xyValues, c('RasterLayer', 'matrix')) would be my choice. Thanks, that's a discovery. I guess I misunderstood the help pages on that one. I don't really have the right experience, but Chamber's 2008

Re: [R] Comparing and Interpreting GAMMs

2010-07-02 Thread Simon Wood
On Friday 28 May 2010 08:23, Andrea Meyer wrote: Dear R users I have a question related to the interpretation of results based on GAMMs using Simon Woods package gamm4. I have repeated measurements (hours24) of subjects (vpnr) and one factor with three levels (pred). The outcome (dv) is

[R] Problem with aggregating data across time points

2010-07-02 Thread Chris Beeley
Hello- I have a dataset which basically looks like this: Location Sex Date Time VerbalSelf harm Violence_objects Violence A 1 1-4-2007 1800 3 0 1 3 A 1 1-4-2007 1230

Re: [R] Some questions about R's modelling algebra

2010-07-02 Thread RICHARD M. HEIBERGER
Hadley, The S language modeling language was designed with Wilkinson and Rogers in mind. The notation was changed from their paper to retain consistency with the parsing rules for ordinary algebra in S. I think of : as an indicator of an indexing system into the dummy variables. It is not an

Re: [R] S4 classes and debugging - Is there a summary?

2010-07-02 Thread Martin Maechler
Martin Morgan mtmor...@fhcrc.org on Fri, 02 Jul 2010 07:01:54 -0700 writes: On 07/02/2010 05:05 AM, Joris Meys wrote: Dear all, I'm getting more and more frustrated with the whole S4 thing and I'm looking for a more advanced summary on how to deal with them. Often

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

2010-07-02 Thread Michael Friendly
I have two questions related to plotting predicted values for a linear mixed model using xyplot: 1: With a groups= argument, I can't seem to get the key to appear inside the xyplot. (I have the Lattice book, but don't find an example that actually does this.) 2: With lme(), how can I

Re: [R] parameterization of glm nested design

2010-07-02 Thread Charles C. Berry
Comments inline, below. On Thu, 1 Jul 2010, Huso, Manuela wrote: Chuck, Thank you for welcoming me to the list and thank you for taking the time to address my question. Pointing me to the pivot and qr components of my model object was very useful. But I still don't understand how R

[R] Odp: Problem with aggregating data across time points

2010-07-02 Thread Petr PIKAL
Hi did you try aggregate? aggregate(data[, 5:8],list(data$Date), sum, na.rm=T) Group.1 verbal self.harm violence_objects violence 1 0 000 2 01/04/07 251539 3 02/04/07 24 68 13

[R] merging plot labels in a lattice plot

2010-07-02 Thread Rajarshi Guha
Hi, I have a lattice lot conditioned on two variables. Example code is: library(lattice) x - data.frame(d=runif(100), f1=sample(c('yes', 'no'),100,replace=TRUE), f2=c(rep('Run1',30),rep('Run2',30),rep('Run3',40))) histogram(~d | f1 + f2, x) In the plot, for a

Re: [R] Double Integration

2010-07-02 Thread Christos Argyropoulos
Function adapt in package integrate maybe? Date: Thu, 1 Jul 2010 05:30:25 -0700 From: sarah_sanche...@yahoo.com To: r-help@r-project.org Subject: [R] Double Integration Dear R helpers I am working on the Bi-variate Normal distribution probabilities. I need to double integrate the

Re: [R] Why doesn't my change of Renviron.site work on windows 7 ?

2010-07-02 Thread Joshua Wiley
On Fri, Jul 2, 2010 at 2:14 AM, Tal Galili tal.gal...@gmail.com wrote: Thank you Joshua and Allan, I tried it now, and indeed it works. It remember for some reason that it didn't work for windows XP, was I wrong, or did something change ? I just tried it out on an XP machine and it did not

Re: [R] Combining several plots besides a dendrogram?

2010-07-02 Thread Carlos Ortega
Hello, Check for function layout. With it you can create separate plotting regions one for each graph. And these regions can be customized depending on the size of the objects you want to plot. Regards, Carlos. On Fri, Jul 2, 2010 at 1:34 PM, Tal Galili tal.gal...@gmail.com wrote: Hello all,

Re: [R] problems loading the twitteR package

2010-07-02 Thread lyolya
Ok thank you. Shall I re-write my forum entry (I mean, if it will be deleted as it is now)? On 2 July 2010 10:27, cybaea [via R] ml-node+2276136-1795582731-304...@n4.nabble.comml-node%2b2276136-1795582731-304...@n4.nabble.com wrote: If you are on a Unix machine try your package manager

Re: [R] Why doesn't my change of Renviron.site work on windows 7 ?

2010-07-02 Thread Bert Gunter
I think you need to enclose the path in quotes. -- Bert Bert Gunter Nonclinical Biostatistics 467-7374 http://devo.gene.com/groups/devo/depts/ncb/home.shtml -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Tal Galili Sent:

[R] GMM with covariance moment condicion

2010-07-02 Thread peterko
hello I have covariance stacionary proces, and i want to estimate some parameter of this proces via gmm. My problem is with write g -function. 0 order autocovariance is not problem 1 and higher order autocavariance are problem, because add order from 0 mean that I loose one observacion if I

[R] unable to get bigglm working, ATTN: Thomas Lumley

2010-07-02 Thread stephenb
I am using an example posted in this help forum to work with a file. the head of the file looks like: 97 2007-03-05 2007-06-01 90 3 5.450 205500.00 999.00 999.000 0.000 0 0 97 2007-03-06 2007-06-01 90 3 5.450 205500.00 999.00 999.000 0.000 1 0 97 2007-03-07 2007-06-01 90 3 5.450

Re: [R] Combining several plots besides a dendrogram?

2010-07-02 Thread Dylan Beaudette
Hi, The 'aqp' package originally used the subplot function to add images to to a dendrogram plot. I have since changed to use base graphics primitives, as the results tend to scale better. I avoided the use of 'layout' because sometimes it is convenient to add further embellishments that

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

2010-07-02 Thread Peter Ehlers
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] Problem with aggregating data across time points

2010-07-02 Thread Allan Engelhardt
On 02/07/10 16:21, Chris Beeley wrote: Hello- I have a dataset which basically looks like this: Location Sex Date Time VerbalSelf harm Violence_objects Violence A 1 1-4-2007 1800 3 0 1 3 A

Re: [R] Odp: Problem with aggregating data across time points

2010-07-02 Thread David Winsemius
On Jul 2, 2010, at 11:55 AM, Petr PIKAL wrote: Hi did you try aggregate? aggregate(data[, 5:8],list(data$Date), sum, na.rm=T) Group.1 verbal self.harm violence_objects violence 1 0 000 2 01/04/07 251539 3

Re: [R] Some questions about R's modelling algebra

2010-07-02 Thread Thomas Lumley
On Fri, 2 Jul 2010, Hadley Wickham wrote: Hi all, In preparation for teaching a class next week, I've been reviewing R's standard modelling algebra. I've used it for a long time and have a pretty good intuitive feel for how it works, but would like to understand more of the technical details.

Re: [R] unable to get bigglm working, ATTN: Thomas Lumley

2010-07-02 Thread Thomas Lumley
Actually, I think the problem *is* reading in the data If I try reading in your supplied lines of data with the read.table arguments() in your make.data() function I get your error message. data-read.table(tmp-textConnection( + 97 2007-03-05 2007-06-01 90 3 5.450 205500.00 999.00

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

2010-07-02 Thread Peter Ehlers
On 2010-07-02 9:37, Michael Friendly wrote: I have two questions related to plotting predicted values for a linear mixed model using xyplot: 1: With a groups= argument, I can't seem to get the key to appear inside the xyplot. (I have the Lattice book, but don't find an example that actually

Re: [R] Why doesn't my change of Renviron.site work on windows 7 ?

2010-07-02 Thread Tal Galili
Thanks Bert, I did and it worked. Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com (English)

Re: [R] Visualization of coefficients

2010-07-02 Thread Tal Galili
Specifically this link: http://tables2graphs.com/doku.php?id=04_regression_coefficients Great reference Bernd, thank you. Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com

Re: [R] Combining several plots besides a dendrogram?

2010-07-02 Thread Tal Galili
Wonderful Dylan, exactly the type of answer I was hoping for! :) (p.s: I would still be glad to read more replies if anyone has anything of value to add) Thanks, Tal Contact Details:--- Contact me: tal.gal...@gmail.com |

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

2010-07-02 Thread Peter Ehlers
On 2010-07-02 9:57, Rajarshi Guha wrote: Hi, I have a lattice lot conditioned on two variables. Example code is: library(lattice) x- data.frame(d=runif(100), f1=sample(c('yes', 'no'),100,replace=TRUE), f2=c(rep('Run1',30),rep('Run2',30),rep('Run3',40)))

[R] how to save summary(lm) and anova (lm) in format?

2010-07-02 Thread Yi
Hi, folks, I would like to copy the output of summary(lm) and anova (lm) in R to my word file. But the output will be a mess if I just copy after I call summary and anova. # x=rnorm(10) y=rnorm(10,mean=3) lm=lm(y~x) summary(lm) Call: lm(formula = y ~ x) Residuals: Min

Re: [R] Visualization of coefficients

2010-07-02 Thread Achim Zeileis
I've thought about adding a plot() method for the coeftest() function in the lmtest package. Essentially, it relies on a coef() and a vcov() method being available - and that a central limit theorem holds. For releasing it as a general function in the package the code is still too raw, but

Re: [R] how to save summary(lm) and anova (lm) in format?

2010-07-02 Thread David Winsemius
On Jul 2, 2010, at 1:34 PM, Yi wrote: Hi, folks, I would like to copy the output of summary(lm) and anova (lm) in R to my word file. But the output will be a mess if I just copy after I call summary and anova. A mess in what way exactly? # x=rnorm(10)

Re: [R] how to save summary(lm) and anova (lm) in format?

2010-07-02 Thread Yi
Hi, I attached a picture to compare the 'mess' to what I what. Maybe there is difference between computers? Thanks. On Fri, Jul 2, 2010 at 10:52 AM, David Winsemius dwinsem...@comcast.netwrote: On Jul 2, 2010, at 1:34 PM, Yi wrote: Hi, folks, I would like to copy the output of

Re: [R] how to save summary(lm) and anova (lm) in format?

2010-07-02 Thread David Winsemius
Nothing survived the trip through the list server. On Jul 2, 2010, at 1:59 PM, Yi wrote: Hi, I attached a picture to compare the 'mess' to what I what. Maybe there is difference between computers? Thanks. On Fri, Jul 2, 2010 at 10:52 AM, David Winsemius dwinsem...@comcast.net wrote:

Re: [R] Visualization of coefficients

2010-07-02 Thread Tal Galili
BTW, another visualization that might be useful in your case is Nomogramhttp://en.wikipedia.org/wiki/Nomogram : http://lib.stat.cmu.edu/S/Harrell/help/Design/html/nomogram.html (I remember first encountering it on a lecture by Frank Harrell lecture and being very happy for the discovery) Tal

Re: [R] Visualization of coefficients

2010-07-02 Thread Frank E Harrell Jr
Please note that Statlib is about 10 years out of date with respect to my software. See http://biostat.mc.vanderbilt.edu/Rrms Frank On 07/02/2010 01:12 PM, Tal Galili wrote: BTW, another visualization that might be useful in your case is Nomogramhttp://en.wikipedia.org/wiki/Nomogram :

Re: [R] how to save summary(lm) and anova (lm) in format?

2010-07-02 Thread Joshua Wiley
On Fri, Jul 2, 2010 at 10:59 AM, Yi liuyi.fe...@gmail.com wrote: Hi, I attached a picture to compare the 'mess' to what I what. Maybe there is difference between computers? Thanks. On Fri, Jul 2, 2010 at 10:52 AM, David Winsemius dwinsem...@comcast.netwrote: On Jul 2, 2010, at 1:34

[R] Producing residual plots by time for lme object

2010-07-02 Thread Gregory Gentlemen
Fellow R-users, I have a longitudinal data set with missing values in it. I would like to produce a residual plot for each time using panel.xyplot function but I get an error message. Here's a simple example, library(nlme) set.seed(1544) longdata - data.frame(ID=gl(10,1,50), y=rnorm(50), time

[R] finding subroutines

2010-07-02 Thread Erin Hodgess
Dear R People: I have found the starma.c program in R. But now I need to find the R_setup_starma and R_free_starma subroutines as well. Where would I go about finding them, please? Thanks for any help. Have a great weekend. Sincerely, Erin -- Erin Hodgess Associate Professor Department of

[R] editing/formatting stacked time series panels

2010-07-02 Thread Phil Morefield
Hello R-philes, I'm very new to R and hope this is a simple question to answer: I reading in 10 univariate time series and then using plot(), which produces two columns of stacked plots showing the data. Perfect. Next, I want to add things using lines(), but it's not working correctly. How do

Re: [R] Visualization of coefficients

2010-07-02 Thread Thomas Lumley
On Fri, 2 Jul 2010, Achim Zeileis wrote: I've thought about adding a plot() method for the coeftest() function in the lmtest package. Essentially, it relies on a coef() and a vcov() method being available - and that a central limit theorem holds. For releasing it as a general function in the

[R] generating list of all arguments that function was called with

2010-07-02 Thread Erik Iverson
Hello, Consider: f1 - function(a, b, c, d, ...) { c(list(a = a, b = b, c = c, d = d), list(...)) } f1(a = 1, b = 2, c = 3, d = 4, more = 5) $a [1] 1 $b [1] 2 $c [1] 3 $d [1] 4 $more [1] 5 Question: I'm guessing there exists a function such that I don't have to list each bound

Re: [R] Visualization of coefficients

2010-07-02 Thread Achim Zeileis
On Fri, 2 Jul 2010, Thomas Lumley wrote: On Fri, 2 Jul 2010, Achim Zeileis wrote: I've thought about adding a plot() method for the coeftest() function in the lmtest package. Essentially, it relies on a coef() and a vcov() method being available - and that a central limit theorem holds. For

Re: [R] unable to get bigglm working, ATTN: Thomas Lumley

2010-07-02 Thread stephenb
Sincere apologies and many thanks. The pasted code from the help forum has the error and I did not see it. When I wrote my own import I used c as I should. See: http://r.789695.n4.nabble.com/R-Example-function-for-bigglm-biglm-data-input-from-file-td816496.html#a816496 Stephen Bond | Senior

Re: [R] unable to get bigglm working, ATTN: Thomas Lumley

2010-07-02 Thread stephenb
after correcting the error spotted by Thomas, I tried again and it goes to work, but there is no result after 1 hour. is there anything I can do to debug? the dataset has 12mln rows and SAS on a server will produce results in 15 secs. I am running this on PC with XP (hence mem limit of 2.5 gigs)

Re: [R] generating list of all arguments that function was called with

2010-07-02 Thread Marc Schwartz
On Jul 2, 2010, at 3:01 PM, Erik Iverson wrote: Hello, Consider: f1 - function(a, b, c, d, ...) { c(list(a = a, b = b, c = c, d = d), list(...)) } f1(a = 1, b = 2, c = 3, d = 4, more = 5) $a [1] 1 $b [1] 2 $c [1] 3 $d [1] 4 $more [1] 5 Question: I'm guessing

Re: [R] Double Integration

2010-07-02 Thread Ravi Varadhan
There is no package called `integrate', but there is a function called `adaptIntegrate' in the cubature package. Ravi. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Christos Argyropoulos Sent: Friday, July 02, 2010 8:41 AM To:

Re: [R] Some questions about R's modelling algebra

2010-07-02 Thread Adrian Waddell
Hadley Wickham hadley at rice.edu Where is %in% documented within R? I'm pretty sure it's a different action to ?%in%, and it's not mentioned in ?formula You find the documentation for operators like -, %in%, if, etc by putting the operators between qoutes ?%in% ?- ?if Regards,

[R] strange data set output

2010-07-02 Thread Changbin Du
Hi, Dear Community, My data set logit.pred contains 2 columns and 1400 rows. When I want to use the first column, it is very strange. Where the $ come out? Thanks so much! dim(logit.pred) [1] 14002 head(logit.pred) tree.pred valid.out 754 0.6550606 1 1080 0.6353524

Re: [R] XML and RCurl: problem with encoding (htmlTreeParse)

2010-07-02 Thread Ryusuke Kenji
Hi All, First method:- library(XML) theurl - http://home.sina.com; download.file(theurl, tmp.html) txt - readLines(tmp.html) txt - htmlTreeParse(txt, error=function(...){}, useInternalNodes = TRUE) g - xpathSApply(txt, //p, function(x) xmlValue(x)) head(grep( , g, value=T)) [1] 繁體 |

Re: [R] strange data set output

2010-07-02 Thread David Winsemius
The str function would have given you a better one-shot look at logit.pred. -- David. On Jul 2, 2010, at 8:07 PM, Changbin Du wrote: Hi, Dear Community, My data set logit.pred contains 2 columns and 1400 rows. When I want to use the first column, it is very strange. Where the $ come

Re: [R] Double Integration

2010-07-02 Thread David Winsemius
And an adapt() in fCopulae. -- David. On Jul 2, 2010, at 7:06 PM, Ravi Varadhan wrote: There is no package called `integrate', but there is a function called `adaptIntegrate' in the cubature package. Ravi. -Original Message- From: r-help-boun...@r-project.org

[R] logistic regression - glm() - example in Dalgaard's book ISwR

2010-07-02 Thread Paulo Barata
Dear R-list members, I would like to pose a question about the use and results of the glm() function for logistic regression calculations. The question is based on an example provided on p. 229 in P. Dalgaard, Introductory Statistics with R, 2nd. edition, Springer, 2008. By means of this