[R] [R-pkgs] caret version 4.06 released

2009-01-26 Thread Max Kuhn
Version 4.06 of the caret package was sent to CRAN. caret can be used to tune the parameters of predictive models using resampling, estimate variable importance and visualize the results. There are also various modeling and helper functions that can be useful for training models. caret has

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

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

[R] How to analyse and model 2 time series, when one series needs to be differenced?

2009-01-26 Thread Andreas Klein
Hello. How can I analyse the cross-correlation between two time series with ccf, if one of the time series need to be differenced, so it is stationary? The two time series differ when in length and maybe ccf produces not the correct cross-correlation?! Another problem: How can I model the two

[R] HMISC package: wtd.table()

2009-01-26 Thread Norbert NEUWIRTH
Hi useRs developeRs, I got stuck within a function of the Hmisc package. Sounds easy, hope it is: I got 2 items (FamTyp.kurz, HGEW) of same length, no missings. length(FamTyp.kurz);summary(FamTyp.kurz) [1] 14883 Min. 1st Qu. MedianMean 3rd Qu.Max. 10.00 20.00 21.00 21.66

[R] heatmap with levelplot

2009-01-26 Thread Antje
Hi there, I'd like to create a heatmap from my matrix with a) a defined color range (lets say from yellow to red) b) using striking colors above and below a certain threshold (above = green, below = blue) Example matrix (there should be a few outliers generated...) + simple levelplot without

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

2009-01-26 Thread Norbert NEUWIRTH
oops, that really was easy: (wtd.table(FamTyp.kurz,HGEW,normwt=FALSE,na.rm=TRUE)) instead of (wtd.table(FamTyp.kurz,weigths=HGEW,normwt=FALSE,na.rm=TRUE)) sorry for that question ... Am 26.01.2009, 10:30 Uhr, schrieb Norbert NEUWIRTH norbert.s.neuwi...@univie.ac.at: Hi useRs developeRs,

[R] Help with clustering

2009-01-26 Thread mauede
I am going to try out a tentative clustering of some feature vectors. The range of values spanned by the three items making up the features vector is quite different: Item-1 goes roughly from 70 to 525 (integer numbers only) Item-2 is in-between 0 and 1 (all real numbers between 0 and 1) Item-3

Re: [R] infer haplotypes phasing trios tdthap

2009-01-26 Thread David Clayton
tdthap wassn't intended to solve that problem and it has been removed from my own web site since I no longer consider it important enough to support. DC -Original Message- From: Tiago R Magalhães [mailto:tiag...@gmail.com] Sent: 22 January 2009 11:10 To: r-help@R-project.org

[R] Meaning of Inner Product (%*%) Between Slot and Vector

2009-01-26 Thread Gundala Viswanath
Dear all, I have the following object and vector: print(alpha) Slot ra: [1] 0.994704478 0.002647761 0.000882587 0.000882587 0.000882587 0.989459074 [7] 0.005270463 0.002635231 0.002635231 0.994717023 0.005282977 1.0 [13] 1.0 Slot ja: [1] 1 5 2 3 4 2 1 3 4 3 3 4 5 Slot ia:

[R] how to modify an R built-in function?

2009-01-26 Thread diego Diego
Hello R experts! Last week I run in to a lot a problems triyng to fit an ARIMA model to a time series. The problem is that the internal process of the arima function call function optim to estimate the model parameters, so far so good... but my data presents a problem with the default method BFGS

Re: [R] Help with clustering

2009-01-26 Thread Christian Hennig
Generally, how to scale different variables when aggregating them in a dissimilarity measure is strongly dependent on the subject matter, what the aim of clustering and your cluster comncept is. This cannot be answered properly on such a mailing list. A standard transformation before

Re: [R] XML package help

2009-01-26 Thread Skewes,Aaron
Thanks! Works like a charm. -Aaron From: Duncan Temple Lang [dun...@wald.ucdavis.edu] Sent: Friday, January 23, 2009 6:48 PM To: Skewes,Aaron Cc: r-help@r-project.org Subject: Re: [R] XML package help Skewes,Aaron wrote: Please consider this: Manifest

[R] Mode (statistics) in R?

2009-01-26 Thread Jason Rupert
Hopefully this is a pretty simple question:   Is there a function in R that calculates the mode of a sample?   That is, I would like to be able to determine the value that occurs the most frequently in a data set.   I tried the default R mode function, but it appears to provide a storage type

[R] glm StepAIC with all interactions and update to remove a term vs. glm specifying all but a few terms and stepAIC

2009-01-26 Thread Robert Michael Inman
Problem: I am sorting through model selection process for first time and want to make sure that I have used glm, stepAIC, and update correctly. Something is strange because I get a different result between: 1) a glm of 12 predictor variables followed by a stepAIC where all interactions are

Re: [R] how to modify an R built-in function?

2009-01-26 Thread stephen sefick
type the name of the function in an R session. The source code will result- have fun. On Mon, Jan 26, 2009 at 7:36 AM, diego Diego dhab...@gmail.com wrote: Hello R experts! Last week I run in to a lot a problems triyng to fit an ARIMA model to a time series. The problem is that the internal

Re: [R] Mode (statistics) in R?

2009-01-26 Thread Carlos J. Gil Bellosta
Hello, You can try ?table. Best regards, Carlos J. Gil Bellosta http://www.datanaytics.com On Mon, 2009-01-26 at 05:28 -0800, Jason Rupert wrote: Hopefully this is a pretty simple question: Is there a function in R that calculates the mode of a sample? That is, I would like to be able

Re: [R] Mode (statistics) in R?

2009-01-26 Thread Mike Lawrence
Here's a rather convoluted way of finding the mode (or, at least, the first mode): x = round(rnorm(100,sd=5)) my_mode = as.numeric(names(table(x))[which.max(table(x))]) On Mon, Jan 26, 2009 at 9:28 AM, Jason Rupert jasonkrup...@yahoo.com wrote: Hopefully this is a pretty simple question:

Re: [R] Mode (statistics) in R?

2009-01-26 Thread Jason Rupert
Thanks.   I ended up breaking it up into two steps:   table_data-table(data) subset(table_data, table_data==max(table_data))   Thanks again. --- On Mon, 1/26/09, Mike Lawrence m...@thatmike.com wrote: From: Mike Lawrence m...@thatmike.com Subject: Re: [R] Mode (statistics) in R? To:

Re: [R] how to modify an R built-in function?

2009-01-26 Thread Gabor Grothendieck
Unless there is some real reason you need an arima model perhaps you could just try an ar model instead. ?ar On Mon, Jan 26, 2009 at 7:36 AM, diego Diego dhab...@gmail.com wrote: Hello R experts! Last week I run in to a lot a problems triyng to fit an ARIMA model to a time series. The

Re: [R] Mode (statistics) in R?

2009-01-26 Thread Marc Schwartz
on 01/26/2009 07:28 AM Jason Rupert wrote: Hopefully this is a pretty simple question: � Is there a function in R that calculates the mode of a sample?�� That is, I would like to be able to determine the value that occurs the most frequently in a data set. � I tried the default R mode

Re: [R] text vector clustering

2009-01-26 Thread San Miguel Martín , Eduardo
Dear srinivas, You can try using trigrams, a special case of N-grams, often used in Natural Language Processing. I am interested in grouping/cluster these names as those which are similar letter to letter. Are there any text clustering algorithm in R which can group names of similar type

Re: [R] heatmap with levelplot

2009-01-26 Thread Antje
I played a little bit around and got the following solution which works for now, though it seems to be too complicated to me. If anybody else know another solution - please let me know!!! library(lattice) my.mat - matrix(rnorm(800), nrow = 40) colorFun - colorRampPalette(c(yellow,red)) b -

Re: [R] Mode (statistics) in R?

2009-01-26 Thread patricia garcía gonzález
Hello, I think this will work: names( sort( -table( x ) ) )[1] Regards Patricia García From: c...@datanalytics.com To: jasonkrup...@yahoo.com Date: Mon, 26 Jan 2009 18:34:00 +0500 CC: r-help@r-project.org Subject: Re: [R] Mode (statistics) in R? Hello, You can try ?table. Best

Re: [R] Plotting graph for Missing values

2009-01-26 Thread bartjoosen
I added patientinformation1 variable and then I gave the command for tapply but its giving me the following error: Error in tapply(pat1, format(dos, %Y%m), function(x) sum(x == 0)) : arguments must have same length seems like you added patientinformation1, but still use pat1 in the

[R] Getting data from a PDF-file into R

2009-01-26 Thread joe1985
Hello I have around 200 PDF-documents, containing data i want organized in R as a dataframe. The PDF-documents look like this; http://www.nabble.com/file/p21667074/PRRS-billede%2Bmed%2Bfarver.jpeg or like this; http://www.nabble.com/file/p21667074/PRRS-billede%2Bmed%2Bfarver%2B2.jpeg So

Re: [R] RCurl unable to download a particular web page -- what is so special about this web page?

2009-01-26 Thread Tony Breyal
Hi, i ran your getURL example and had the same problem with downloading the file. ## R Start.. library(RCurl) toString(getURL(http://www.nytimes.com/2009/01/07/technology/business-computing/07program.html?_r=2;)) [1] ## R end. However, if it is interesting that if you manually save the page

Re: [R] Plotting graph for Missing values

2009-01-26 Thread Petr PIKAL
Hi Jim r-help-boun...@r-project.org napsal dne 26.01.2009 15:44:32: From your original posting: I tried the code which u provided. In place of dos in command pat1 - rbinom(length(dos), 1, .5) # generate some data I added patientinformation1 variable and then I gave the command for

[R] Help with sas.get

2009-01-26 Thread Sebastien Bihorel
Dear R-users, I am seeking advises on the sas.get function from the Hmisc package. I have tried to import some of our SAS files using the syntax presented in the function help example but the importation always failed. The function does not seem to recognize our sas files and complains about

Re: [R] Getting data from a PDF-file into R

2009-01-26 Thread Peter Dalgaard
joe1985 wrote: Hello I have around 200 PDF-documents, containing data i want organized in R as a dataframe. The PDF-documents look like this; http://www.nabble.com/file/p21667074/PRRS-billede%2Bmed%2Bfarver.jpeg or like this;

Re: [R] Getting data from a PDF-file into R

2009-01-26 Thread hadley wickham
On Mon, Jan 26, 2009 at 9:40 AM, Peter Dalgaard p.dalga...@biostat.ku.dk wrote: joe1985 wrote: Hello I have around 200 PDF-documents, containing data i want organized in R as a dataframe. The PDF-documents look like this;

[R] PCALG Package

2009-01-26 Thread Tibert, Brock
Hi all, Can anyone help me setup this package so I can use it. I am getting errors with the Rgraphviz package and have tried a number of ways to get this to work. Any help will be greatly appreciated! I am sorta new to R but have been actively trying to get into using it as my main

Re: [R] list.files changed in 2.7.0

2009-01-26 Thread davidr
Hmm. I get exactly the same files and directories with C: and C:/, except for the double slashes now. Previously the two calls to list.files gave exactly the same results. My current directory (getwd()) is not C:. I'm puzzled by your output. -- David -Original Message- From:

Re: [R] Getting data from a PDF-file into R

2009-01-26 Thread Henrique Dallazuanna
You can convert the pdf to text, then manipulate the output to read only the data. In linux has pdftotext function, in linux you can download the xpdf zip, that contais such function. Best On 1/26/09, joe1985 johan...@dsr.life.ku.dk wrote: Hello I have around 200 PDF-documents, containing

[R] RExcel foreground and background server

2009-01-26 Thread Irina Ursachi
Dear all, I have a question regarding background and foreground server in RExcel: Can somebody explain the main difference between them? As far as I understood from the RExcel webpage, for both of them one needs rights for access to Windows registries. The only difference that I can see so far

Re: [R] RCurl unable to download a particular web page -- what is so special about this web page?

2009-01-26 Thread Duncan Temple Lang
clair.crossup...@googlemail.com wrote: Dear R-help, There seems to be a web page I am unable to download using RCurl. I don't understand why it won't download: library(RCurl) my.url - http://www.nytimes.com/2009/01/07/technology/business-computing/07program.html?_r=2; getURL(my.url) [1]

[R] ANOVA with subsampling question

2009-01-26 Thread Wade Wall
Hi all, I am trying to analyze an experiment I ran, but not sure how to code in R. I have germinated seeds in petri dishes at 3 different temperatures (call it low, med, and high) and 2 different light levels (light and dark). For each seed I have recorded time to germination (not counting

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

2009-01-26 Thread Frank E Harrell Jr
Norbert NEUWIRTH wrote: oops, that really was easy: (wtd.table(FamTyp.kurz,HGEW,normwt=FALSE,na.rm=TRUE)) instead of (wtd.table(FamTyp.kurz,weigths=HGEW,normwt=FALSE,na.rm=TRUE)) That is one solution. The other is to spell 'weights' correctly :-) Frank sorry for that question ... Am

Re: [R] commercially supported version of R for 64 -bit Windows?

2009-01-26 Thread David M Smith
That's correct: REvolution Computing (whom I work for) is in the process of porting R and packages to 64-bit Windows. The development process has been underway for several months and is near completion. There will be a beta test in February and we expect to release in March. When the beta

[R] reshape problem: id and variable names not being recognized

2009-01-26 Thread MW Frost
Hi everyone. Long time listener, first-time caller here. I have a data set that's been melted with the excellent reshape package, but I can't seem to cast it the way I need to. Here's the melted data's structure: str(mdat) 'data.frame': 6978 obs. of 4 variables: $ VehType : Factor w/ 2

Re: [R] RCurl unable to download a particular web page -- what is so special about this web page?

2009-01-26 Thread Jeffrey Horner
Duncan Temple Lang wrote: clair.crossup...@googlemail.com wrote: Dear R-help, There seems to be a web page I am unable to download using RCurl. I don't understand why it won't download: library(RCurl) my.url -

Re: [R] reshape problem: id and variable names not being recognized

2009-01-26 Thread jim holtman
Look at your 'str(mdat)' and you will see that there is not a variable call 'mpg'; it is one of the levels of the 'variable'. On Mon, Jan 26, 2009 at 11:38 AM, MW Frost mwfr...@gmail.com wrote: Hi everyone. Long time listener, first-time caller here. I have a data set that's been melted with

[R] Tinn-R

2009-01-26 Thread Tibert, Brock
Hi Everyone, I was hoping someone could help me with the settings for Tinn-R. I see in the screen shots that it has syntax help, or something similar (tips on what functions, etc). I can not seem to get this to turn on in the program, and I am wondering if I have to set up a few options. I

[R] Goodness of fit for gamma distributions

2009-01-26 Thread Dan31415
I'm looking for goodness of fit tests for gamma distributions with large data sizes. I have a matrix with around 10,000 data values in it and i have fitted a gamma distribution over a histogram of the data. The problem is testing how well that distribution fits. Chi-squared seems to be used

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

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

[R] error managment

2009-01-26 Thread diego Diego
Hello R experts! I'm running a FOR loop in which at every step an arima model is generated. The problem is some series produces numeric problems with optim. My question is if there is a way of telling to R that at every critical error of optim jumps to the next series instead of stopping the

Re: [R] error managment

2009-01-26 Thread Gabor Grothendieck
?try On Mon, Jan 26, 2009 at 1:05 PM, diego Diego dhab...@gmail.com wrote: Hello R experts! I'm running a FOR loop in which at every step an arima model is generated. The problem is some series produces numeric problems with optim. My question is if there is a way of telling to R that at

[R] problem with table()

2009-01-26 Thread Dominik Hattrup
Hey everyone, I am looking for the easiest way to get this table # Table2 Year / 2000 / 2002 / 2004 Julia / 3 / 4 / 1 Peter / 1 / 2 / 4 ... / ... / ... / ... out of this one? # Table1 name / year / cases Julia / 2000 / 1 Julia / 2000 / 2 Julia / 2002

Re: [R] problem with table()

2009-01-26 Thread Marc Schwartz
on 01/26/2009 12:23 PM Dominik Hattrup wrote: Hey everyone, I am looking for the easiest way to get this table # Table2 Year / 2000 / 2002 / 2004 Julia / 3 / 4 / 1 Peter / 1 / 2 / 4 ... / ... / ... / ... out of this one? # Table1 name / year / cases Julia

[R] Large regular expressions

2009-01-26 Thread Stavros Macrakis
Given a vector of reference strings Ref and a vector of test strings Test, I would like to find elements of Test which do not contain elements of Ref as \b-delimited substrings. This can be done straightforwardly for length(Ref) 6000 or so (R 2.8.1 Windows) by constructing a pattern like

[R] name scoping within dataframe index

2009-01-26 Thread Alexy Khrabrov
Every time I have to prefix a dataframe column inside the indexing brackets with the dataframe name, e.g. df[df$colname==value,] -- I am wondering, why isn't there an R scoping rule that search starts with the dataframe names, as if we'd said with(df, df[colname==value,]) -- wouldn't

Re: [R] problem with table()

2009-01-26 Thread Dominik Hattrup
Marc Schwartz schrieb: on 01/26/2009 12:23 PM Dominik Hattrup wrote: Hey everyone, I am looking for the easiest way to get this table # Table2 Year / 2000 / 2002 / 2004 Julia / 3 / 4 / 1 Peter / 1 / 2 / 4 ... / ... / ... / ... out of this one? # Table1 name /

Re: [R] name scoping within dataframe index

2009-01-26 Thread Duncan Murdoch
On 1/26/2009 1:46 PM, Alexy Khrabrov wrote: Every time I have to prefix a dataframe column inside the indexing brackets with the dataframe name, e.g. df[df$colname==value,] -- I am wondering, why isn't there an R scoping rule that search starts with the dataframe names, as if we'd said

Re: [R] Large regular expressions

2009-01-26 Thread Gabor Grothendieck
I am using R.version.string # Vista [1] R version 2.8.1 Patched (2008-12-26 r47350) and it also caused R to actually crash for me. On Mon, Jan 26, 2009 at 1:38 PM, Stavros Macrakis macra...@alum.mit.edu wrote: Given a vector of reference strings Ref and a vector of test strings Test, I would

Re: [R] name scoping within dataframe index

2009-01-26 Thread Gabor Grothendieck
Try: subset(df, colname == value) On Mon, Jan 26, 2009 at 1:46 PM, Alexy Khrabrov delivera...@gmail.com wrote: Every time I have to prefix a dataframe column inside the indexing brackets with the dataframe name, e.g. df[df$colname==value,] -- I am wondering, why isn't there an R scoping

Re: [R] Problem with colormodel in pdf driver

2009-01-26 Thread Luis Torgo
Greg Snow wrote: You may want to consider a dotchart instead of a barplot. Then you can distinguish between groups by using symbols, grouping, and labels rather than depending on colors/shades of grey. Thanks Greg. The only problem is that I was trying to illustrate the use of barplot()

Re: [R] name scoping within dataframe index

2009-01-26 Thread Alexy Khrabrov
On 1/26/2009 1:46 PM, Alexy Khrabrov wrote: Every time I have to prefix a dataframe column inside the indexing brackets with the dataframe name, e.g. df[df$colname==value,] -- I am wondering, why isn't there an R scoping rule that search starts with the dataframe names, as if we'd said

Re: [R] name scoping within dataframe index

2009-01-26 Thread Duncan Murdoch
On 1/26/2009 2:01 PM, Alexy Khrabrov wrote: On 1/26/2009 1:46 PM, Alexy Khrabrov wrote: Every time I have to prefix a dataframe column inside the indexing brackets with the dataframe name, e.g. df[df$colname==value,] -- I am wondering, why isn't there an R scoping rule that search starts

Re: [R] Tinn-R

2009-01-26 Thread anna freni sterrantino
Hi, if you are on linux Emacs + ESS is quite popular too. Cheers Anna Anna Freni Sterrantino Ph.D Student Department of Statistics University of Bologna, Italy via Belle Arti 41, 40124 BO. Da: Tibert, Brock btib...@bentley.edu A: r-help@r-project.org

Re: [R] Tinn-R

2009-01-26 Thread Thomas Roth (geb. Kaliwe)
First thing you need to do is, save the file with an .r ending. Tinn-R will then come up with syntax highlighting (as far as i remember). Options-Main-Application leads you to the r-configuration... HTH Thomas Tibert, Brock schrieb: Tibert, Brock schrieb: Hi Everyone, I was hoping

Re: [R] name scoping within dataframe index

2009-01-26 Thread Alexy Khrabrov
On Jan 26, 2009, at 2:12 PM, Duncan Murdoch wrote: df[get(colname, parent.frame()) == value,] Actually, what I propose is a special search rule which simply looks at the enclosing dataframe.name[...] outside the brackets and looks up the columns first. Yes, I understood that, and I

[R] Spectral analysis with mtm-svd Multi-Taper Method Combined with Singular Value Decomposition

2009-01-26 Thread Yasir Kaheil
Hi list, Does anyone know if there is a library in R that does MTM-SVD method for spectral analysis? Thanks - Yasir H. Kaheil Columbia University -- View this message in context:

[R] randomSurvivalForest plotting

2009-01-26 Thread A Van Dyke
i would like to plot a subset of variables with the highest variable importance measures (say the top 20) instead of plotting all of the variables included in the analysis (~75). i tried arguments that work to restrict the number of variables displayed in the plot in randomForest as follows:

Re: [R] name scoping within dataframe index

2009-01-26 Thread Duncan Murdoch
On 1/26/2009 2:20 PM, Alexy Khrabrov wrote: On Jan 26, 2009, at 2:12 PM, Duncan Murdoch wrote: df[get(colname, parent.frame()) == value,] Actually, what I propose is a special search rule which simply looks at the enclosing dataframe.name[...] outside the brackets and looks up the

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

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

Re: [R] Problem with colormodel in pdf driver

2009-01-26 Thread Prof Brian Ripley
On Mon, 26 Jan 2009, Luis Torgo wrote: Greg Snow wrote: You may want to consider a dotchart instead of a barplot. Then you can distinguish between groups by using symbols, grouping, and labels rather than depending on colors/shades of grey. Thanks Greg. The only problem is that I was

Re: [R] [stats-rosuda-devel] Problem with JGR. Was: Re: Using help()

2009-01-26 Thread Simon Urbanek
On Jan 25, 2009, at 9:35 , Michael Kubovy wrote: Dear Friends, Thanks to Rolf Turner, Brian Ripley and Patrick Burns for their answers.They don't quite resolve the problem, which I now realize is due to non-standard behavior of JGR, at least on my machine (I verified that Mac GUI works

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

2009-01-26 Thread Peter Dalgaard
Frank E Harrell Jr wrote: I'm still learning english myself. Including capitalization rules? -- O__ Peter Dalgaard Øster Farimagsgade 5, Entr.B c/ /'_ --- Dept. of Biostatistics PO Box 2099, 1014 Cph. K (*) \(*) -- University of Copenhagen Denmark Ph:

[R] Power analysis for MANOVA?

2009-01-26 Thread Adam D. I. Kramer
Hello, I have searched and failed for a program or script or method to conduct a power analysis for a MANOVA. My interest is a fairly simple case of 5 dependent variables and a single two-level categorical predictor (though the categories aren't balanced). If anybody happens to

Re: [R] Power analysis for MANOVA?

2009-01-26 Thread Mitchell Maltenfort
http://www.amazon.com/Statistical-Power-Analysis-Behavioral-Sciences/dp/0805802835 Cohen's book was in fact the basis for the pwr package at CRAN. And it does have a MANOVA power analysis, which was left out of the pwr package. On Mon, Jan 26, 2009 at 4:12 PM, Adam D. I. Kramer

Re: [R] Power analysis for MANOVA?

2009-01-26 Thread Stephan Kolassa
Hi Adam, My (and, judging from previous traffic on R-help about power analyses, also some other people's) preferred approach is to simply simulate an effect size you would like to detect a couple of thousand times, run your proposed analysis and look how often you get significance. In your

[R] Sweave'ing Danish characters

2009-01-26 Thread Peter Jepsen
Hi, I am writing an Sweave document and am using 'xtable' to make frequency tables of diagnoses of people undergoing cholecystectomy. Some of these diagnoses contain Danish characters (æ, ø, and å), and these characters are all garbled in the Latex document after I run Sweave. The odd thing

[R] suppressing time shift in plot of POSIXct object?

2009-01-26 Thread Jim Porzak
Friends, I have a POSIXct vector located in EST timezone. When I plot against it here in PST, the time axis is shifted 3 hours back in time. IOW, plot adjusts for time zone difference. Now that's really great, if that's what one wants. However, I want time axis to use actual times in object

Re: [R] Power analysis for MANOVA?

2009-01-26 Thread Adam D. I. Kramer
On Mon, 26 Jan 2009, Stephan Kolassa wrote: My (and, judging from previous traffic on R-help about power analyses, also some other people's) preferred approach is to simply simulate an effect size you would like to detect a couple of thousand times, run your proposed analysis and look how

Re: [R] Problem with colormodel in pdf driver

2009-01-26 Thread Luis Torgo
Prof Brian Ripley wrote: On Mon, 26 Jan 2009, Luis Torgo wrote: Greg Snow wrote: You may want to consider a dotchart instead of a barplot. Then you can distinguish between groups by using symbols, grouping, and labels rather than depending on colors/shades of grey. Thanks Greg. The only

Re: [R] Power analysis for MANOVA?

2009-01-26 Thread Charles C. Berry
If you know what a 'general linear hypothesis test' is see http://cran.r-project.org/src/contrib/Archive/hpower/hpower_0.1-0.tar.gz HTH, Chuck On Mon, 26 Jan 2009, Adam D. I. Kramer wrote: On Mon, 26 Jan 2009, Stephan Kolassa wrote: My (and, judging from previous traffic on

[R] Error in segmented() output from segmented package

2009-01-26 Thread tsippel
Hi- I'm gettting the following error message when trying to use the segmented function to look for breakpoints in my data. Error in segmented.glm(glm, seg.Z = ~segmentdist, psi = 2, control = seg.control(display = F), : (Some) estimated psi out of its range Here are some real data and the

[R] why two diff. se in nlsList?

2009-01-26 Thread Tao Shi
Hi list, In the object returned by summary.nlsList, what's the difference between coefficients and parameters? The have the same Estimate, different se (therefore t value), but same p values. R.2.8.0 on winxp with nlme_3.1-89 Thanks, ...Tao

[R] Error in Surv(time, status) : Time variable is not numeric

2009-01-26 Thread Braem M
Dear, I want to analyze two-level survival data using a shared frailty model, for which I want to use the R package 'Frailtypack, proposed by Rondeau et al. The dataset was built using SAS software. I also tried to change the format using SPSS and Excell. My (reduced) dataset has following

Re: [R] Power analysis for MANOVA?

2009-01-26 Thread Adam D. I. Kramer
On Mon, 26 Jan 2009, Charles C. Berry wrote: If you know what a 'general linear hypothesis test' is see http://cran.r-project.org/src/contrib/Archive/hpower/hpower_0.1-0.tar.gz I do, and am quite interested, however this package will not install on R 2.8.1: First, it said that

Re: [R] suppressing time shift in plot of POSIXct object?

2009-01-26 Thread jim holtman
Try: Sys.setenv(TZ=EST) plot(y ~ t, type = l) You can save TZ before you set it and then restore it. On Mon, Jan 26, 2009 at 5:47 PM, Jim Porzak jpor...@gmail.com wrote: Friends, I have a POSIXct vector located in EST timezone. When I plot against it here in PST, the time axis is shifted 3

[R] working with tables -- was Re: Mode (statistics) in R?

2009-01-26 Thread Carl Witthoft
Ok, so I'm slowly figuring out what a factor is, and was able to follow the related thread about finding a mode by using constructs like my_mode = as.numeric(names(table(x))[which.max(table(x))]) Now, suppose I want to keep looking for other modes? For example, Rgames

Re: [R] The Quality Accuracy of R

2009-01-26 Thread Gabor Grothendieck
It would be possible to develop tools to develop code coverage statistics quantifying the percent of the code that the tests exercise. On Fri, Jan 23, 2009 at 10:04 AM, Muenchen, Robert A (Bob) muenc...@utk.edu wrote: Hi All, We have all had to face skeptical colleagues asking if software

Re: [R] The Quality Accuracy of R

2009-01-26 Thread Muenchen, Robert A (Bob)
That's a great idea. I know of no commercial vendors who provide such detailed info. Bob -Original Message- From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com] Sent: Monday, January 26, 2009 7:52 PM To: Muenchen, Robert A (Bob) Cc: R-help@r-project.org Subject: Re: [R] The Quality

[R] Pausing processing into an interactive session

2009-01-26 Thread Zhou Fang
Hi all, As a possibly silly request, is it possible to interactively pause a R-calculation and do a browser(), say, without browser or other debug handlers being explicitly included in the code? Imagine the following situation: You write up a big calculation for R to calculate. We are talking

[R] WhisperStation R

2009-01-26 Thread zerfetzen
What do you think of this: http://www.microway.com/whisperstation/whisperstation-r.html I'm considering ditching my Windows Vista 2 GB RAM computer for WhisperStation R using Debian 64-bit Linux with 32 GB RAM and setting the whole thing up for R and WinBUGS. I put in a price request, but I

Re: [R] Running R under Sun Grid Engine with OpenMPI tight integration

2009-01-26 Thread Peter Waltman
Hi - I saw your posting on the R-help mailing list. Were you ever able to get this working? did you end up switching to use the rsge library? I'm trying to do the same, and not having very much luck getting it going. Thanks! Peter Waltman [[alternative HTML version deleted]]

[R] Problem with loading RMySQL under sge/qsub

2009-01-26 Thread Peter Waltman
Hi - I'm trying to set up a parallelized batch job that is run under rmpi and managed by sge, using qsub, but it reports that it can't load RMySQL because it can't find the libmysqlclient.so.15 file. Note, when I run R interactively, and manually load the RMySQL library, it works without a

Re: [R] Power analysis for MANOVA?

2009-01-26 Thread Charles C. Berry
On Mon, 26 Jan 2009, Adam D. I. Kramer wrote: On Mon, 26 Jan 2009, Charles C. Berry wrote: If you know what a 'general linear hypothesis test' is see http://cran.r-project.org/src/contrib/Archive/hpower/hpower_0.1-0.tar.gz I do, and am quite interested, however this package will not

Re: [R] WhisperStation R

2009-01-26 Thread andrew
any idea why DDR2 667 MHz RAM isn't used instead of DDR? I thought that DDR 400MHz was almost finished in production... On Jan 27, 1:01 pm, zerfetzen zerfet...@yahoo.com wrote: What do you think of this: http://www.microway.com/whisperstation/whisperstation-r.html I'm considering ditching

Re: [R] Running R under Sun Grid Engine with OpenMPI tight integration

2009-01-26 Thread Rainer M Krug
On Tue, Jan 27, 2009 at 2:30 AM, Peter Waltman peter.walt...@gmail.com wrote: Hi - I saw your posting on the R-help mailing list. Were you ever able to get this working? did you end up switching to use the rsge library? Yes - that is exactly what I did - I am using rsge or, which is in most

[R] Control of Quartz Window Location

2009-01-26 Thread Larry Weldon
If I use plot(1:10) quartz() plot(1:10) I get the second graph window almost on top of the first graph window. How can I control the location of the quartz window? Larry Weldon Simon Fraser University wel...@sfu.ca www.stat.sfu.ca/~weldon __

[R] How do you specify font family in png output; png cross-platform issues

2009-01-26 Thread Paul Johnson
For teaching purposes, I prepared a little R program. I want to give this to students who can run it and dump out many formats and then compare their use in LaTeX documents. I do not have too much trouble with xfig or postscript format, but I've really run into a roadblock where png files are

Re: [R] Getting data from a PDF-file into R

2009-01-26 Thread joe1985
Peter Dalgaard wrote: joe1985 wrote: Hello I have around 200 PDF-documents, containing data i want organized in R as a dataframe. The PDF-documents look like this; http://www.nabble.com/file/p21667074/PRRS-billede%2Bmed%2Bfarver.jpeg or like this;