Re: [R] newton method

2009-03-24 Thread Yihui Xie
I'm not sure what you meant by a topic on newton's method (algorithm? demo?), but the demonstration in the package 'animation' might help: install.packages('animation') par(pch = 20) ani.options(nmax = 50) newton.method(function(x) 5 * x^3 - 7 * x^2 - 40 * x + 100, 7.15, c(-6.2, 7.1))

Re: [R] Legend containing maths symbol and values of variables

2009-03-24 Thread Daren Tan
I tried sprintf and paste, which solves one of my problem i.e. replacing the variable name by its value. However, = is not renderend correctly as the inequality symbol. On Tue, Mar 24, 2009 at 12:24 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: For the legend argument try this: leg =

[R] Calculating percentage Missing value for variables using one object

2009-03-24 Thread Shreyasee
Hi, I have a dataset in which there are in all 250 variables and for each variable the data is entered over the months. I need to calculate the percentage of missing values for each variable over each month and then plot a graph for that. I am running the following code for doing the same *ds -

[R] Barplot With Selected X-Axis Region

2009-03-24 Thread Gundala Viswanath
Dear all, I have a data that looks like this: print(dat) V1 V2 1 -43342073 14 2 -433377304 3 -432846761 11372 75188572 11 11373 752061656 11374 75262720 24 What I want to do is to have a barplot where x-axis is taken from V1 and y-axis

Re: [R] If statement generates two outputs

2009-03-24 Thread Wacek Kusnierczyk
Berwin A Turlach wrote: G'day Carl, On Mon, 23 Mar 2009 20:11:19 -0400 Carl Witthoft c...@witthoft.com wrote: From: Wacek Kusnierczyk Waclaw.Marcin.Kusnierczyk_at_idi.ntnu.no Date: Sun, 22 Mar 2009 22:58:49 +0100 just for fun, you could do this with multiassignment, e.g., using

Re: [R] How to set up a function for Central Limit Theorem

2009-03-24 Thread Thomas Lumley
On Mon, 23 Mar 2009, pfc_ivan wrote: I tried using the for (i..) to make 1000 differents sets of numbers, but then I dont know how to get the mean value of all of them... because I dont even think 1000 different sets of numbers were made, because when i print it it always shows me the same

Re: [R] read in large data file (tsv) with inline filter?

2009-03-24 Thread Thomas Lumley
On Mon, 23 Mar 2009, David Reiss wrote: I have a very large tab-delimited file, too big to store in memory via readLines() or read.delim(). Turns out I only need a few hundred of those lines to be read in. If it were not so large, I could read the entire file in and grep the lines I need. For

Re: [R] confidence interval or error of x intercept of a linear regression

2009-03-24 Thread Thomas Lumley
On Mon, 23 Mar 2009, Kevin J Emerson wrote: Now, it is simple enough to calculate the x-intercept itself ( - intercept / slope ), but it is a whole separate process to generate the confidence interval of it. I can't figure out how to propagate the error of the slope and intercept into

Re: [R] confidence interval or error of x intercept of a linear regression

2009-03-24 Thread Thomas Lumley
I should probably also add the health warning that the delta-method approach to ratios of coefficients is REALLY BAD unless the denominator is well away from zero. When the denominator is near zero nothing really works (Fieller's method is valid but only because 'valid' means less than you

Re: [R] If statement generates two outputs

2009-03-24 Thread Wacek Kusnierczyk
Berwin A Turlach wrote: G'day Carl, On Mon, 23 Mar 2009 20:11:19 -0400 Carl Witthoft c...@witthoft.com wrote: From: Wacek Kusnierczyk Waclaw.Marcin.Kusnierczyk_at_idi.ntnu.no Date: Sun, 22 Mar 2009 22:58:49 +0100 just for fun, you could do this with multiassignment, e.g., using

[R] How to separate huge dataset into chunks

2009-03-24 Thread Guillaume Filteau
Hello all, I’m trying to take a huge dataset (1.5 GB) and separate it into smaller chunks with R. So far I had nothing but problems. I cannot load the whole dataset in R due to memory problems. So, I instead try to load a few (10) lines at a time (with read.table). However, R kept

[R] Which features would you like to see on the crantastic.org community portal?

2009-03-24 Thread Bjørn Arild Mæland
Hi, As some of you know, one of the R Projects ideas for GSOC 2009 is to expand the crantastic.org portal. Some ideas for what should be done is mentioned in Hadley Wickham's project description [1], others have been briefly mentioned on this mailing list before. Undoubtedly, however, there are

Re: [R] If statement generates two outputs

2009-03-24 Thread Wacek Kusnierczyk
Berwin A Turlach wrote: G'day Carl, On Mon, 23 Mar 2009 20:11:19 -0400 Carl Witthoft c...@witthoft.com wrote: But seriously: can someone explain to me what's going on in the rvalues.r code? I tried a simple experiment, replacing := with a colec in the code, and of course the line

Re: [R] If statement generates two outputs

2009-03-24 Thread Romain Francois
Wacek Kusnierczyk wrote: Berwin A Turlach wrote: G'day Carl, On Mon, 23 Mar 2009 20:11:19 -0400 Carl Witthoft c...@witthoft.com wrote: But seriously: can someone explain to me what's going on in the rvalues.r code? I tried a simple experiment, replacing := with a colec in the

Re: [R] If statement generates two outputs

2009-03-24 Thread Wacek Kusnierczyk
Romain Francois wrote: This is probably due to that in the gram.y file : case ':': if (nextchar(':')) { if (nextchar(':')) { yylval = install(:::); return NS_GET_INT; } else { yylval = install(::); return NS_GET; } } if

Re: [R] If statement generates two outputs

2009-03-24 Thread Thomas Lumley
On Tue, 24 Mar 2009, Wacek Kusnierczyk wrote: thanks. so it seems to be intentionally parsable, though i wouldn't say that this gives a meaning to ':=' -- the operator has a syntactic category, but no semantics. the syntactic category does not imply any semantics, as in '-' = function(a,

Re: [R] How to separate huge dataset into chunks

2009-03-24 Thread Thomas Lumley
On Tue, 24 Mar 2009, Guillaume Filteau wrote: Hello all, I’m trying to take a huge dataset (1.5 GB) and separate it into smaller chunks with R. So far I had nothing but problems. I cannot load the whole dataset in R due to memory problems. So, I instead try to load a few (10) lines at

Re: [R] If statement generates two outputs

2009-03-24 Thread Wacek Kusnierczyk
Thomas Lumley wrote: On Tue, 24 Mar 2009, Wacek Kusnierczyk wrote: thanks. so it seems to be intentionally parsable, though i wouldn't say that this gives a meaning to ':=' -- the operator has a syntactic category, but no semantics. the syntactic category does not imply any semantics, as

Re: [R] If statement generates two outputs

2009-03-24 Thread Thomas Lumley
On Tue, 24 Mar 2009, Wacek Kusnierczyk wrote: thanks for the explanation. for the sake of rvalues, please do not fix this bug. i find it useful to avoid messing with = and -. If you insist on calling it a bug, the chances of its getting fixed are higher. -thomas Thomas Lumley

[R] Random sampling based on the observations

2009-03-24 Thread Arup
Hello!I am having a problem with Random sampling in R. I have used a syntax: mydata.sub=sample(mydata,7,replace=FALSE,prob=NULL) which allows me to choose a random sample based on the variables(correct me if I am wrong!). Suppose I have 10 variable and if I use the above mentioned command then it

[R] Unit root

2009-03-24 Thread Brajkovic J.
I am confused by obtaining different results when testing for unit root when using different packages. I have 2625 price entries for which I want to determine whether they exhibit unit root. First I test using adf.test from tseries package by running: adf.test(P, k=30) Augmented Dickey-Fuller

[R] cross-validation

2009-03-24 Thread per243
I have reviewed all the scripts that appear http://cran.es.r-project.org/ and I cann´t find any suitable for cross-validation with a model of the form y = aX^(b). exp(cZ). Please can someone help me? Thanks, a lot of!! -- View this message in context:

Re: [R] help: what are the basis functions in {mgcv}: gam?

2009-03-24 Thread oliviax
Thanks for Andrew Robinson's suggestion. I am sorry that I have been looking for the book Wood (2006) Generalized Additive Models: Introduction with R for a long time in my univerity and book stores, but I failed to find one. I think it may take a long time to buy one from the Internet, and we

[R] Change teh datatype from factor to double

2009-03-24 Thread Maxl18
Hi all, I want to read a .txt dataset, but when I test the type, some of the data is a factor. How can I avoid them to get a factor or how can I unfactor them without loosing their values? Thanks, Max -- View this message in context:

Re: [R] Random sampling based on the observations

2009-03-24 Thread Adrian Dusa
On Tuesday 24 March 2009, you wrote: Hello!I am having a problem with Random sampling in R. I have used a syntax: mydata.sub=sample(mydata,7,replace=FALSE,prob=NULL) which allows me to choose a random sample based on the variables(correct me if I am wrong!). Actually, no. You are sampling

Re: [R] confidence interval or error of x intercept of a linear

2009-03-24 Thread Ted Harding
On 24-Mar-09 03:31:32, Kevin J Emerson wrote: Hello all, This is something that I am sure has a really suave solution in R, but I can't quite figure out the best (or even a basic) way to do it. I have a simple linear regression that is fit with lm for which I would like to estimate the x

Re: [R] If statement generates two outputs

2009-03-24 Thread Wacek Kusnierczyk
Thomas Lumley wrote: On Tue, 24 Mar 2009, Wacek Kusnierczyk wrote: thanks for the explanation. for the sake of rvalues, please do not fix this bug. i find it useful to avoid messing with = and -. If you insist on calling it a bug, the chances of its getting fixed are higher. no no,

Re: [R] help: what are the basis functions in {mgcv}: gam?

2009-03-24 Thread Simon Wood
If you can't get hold of the book, try appendix A of http://www.maths.bath.ac.uk/~sw283/simon/papers/psm.pdf -- what's described there is the cr basis best, Simon On Tuesday 24 March 2009 05:01, oliviax wrote: Thanks for Andrew Robinson's suggestion. I am sorry that I have been looking for

[R] install.package(TinnR) - there is no package called 'TinnR'

2009-03-24 Thread Martin Tomko
I have troubles make TinnR 2.2.0.2 work, it seems that the dependency on the package TinnR that cannot be found (I tried also manual downloads, but I cannot find the package anywhere on any CRAN mirror). I even set a default cran mirror in the Rprofile.site file, so that the later command can

[R] r online

2009-03-24 Thread Thomas Steiner
Hi, I'd like to execute simple commands and functions in R through a website, is there any service like this somewhere? I only found http://www.osvisions.com/r-online/ but it does not work (last update 2003) and the links to releated websites only give errors (if I calculate 7+3). Thanks for help

Re: [R] Legend containing maths symbol and values of variables

2009-03-24 Thread Gabor Grothendieck
Try this: L - list(bquote(Up ( = .(threshold) ~ )), Normal, Down, NA) legend(top, leg = L) On Tue, Mar 24, 2009 at 2:52 AM, Daren Tan darenta...@gmail.com wrote: I tried sprintf and paste, which solves one of my problem i.e. replacing the variable name by its value. However, = is not

Re: [R] Legend containing maths symbol and values of variables

2009-03-24 Thread Uwe Ligges
So what about threshold - 0.5 plot(NA, xlab=, ylab=, main=, axes=F, xlim=c(0,1), ylim=c(0,1), xaxs=i, yaxs=i) legend1 - substitute(Up ( * phantom() = threshold * ), list(threshold=threshold)) legend(x=0, y=1, fill=c(orange, white, green, gray90), cex=0.8, pt.cex=1.8, x.intersp=0.6, bg

[R] How to find the path or the current file?

2009-03-24 Thread Marie Sivertsen
Dear useRs, I have a collection of source file and some of these call others. The files are distribute among a number of directories, and to know how to call some other file they need to know what file is currently executed. As example, I have a file 'search.R' located in directory 'bin' which

[R] Variogram with Gstat

2009-03-24 Thread giuseppef...@libero.it
Dear all can you help me? i have this problem: i have a dataset in a text file in a matrix of 3 columns: x, y, z where x and y are the coordinates and Z are the mesurements. How can i do a variogram with R? __ R-help@r-project.org mailing list

Re: [R] read in large data file (tsv) with inline filter?

2009-03-24 Thread Gabor Grothendieck
On Mon, Mar 23, 2009 at 5:53 PM, David Reiss dre...@systemsbiology.org wrote: I have a very large tab-delimited file, too big to store in memory via readLines() or read.delim(). Turns out I only need a few hundred of those lines to be read in. If it were not so large, I could read the entire

Re: [R] How to find the path or the current file?

2009-03-24 Thread Duncan Murdoch
On 24/03/2009 7:16 AM, Marie Sivertsen wrote: Dear useRs, I have a collection of source file and some of these call others. The files are distribute among a number of directories, and to know how to call some other file they need to know what file is currently executed. As example, I have a

Re: [R] Legend containing maths symbol and values of variables

2009-03-24 Thread Gabor Grothendieck
that should have been: L - list(bquote(Up ( = .(threshold) ~ )), Normal, Down, NA) legend(top, leg = as.expression(L)) On Tue, Mar 24, 2009 at 7:18 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: Try this: L - list(bquote(Up ( = .(threshold) ~ )), Normal, Down, NA) legend(top, leg =

Re: [R] r online

2009-03-24 Thread Gerhard Schön
try: http://bayes.math.montana.edu/Rweb/Rweb.general.html Gerhard Thomas Steiner schrieb am 24.03.2009 12:06: Hi, I'd like to execute simple commands and functions in R through a website, is there any service like this somewhere? I only found http://www.osvisions.com/r-online/ but it does not

Re: [R] Change teh datatype from factor to double

2009-03-24 Thread David Winsemius
?read.table Consult the help file for the parameters: as.is =or stringsAsFactors = -- David Winsemius On Mar 24, 2009, at 4:22 AM, Maxl18 wrote: Hi all, I want to read a .txt dataset, but when I test the type, some of the data is a factor. How can I avoid them to get a factor or how

Re: [R] r online

2009-03-24 Thread Dirk Eddelbuettel
On 24 March 2009 at 12:06, Thomas Steiner wrote: | Hi, | I'd like to execute simple commands and functions in R through a | website, is there any service like this somewhere? | I only found http://www.osvisions.com/r-online/ but it does not work | (last update 2003) and the links to releated

Re: [R] How to find the path or the current file?

2009-03-24 Thread Gabor Grothendieck
See: https://stat.ethz.ch/pipermail/r-help/2009-January/184745.html On Tue, Mar 24, 2009 at 7:16 AM, Marie Sivertsen mariesiv...@gmail.com wrote: Dear useRs, I have a collection of source file and some of these call others.  The files are distribute among a number of directories, and to know

Re: [R] Legend containing maths symbol and values of variables

2009-03-24 Thread Daren Tan
Thanks, that works beautifully. Without your help, I will be struggling for a long time. On Tue, Mar 24, 2009 at 7:46 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: that should have been: L - list(bquote(Up ( = .(threshold) ~ )), Normal, Down, NA) legend(top, leg = as.expression(L))

Re: [R] read in large data file (tsv) with inline filter?

2009-03-24 Thread Dirk Eddelbuettel
Hi David, On 23 March 2009 at 15:09, Dylan Beaudette wrote: | On Monday 23 March 2009, David Reiss wrote: | I have a very large tab-delimited file, too big to store in memory via | readLines() or read.delim(). Turns out I only need a few hundred of those | lines to be read in. If it were not

Re: [R] mvpart error

2009-03-24 Thread Gavin Simpson
On Sun, 2009-03-22 at 22:45 -0400, Josh Stumpf wrote: Hello all, When attempting a classification tree using mvpart, I get the following error: thesis2.mvp=mvpart(bat_sp~., data=alltrees.df) Error in all(keep) : unused argument(s) (c(TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE,

Re: [R] install.package(TinnR) - there is no package called 'TinnR'

2009-03-24 Thread David Winsemius
Tinn-R is not an R package. It is a standalone text editor: http://www.lmgtfy.com/?q=tinn-r -- David Winsemius On Mar 24, 2009, at 7:05 AM, Martin Tomko wrote: I have troubles make TinnR 2.2.0.2 work, it seems that the dependency on the package TinnR that cannot be found (I tried also

Re: [R] install.package(TinnR) - there is no package called 'TinnR'

2009-03-24 Thread Martin Tomko
David, that is not helpful. I KNOW that TinnR is a standalone editor. If you had a look at the Rprofile.site required by TinnR, you would notice the part of the code I send earlier: # check necesary packages necessary = c('TinnR', 'svSocket') installed = necessary %in% installed.packages()[,

Re: [R] r online

2009-03-24 Thread David Winsemius
Google with rweb -- David Winsemius On Mar 24, 2009, at 7:06 AM, Thomas Steiner wrote: Hi, I'd like to execute simple commands and functions in R through a website, is there any service like this somewhere? I only found http://www.osvisions.com/r-online/ but it does not work (last update

Re: [R] install.package(TinnR) - there is no package called 'TinnR'

2009-03-24 Thread Richardson, Patrick
Hi Martin, If all else fails, you could download the package from http://cran.r-project.org/web/packages/TinnR/index.html and install as a .zip file from within R. That is bizarre that R cannot find the package. I've had no problems downloading and installing. Best regards, Patrick

[R] CONFIDENCE INTERVAL FOR GLMER MODEL

2009-03-24 Thread LiatL
I've built a poisson regression model for multiple subjects by using the GLMER function. I've also developed some curves for defining its limits but I did not succeed in developing confidence interval for the model's curve (confint or predict does not work - only for glm). Does anyone know how

Re: [R] Variogram with Gstat

2009-03-24 Thread ONKELINX, Thierry
Have a look at the examples in the helpfile for variogram from the gstat package. data(meuse) # no trend: coordinates(meuse) = ~x+y variogram(log(zinc)~1, meuse) # residual variogram w.r.t. a linear trend: variogram(log(zinc)~x+y, meuse) # directional variogram: variogram(log(zinc)~x+y, meuse,

Re: [R] install.package(TinnR) - there is no package called 'TinnR'

2009-03-24 Thread Martin Tomko
Patrick, Romain, thank you very much for your help! I have found the site of the package at http://r-forge.r-project.org/projects/tinnr/ , as Romain suggests, but when you try to download, iyou find that the package is actually not contributed to the repository, it is just a space holder for it!

Re: [R] install.package(TinnR) - there is no package called 'TinnR'

2009-03-24 Thread Duncan Murdoch
On 3/24/2009 7:05 AM, Martin Tomko wrote: I have troubles make TinnR 2.2.0.2 work, it seems that the dependency on the package TinnR that cannot be found (I tried also manual downloads, but I cannot find the package anywhere on any CRAN mirror). What R version are you using, on what platform?

Re: [R] Calculating percentage Missing value for variables using one object

2009-03-24 Thread David Winsemius
It looks to me that you should be using the table or the xtabs function. You have apparently already decided not to use NA for missing values, so the instances in which variable1 == you should get counts with those functions: dft - data.frame(var1 = sample(c(, this, that, and), 120,

Re: [R] install.package(TinnR) - there is no package called 'TinnR' (RESOLVED)

2009-03-24 Thread Martin Tomko
Hi Duncan, it is Win Xp with the 2.6.0 R-project version. Sorry, I should have included this before. Installing the TinnR package manually from a local zip file downloaded from CRAN helped. I am still not sure why the package was not picked in the repositories. Can anyone please see if the

Re: [R] help: what are the basis functions in {mgcv}: gam?

2009-03-24 Thread oliviax
Thank you very much. I have downloaded your paper. Best regards. Olivia Simon Wood-4 wrote: If you can't get hold of the book, try appendix A of http://www.maths.bath.ac.uk/~sw283/simon/papers/psm.pdf -- what's described there is the cr basis best, Simon On Tuesday 24 March

[R] RES: Plot Means Line with Standard Deviation as Whiskers

2009-03-24 Thread Rodrigo Aluizio
Thanks for the answer, I'll take a look at these, meanwhile I forged some data as an example of what I'll have soon. Sam.Vol. A1R179,40 A1R2102,09 A1R3102,09 A2R1124,78 A2R2158,81 A2R3153,13 A3R1226,87 A3R2164,48 A3R3153,13 A4R1181,49 A4R2175,82

Re: [R] install.package(TinnR) - there is no package called 'TinnR' (RESOLVED)

2009-03-24 Thread Uwe Ligges
Martin Tomko wrote: Hi Duncan, it is Win Xp with the 2.6.0 R-project version. Sorry, I should have included this before. Installing the TinnR package manually from a local zip file downloaded from CRAN helped. I am still not sure why the package was not picked in the repositories. Can

Re: [R] CONFIDENCE INTERVAL FOR GLMER MODEL

2009-03-24 Thread Douglas Bates
On Tue, Mar 24, 2009 at 5:25 AM, LiatL liat.lam...@gmail.com wrote: I've built a poisson regression model for multiple subjects by using the GLMER function. I've also developed some curves for defining its limits but I did not succeed in developing confidence interval for the model's curve

Re: [R] Calculating percentage Missing value for variables using one object

2009-03-24 Thread Gabor Grothendieck
Read in the data, aggregate it by month and then turn it into a monthly zoo object and plot using a custom X axis: Lines - 'dos,variable1,variable2 May-06,1, May-06,2, June-06,,2 June-06,1,4 July-06,1,4 July-06,1,4 August-06,1,4 August-06,1,4' DF - read.table(textConnection(Lines), header = TRUE,

Re: [R] install.package(TinnR) - there is no package called 'TinnR' (RESOLVED)

2009-03-24 Thread Duncan Murdoch
On 3/24/2009 8:55 AM, Martin Tomko wrote: Hi Duncan, it is Win Xp with the 2.6.0 R-project version. Sorry, I should have included this before. Installing the TinnR package manually from a local zip file downloaded from CRAN helped. I am still not sure why the package was not picked in the

Re: [R] install.package(TinnR) - there is no package called 'TinnR' (RESOLVED)

2009-03-24 Thread Martin Tomko
Yes, will do. I just was not in need to upgrade as everything worked... An R is not my primary development environment - i.e., I need it a couple of times a year. But you are right, it is good to be on an updated version. Cheers Martin Duncan Murdoch wrote: On 3/24/2009 8:55 AM, Martin Tomko

Re: [R] Barplot With Selected X-Axis Region

2009-03-24 Thread Uwe Ligges
Gundala Viswanath wrote: Dear all, I have a data that looks like this: print(dat) V1 V2 1 -43342073 14 2 -433377304 3 -432846761 11372 75188572 11 11373 752061656 11374 75262720 24 What I want to do is to have a barplot where x-axis is

[R] Do we have to control for block in block designs if it is insignificant?

2009-03-24 Thread J S
I am wondering if in a block experimental design (ex. triple square lattice), the block effect was not significant, is it all right not to include the block effect in an empirical model (even though the sampling was done from different blocks)? Or we are forced to control for the block

Re: [R] How to find the path or the current file?

2009-03-24 Thread Marie Sivertsen
Thank you Gabor and Duncan for your replys. Mvh. Marie On Tue, Mar 24, 2009 at 1:00 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: See: https://stat.ethz.ch/pipermail/r-help/2009-January/184745.html On Tue, Mar 24, 2009 at 7:16 AM, Marie Sivertsen mariesiv...@gmail.com wrote:

[R] as.list cannot handle an environment

2009-03-24 Thread Wacek Kusnierczyk
?as.list states that as.list can be used with environments as arguments (what is referred to there as the 'environment method for as.list). indeed, e = globalenv() as.list(e) # $e # environment: R_GlobalEnv however, it seems to have problems with environment objects with a class

[R] Help with lattice

2009-03-24 Thread Steve_Friedman
Hello. I'm working on a windows XP machine with R 2.8.1 I'm working with lattice and I don't understand how to set the size of the text for the labels. Can anyone suggest a solution? I'm establishing a melted dataframe first and then setting a key and calling the plot as follows;

[R] using residuals of binomial GLM

2009-03-24 Thread Yuval Sapir
Hi all, This is more a question in statistics, but I hope to get also the R practice for my question: I have an ancova model where the response variable is flowering (plant has a flower = 1, no flower = 0). The explanatory variables are leaf length, leaf thick (both continuous variables), and

Re: [R] Random sampling based on the observations

2009-03-24 Thread Arup
Thank you Adrian. Its working.Thank you so much. Take care =^D=) Cheers Arup Pramanik Business Analyst Redwood Associates Bangalore,India. Adrian Dusa wrote: On Tuesday 24 March 2009, you wrote: Hello!I am having a problem with Random sampling in R. I have used a syntax:

[R] Has someone connected R with ITK?

2009-03-24 Thread Brandon . J . Whitcher
Dear R users, I am curious to know if someone has connected R with the Insight Segmentation and Registration Toolkit (www.itk.org). From the website... ITK is an open-source, cross-platform system that provides developers with an extensive suite of software tools for image analysis.

[R] Odp: error statement: missing value where TRUE/FALSE needed

2009-03-24 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 23.03.2009 15:38:49: Hi list, I want to try Gibbs sampling as a method of estimating a markov-switching model of a mean-deviating, pth-order autoregressive process with time varying transition probabilities via R and am using a code originally

Re: [R] Do we have to control for block in block designs if it is insignificant?

2009-03-24 Thread Liaw, Andy
The short answer is no (meaning to leave the blocks in the model). As Frank Harrell said, you've spent your degrees of freedom. Go home and be happy. Best, Andy From: J S Sent: Tuesday, March 24, 2009 9:49 AM To: r-help@r-project.org Subject: [R] Do we have to control for block in block

[R] how to do trend.test for subset data

2009-03-24 Thread Xin Shi
Dear: I try to do the trend.test {pastecs} for a subset data. Anyone have this experience? Many thanks! Xin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Hotteling Test

2009-03-24 Thread MarcioRibeiro
Hi listers, I would like to do a graphic of the power function of the Hotteling test (T2). T2~F(p,n-p;lambda), lambda is my paremeter of non-centrality. So I would like to do a graphic of the power of the test as a function of the parameter lambda. I also would like to define lambda varying

Re: [R] How to find the path or the current file?

2009-03-24 Thread William Dunlap
I sometimes deal with this problem by putting all the scripts into a package (in the inst/ directory if they don't fit elsewhere, perhaps in demo/) and then have them all use system.file(package=myPkg,...) to locate the related files. E.g., terms.conf -

Re: [R] how to do trend.test for subset data

2009-03-24 Thread David Winsemius
That question seems too vague to be answered. There should be no problem in doing such. Do you have a problem that you are encountering? -- David Winsemius On Mar 24, 2009, at 10:52 AM, Xin Shi wrote: Dear: I try to do the trend.test {pastecs} for a subset data. Anyone have this

Re: [R] confidence interval or error of x intercept of a linear

2009-03-24 Thread Peter Dalgaard
(Ted Harding) wrote: On 24-Mar-09 03:31:32, Kevin J Emerson wrote: ... When I have time for it (not today) I'll see if I can implement this neatly in R. It's basically a question of solving (N-2)*(1 - R(X0))/R(X0) = qf(P,1,(N-1)) for X0 (two solutions, maybe one, if any exist). etc. A

Re: [R] R: how to estimate multidimensional spectral measure of coherence

2009-03-24 Thread stephen sefick
Look at the function wco to see if you can use it. I don't have R on my work computer right now to test and see, but AI think you could look at this. I don't know what a slidding window is, but wavelets are time localized so I don't think this would be a problem. Also, you could easily slide

Re: [R] r online

2009-03-24 Thread j verzani
Dirk Eddelbuettel edd at debian.org writes: On 24 March 2009 at 12:06, Thomas Steiner wrote: | Hi, | I'd like to execute simple commands and functions in R through a | website, is there any service like this somewhere? | I only found http://www.osvisions.com/r-online/ but it does not

[R] Display Equation on plot

2009-03-24 Thread Douglas M. Hultstrand
Hello, I am using the lm to fit a linear model to data, I was wondering if there is a way to display the equation on a plot using the extracted lm coefficients? I am using the plot() function to create the plot/.png. Example: lm_mod - lm(data1~data2) intercept = 48.54 slope = 0.4856 I

Re: [R] Display Equation on plot

2009-03-24 Thread Uwe Ligges
Douglas M. Hultstrand wrote: Hello, I am using the lm to fit a linear model to data, I was wondering if there is a way to display the equation on a plot using the extracted lm coefficients? I am using the plot() function to create the plot/.png. Example: lm_mod - lm(data1~data2)

[R] no snowfall on debian after lenny upgrade

2009-03-24 Thread Martin Renner
I feel like I've tried the works on something that should be obvious but I'm stuck. Package snowfall used to run just great under debian etch (amd-64, 2x quad-core). After upgrading to lenny (debian 5, reinstalled from scratch), I haven't been able to get snowfall running again.

Re: [R] Display Equation on plot

2009-03-24 Thread Duncan Murdoch
On 3/24/2009 12:50 PM, Douglas M. Hultstrand wrote: Hello, I am using the lm to fit a linear model to data, I was wondering if there is a way to display the equation on a plot using the extracted lm coefficients? I am using the plot() function to create the plot/.png. Example: lm_mod -

Re: [R] Display Equation on plot

2009-03-24 Thread Matthieu Dubois
Douglas M. Hultstrand dmhultst at metstat.com writes: Hello, I am using the lm to fit a linear model to data, I was wondering if there is a way to display the equation on a plot using the extracted lm coefficients? I am using the plot() function to create the plot/.png. Example:

[R] Testing for symbolic links

2009-03-24 Thread Gregoire Pau
Hello, How can I test in R whether a file is a symbolic link or not ? I haven't found anything relevant in R-help but I may be wrong. Thanks ! Greg --- Gregoire Pau Research Officer EMBL Heidelberg, Germany http://www.ebi.ac.uk/~gpau/ __

[R] generalized Poisson regression

2009-03-24 Thread Arafat Ud Zaman
Dear sir, I want to know about R command for parameter estimation for generalized Poisson regression. yours faithfully Arafat Ud Zaman 4th year(Hons) Applied statistics University of Dhaka Mob- +8801715173176 E-mail: arafat.za...@yahoo.com azam...@isrt.ac.bd

[R] Workspace Image

2009-03-24 Thread Yikan Liu
Hello, I'm wondering how to delete previously saved workspace in R. Thank you. Regards, Yikan Liu [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] Why na.rm=FALSE is the default

2009-03-24 Thread Adam D. I. Kramer
Dear Colleagues, I've been searching for a post or article or something which explains why having na.rm=FALSE or na.action=na.fail as the default is a better choice than TRUE or na.omit. I understand the basic argument: it does not make sense to average a nonexistance into an

Re: [R] generalized Poisson regression

2009-03-24 Thread Rob Denniker
Not 100% sure what you are looking for, but have a look at the Generalized Event Count model in the Zelig package. It will also let you fit a Poisson and other event count models by MLE. -Original Message- From: azam...@isrt.ac.bd Sent: Tue, 24 Mar 2009 23:26:59 +0600 To:

[R] Search Errors

2009-03-24 Thread CE.KA
Hi R users, I have a big program So in Rgui I can t see all the execution of it Is there a way to ask R if there is Errors in my program Sincerely yours -- View this message in context: http://www.nabble.com/Search-%22Errors%22-tp22687080p22687080.html Sent from the R help mailing list archive

Re: [R] Workspace Image

2009-03-24 Thread Dieter Menne
Yikan Liu u4297812 at anu.edu.au writes: I'm wondering how to delete previously saved workspace in R. file.remove(mywork.rdata) Dieter __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Summarizing each row into a frequency table

2009-03-24 Thread William Dunlap
Have you tried table(row(m), m) ? If there is a chance that m doesn't have at least 1 of each of c(-1,0,1) then force them into the table with table(row(m), factor(m, levels=c(-1,0,1))) Bill Dunlap TIBCO Software Inc - Spotfire Division wdunlap tibco.com

Re: [R] Why na.rm=FALSE is the default

2009-03-24 Thread Greg Snow
The only reference that I can think of (a bit subtle/indirect) is: http://www.burns-stat.com/pages/Tutor/spreadsheet_addiction.html (look in the section on Propagation of Blanks). But I think that it really comes down to the following 2 variations on a rule: 1. Important decisions (such as

Re: [R] read in large data file (tsv) with inline filter?

2009-03-24 Thread David Reiss
Hi all, Thanks a lot for your responses. I forgot to mention in my posting that I also want to try to be as cross-platform as possible (and thus to avoid relying on external calls to UNIX programs such as grep). I like the idea from Thomas: You certainly don't want to use repeated reads from the

[R] modelling probabilities instead of binary data with logistic regression

2009-03-24 Thread joris meys
Dear all, I have a dataset where I reduced the dimensionality, and now I have a response variable with probabilities/proportions between 0 and 1. I wanted to do a logistic regression on those, but the function glm refuses to do that with non-integer values in the response. I also tried lrm, but

[R] Probability trees

2009-03-24 Thread Waldir Leôncio
How can I create probability trees in R? I've searched everywhere I could and still have no clue. :-/ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] segfault when running heatmap()

2009-03-24 Thread Andrew Yee
Hi, I was wondering if someone in the mailing list has any insight into this segfault error that I consistently find when running a script containing heatmap() in R 2.8.1 and 2.8.0 on a Linux 64-bit machine. Some points: 1. This occurs when running heatmap(). 2. Interestingly, if I source() the

Re: [R] Testing for symbolic links

2009-03-24 Thread Rolf Turner
On 25/03/2009, at 6:29 AM, Gregoire Pau wrote: Hello, How can I test in R whether a file is a symbolic link or not ? I haven't found anything relevant in R-help but I may be wrong. Thanks ! Since symbolic links are a Unix thing, I assume you are working on a system that is

Re: [R] modelling probabilities instead of binary data with logistic regression

2009-03-24 Thread joris meys
Thank you all for the very fast answers. My proportions come from a factor analysis on a number of binary variables, in order to avoid having to fit 12 logistic regressions on the same dataset. By scaling the obtained scores to 0 and 1, I get weighted averages of the response combinations I'm

[R] Green and Byar (1980) Prostate Cancer Data set from Andrews and Herzberg - Data

2009-03-24 Thread Ravi Varadhan
Hi, I am looking for a data set containing the information from a randomized trial evaluating the effect of DES (diethylsilbestrol) on multiple time-to-event endpoints, prostate cancer, CVD, and other causes. The original source of this data is Green and Byar (1980). This is a popular

[R] Executing an external executable from within R?

2009-03-24 Thread Jason Rupert
For example, from within an R script I would like to be able to launch NotePad.exe. Is there an R call for launching Notepad.exe from within R? I've tried the following, but they don't seem to work: system(paste('C:/WINDOWS/NOTEPAD.exe'), wait = FALSE) system(C:/WINDOWS/NOTEPAD.exe, wait

  1   2   >