[R] rpart.object help

2010-12-13 Thread jagdeesh_mn
Hi, Suppose i have generated an object using the following : fit - rpart(Kyphosis ~ Age + Number + Start, data=kyphosis) And when i print fit, i get the following : n= 81 node), split, n, loss, yval, (yprob) * denotes terminal node 1) root 81 17 absent (0.7901235 0.2098765) 2)

Re: [R] rpart.object help

2010-12-13 Thread Prof Brian Ripley
On Sun, 12 Dec 2010, jagdeesh_mn wrote: Hi, Suppose i have generated an object using the following : fit - rpart(Kyphosis ~ Age + Number + Start, data=kyphosis) And when i print fit, i get the following : n= 81 node), split, n, loss, yval, (yprob) * denotes terminal node 1) root 81 17

[R] descriptive statistics

2010-12-13 Thread effeesse
Hi. In a data set I have a variable that takes values from 1 to 14. For each subgroup of values of this variable, I would like to obtain some descriptive statistics of other variables present in the data set. I've been trying with a for loop but I couldn't get nothing. Could you please suggest me

Re: [R] descriptive statistics

2010-12-13 Thread Ivan Calandra
?aggregate ?doBy::summaryBy Le 12/13/2010 11:04, effeesse a écrit : Hi. In a data set I have a variable that takes values from 1 to 14. For each subgroup of values of this variable, I would like to obtain some descriptive statistics of other variables present in the data set. I've been trying

Re: [R] package sampling

2010-12-13 Thread Uwe Ligges
Which version of R and sampling? Where is the reproducible code? Uwe Ligges On 11.12.2010 16:16, andrija djurovic wrote: Hi R users. I have a problem with function strata in sampling packages. st0 = strata(dom, stratanames=stratas, size=sample.size, method=systematic,pik, FALSE) Error in

Re: [R] Plot's aspect ratio and pty

2010-12-13 Thread Uwe Ligges
It does for me under R-2.12.0 32-bit on Windows with the windows() device, so: Which version of R, which OS, which device do you use? Uwe Ligges On 13.12.2010 07:30, Marcin Kozak wrote: Dear All, I've been playing with pty, and it seems it does not produce square plots as it is expected

Re: [R] descriptive statistics

2010-12-13 Thread Jim Lemon
On 12/13/2010 09:04 PM, effeesse wrote: Hi. In a data set I have a variable that takes values from 1 to 14. For each subgroup of values of this variable, I would like to obtain some descriptive statistics of other variables present in the data set. I've been trying with a for loop but I

Re: [R] Plot's aspect ratio and pty

2010-12-13 Thread Marcin Kozak
R-11.1 on both 32-bit and 64-bit on Windows with the windows() device. Best Marcin 2010/12/13 Uwe Ligges lig...@statistik.tu-dortmund.de: It does for me under R-2.12.0 32-bit on Windows with the windows() device, so: Which version of R, which OS, which device do you use? Uwe Ligges

[R] Re : descriptive statistics

2010-12-13 Thread justin bem
A nice way to obtain summary for data is to use summary.formula in Hmisc package. Justin BEM BP 1917 Yaoundé Tél (237) 76043774 De : Jim Lemon j...@bitwrit.com.au À : effeesse scarpin...@gmail.com Cc : r-help@r-project.org Envoyé le : Lun 13 décembre 2010,

Re: [R] Summary (Re: (S|odf)weave : how to intersperse (\LaTeX{}|odf) comments in source code ? Delayed R evaluation ?)

2010-12-13 Thread Claudia Beleites
Dear Emmanuel and dear list, Therefore, I let this problem to sleep. However, I Cc this answer (with the original question below) to Max Kuhn and Friedrich Leisch, in the (faint) hope that this feature, which does not seem to have been missed by anybody in 8 years, I've been missing it every

[R] How to compare two square matrices

2010-12-13 Thread Sabari
Hi, I have two matrices containing some probabilities score obtained from two different prediction programs. Now, I want to compare these two matrices to measure the difference. Could you please suggest some method to do this in R. Thanks Sabari [[alternative HTML version deleted]]

Re: [R] descriptive statistics

2010-12-13 Thread David Hajage
Another way is the remix function of the remix package. On Monday, December 13, 2010, justin bem justin_...@yahoo.fr wrote: A nice way to obtain summary for data is to use summary.formula in Hmisc package.  Justin BEM BP 1917 Yaoundé Tél (237) 76043774

Re: [R] Plot's aspect ratio and pty

2010-12-13 Thread Uwe Ligges
On 13.12.2010 11:29, Marcin Kozak wrote: R-11.1 This one does not exist. Please try R-2.12.0 (but it also worked with R-2.11.1 if you meant that). My guess is that you are confusing plotting region with device region. In order to get a squared device region, you have to ask the device

[R] Browsing through a dataframe page by page (like with shell command more)

2010-12-13 Thread Alexandre CESARI
Hello, I'm looking for an easy way to display a data.frame (or other variables) page by page, similarly to what is possible on a file using the more command in a standard UNIX shell. Any help would be greatly appreciated. Thanks Alexandre [[alternative HTML version deleted]]

Re: [R] Plot's aspect ratio and pty

2010-12-13 Thread Prof Brian Ripley
On Mon, 13 Dec 2010, Uwe Ligges wrote: It does for me under R-2.12.0 32-bit on Windows with the windows() device, so: Which version of R, which OS, which device do you use? Since (s)he used windows() we know the OS. But I think one possible explanation is on the help page, arguments

Re: [R] overlap different line in a xyplot (lattice)

2010-12-13 Thread Francesco Nutini
From: fe...@nfrac.org Date: Sun, 12 Dec 2010 11:47:55 +1100 Subject: Re: [R] overlap different line in a xyplot (lattice) To: ehl...@ucalgary.ca CC: nutini.france...@gmail.com; r-help@r-project.org On 12 December 2010 00:08, Peter Ehlers ehl...@ucalgary.ca wrote: On 2010-12-11 03:12,

Re: [R] rpart.object help

2010-12-13 Thread jagdeesh_mn
Prof Brian Ripley wrote: On Sun, 12 Dec 2010, jagdeesh_mn wrote: Hi, Suppose i have generated an object using the following : fit - rpart(Kyphosis ~ Age + Number + Start, data=kyphosis) And when i print fit, i get the following : n= 81 node), split, n, loss, yval, (yprob) *

[R] Problem with retrieve.nc of clim.pact

2010-12-13 Thread Pablo85
Goodmorning to everyone, I'm new so sorry for bad english and formulation. I go to the point: I am using since months the retrieve.nc function from clim.pact to extract data from netcdf files. I always had no problem, but some days ago I updated both R and the package to the latest versione

[R] post-hoc test for ANCOVA method

2010-12-13 Thread Francesco Nutini
Dear [R] Users, I have implemented a linear model with this syntax: model- lm (var_dependent ~ var_indipendent + factor + var_indipendent : factor, dataframe) anova (model) Response: var_dependent Df Sum Sq Mean Sq F value Pr(F)

[R] Wrong contrast matrix for nested factors in lm(), rlm(), and lmRob()

2010-12-13 Thread Saul Sternberg
This message also reports wrong estimates produced by lmRob.fit.compute() for nested factors when using the correct contrast matrix. And in these respects, I have found that S-Plus behaves the same way as R. Using the three available contrast types (sum, treatment, helmert) with lm() or

Re: [R] simple plotting question

2010-12-13 Thread Anthony Damico
to add to Michael's response: http://www.statmethods.net/advgraphs/parameters.html On Mon, Dec 13, 2010 at 2:23 AM, Michael Bedward michael.bedw...@gmail.com wrote: Hello Erin, Try this... plot(x, y, type=b, pch=16) Michael On 13 December 2010 18:11, Erin Hodgess

Re: [R] Why do we have to turn factors into characters for various functions?

2010-12-13 Thread Heinz Tuechler
Hello Petr, don't want to convince you. If you like the following: x - factor(1:4, labels=c(one, two, three, four)) y - factor(3:5, labels=c(three, four, five)) data.frame(character=c(as.character(x), as.character(y)), numeric=c(x, y)) character numeric 1 one 1 2 two

Re: [R] Need help on nnet

2010-12-13 Thread Georg Ruß
On 10/12/10 02:56:13, jothy wrote: Am working on neural network. Below is the coding and the output [...] summary (uplift.nn) a 3-3-1 network with 16 weights options were - b-h1 i1-h1 i2-h1 i3-h1 16.646.62 149.932.24 b-h2 i1-h2 i2-h2 i3-h2 -42.79 -17.40 -507.50

Re: [R] Help..Neural Network

2010-12-13 Thread Georg Ruß
On 10/12/10 03:45:46, sadanandan wrote: I am trying to develop a neural network with single target variable and 5 input variables to predict the importance of input variables using R. I used the packages nnet and RSNNS. But unfortunately I could not interpret the out put properly and the

[R] help on SAS Macro in R

2010-12-13 Thread Özgür Asar
Dear Researchers, I am looking for to read a SAS macro in R. Although I searched it on web, I couldn’t find anything. Can you help me or direct me? Thank you for your interest and patience. Best. Ozgur [[alternative HTML version deleted]]

Re: [R] spatial clusters

2010-12-13 Thread Georg Ruß
On 10/12/10 23:26:28, dorina.lazar wrote: I am looking for a clustering method usefull to classify the countries in some clusters taking account of: a) the geographical distance (in km) between countries and b) of some macroeconomic indicators (gdp, life expectancy...). Hi Dorina, before

Re: [R] Browsing through a dataframe page by page (like with shell command more)

2010-12-13 Thread Peter Ehlers
On 2010-12-13 02:49, Alexandre CESARI wrote: Hello, I'm looking for an easy way to display a data.frame (or other variables) page by page, similarly to what is possible on a file using the more command in a standard UNIX shell. Any help would be greatly appreciated. How about View(mydata)?

Re: [R] help on SAS Macro in R

2010-12-13 Thread Erik Iverson
On 12/13/2010 07:14 AM, Özgür Asar wrote: Dear Researchers, I am looking for to read a SAS macro in R. Although I searched it on web, I couldn’t find anything. Are you hoping just to read it in, or to actually have it execute the macro as SAS would? What gives you the idea the latter is

Re: [R] help on SAS Macro in R

2010-12-13 Thread David Winsemius
On Dec 13, 2010, at 8:14 AM, Özgür Asar wrote: Dear Researchers, I am looking for to read a SAS macro in R. Although I searched it on web, I couldn’t find anything. Can you help me or direct me? fortune(SAS) For almost 40 years SAS has been the primary tool for statisticians

Re: [R] overlap different line in a xyplot (lattice)

2010-12-13 Thread Peter Ehlers
On 2010-12-13 03:13, Francesco Nutini wrote: From: fe...@nfrac.org [...snip...] xyplot(y1 + y2 ~ x | sites, DF, type = b) Great Felix! this is what I was looking for! But if y1 and y2 have a different scales? Can I plot, for example y2, on secondary axis? There are probably good

Re: [R] Re : descriptive statistics

2010-12-13 Thread effeesse
I am sorry, but I cannot understand how to use the summary function. Maybe, if I describe my needs, you could sketch a line that could work. In the data set variable V can take values 1 to 14. For the subgroup of individuals where V takes value =1 I want the mean and variance of a certain set of

[R] How does R compute sums of squares?

2010-12-13 Thread Ethan Arenson
Consider the following missing data problem: y = c(1, 2, 2, 2, 3) a = factor(c(1, 1, 1, 2, 2)) b = factor(c(1, 2, 3, 1, 2)) fit = lm(y ~ a + b) anova(fit) Analysis of Variance Table Response: y Df Sum Sq Mean SqF valuePr(F) a 1 0.8 0.8 1.3637e+33 2.2e-16

Re: [R] Re : descriptive statistics

2010-12-13 Thread Ivan Calandra
I would suggest what we already suggested to you: ?aggregate ?by ?doBy::summaryBy We could help you more precisely if you could provide a reproducible example, as explained in the posting guide (see link at the end of every email from this list) Ivan Le 12/13/2010 15:14, effeesse a écrit :

Re: [R] How to compare two square matrices

2010-12-13 Thread Arun.stat
Dear Sabari, if you need a single number for comparison then there could be many options. You can calculate smallest absolute eigen value, or may be the determinant (i.e. a measure of volumn of matrices) or may be the smallest element in absolute term, depending on your research need. Thanks, --

Re: [R] help requested

2010-12-13 Thread profaar
hi thanks for your suggestion and reply. let me try it out. With Warm Wishes and Regards A. Abdul Rasheed, M.C.A., M.E., Ph.D., Assistant Professor, Department of Computer Applications, Valliammai Engineering College, SRM Nagar,

[R] Re : Re : descriptive statistics

2010-12-13 Thread justin bem
With summary do this my.summary-function(x) c(mean(x),var(x)) summary(v1~V, fun=my.summary,data=df) summary(v2~V, fun=my.summary,data=df) summary(v3~V, fun=my.summary,data=df) summary(v4~V, fun=my.summary,data=df) summary(v5~V, fun=my.summary,data=df) If you want you get the mean of all

[R] check for item in vector

2010-12-13 Thread C.H.
Dear R users, Suppose I have an vector like this: animal - c(Tiger,Panda) I would like to know is there any function that check for the existence of certain item in a vector. e.g. func(Tiger,animal) # check for the existence of Tiger TRUE func(Acacia,animal) #Acacia is not an item of the

Re: [R] check for item in vector

2010-12-13 Thread Ivan Calandra
Hi, See ?%in% or ?match animal - c(Tiger,Panda) Tiger %in% animal [1] TRUE Acacia %in% animal [1] FALSE Panda %in% animal [1] TRUE HTH, Ivan Le 12/13/2010 15:48, C.H. a écrit : Dear R users, Suppose I have an vector like this: animal- c(Tiger,Panda) I would like to know is there any

[R] Reg : null values in kmeans

2010-12-13 Thread Raji
Hi, I am using k means algorithm for clustering.My data contains a few null/NA values.kmeans doesnt cluster with those values.Are there any option like na.omit which can avoid these null values and cluster the remaining values? Thanks, Raji -- View this message in context:

Re: [R] check for item in vector

2010-12-13 Thread Jorge Ivan Velez
Hi CH, Check ?is.element ?%in% HTH, Jorge On Mon, Dec 13, 2010 at 9:48 AM, C.H. wrote: Dear R users, Suppose I have an vector like this: animal - c(Tiger,Panda) I would like to know is there any function that check for the existence of certain item in a vector. e.g.

Re: [R] check for item in vector

2010-12-13 Thread David L Lorenz
CH, How about any: any(Tiger == animal) The function which will tell you the index if any match which(Tiger == animal. You should also look at the match funciton. Dave From: C.H. chainsawti...@gmail.com To: R-help r-help@r-project.org Date: 12/13/2010 08:50 AM Subject: [R] check for

[R] booklet on using R for biomedical statistics

2010-12-13 Thread Coghlan, Avril
Dear all, I've written a small booklet on using R for biomedical statistics (mostly focussed on cohort and case-control studies), available here under a Creative Commons Attribution 3.0 License : A Little Book of R for Biomedical Statistics

[R] [R-pkgs] batchfiles 0.6-0

2010-12-13 Thread Gabor Grothendieck
batchfiles is a set of batch, javascript and HTML Application files that are useful for running R and associated programs on Windows. Version 0.6-0 updates them for the new architecture specific directory structure in R 2.12.0 . A few of the lesser used utilities have been dropped. Each

[R] How to plot Ellipsoid like function

2010-12-13 Thread Uwe Wolfram
Dear R-Users, I am currently trying to fit a tensorial function in its principal coorinate system. The function is given by: 1~(x1^2 + x2^2 + x3^2 - chi0*(x1*x2 + x1*x3 + x2*x3))/eps0^2 + (x1 + x2 + x3)/xi0 Where eps0 = 0.0066, chi0 = -0.66 and xi0 = 0.011 are obtained from experimental data

Re: [R] Rapache on windows

2010-12-13 Thread Bos, Roger
I asked Jeff Horner that question a while back and he said it was Linux only. He doesn't have time to create a windows version. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Santosh Srinivas Sent: Saturday, December 11, 2010

Re: [R] How to plot Ellipsoid like function

2010-12-13 Thread Duncan Murdoch
On 13/12/2010 10:13 AM, Uwe Wolfram wrote: Dear R-Users, I am currently trying to fit a tensorial function in its principal coorinate system. The function is given by: 1~(x1^2 + x2^2 + x3^2 - chi0*(x1*x2 + x1*x3 + x2*x3))/eps0^2 + (x1 + x2 + x3)/xi0 Where eps0 = 0.0066, chi0 = -0.66 and xi0 =

Re: [R] Tukey HSD not working

2010-12-13 Thread PGZC
TukeyHSD(fit) Error in TukeyHSD(fit) : object 'fit' not found TukeyHSD(anova) Error in UseMethod(TukeyHSD) : no applicable method for 'TukeyHSD' applied to an object of class function TukeyHSD(group) Error in UseMethod(TukeyHSD) : no applicable method for 'TukeyHSD' applied to an object

[R] Integration with LaTex and LyX

2010-12-13 Thread Eduardo de Oliveira Horta
Hello, Are there any packages which allow for a good integration between R and LaTex / LyX? I'm interested mainly in automatic (automagic?) imports of plots/graphics. Thanks in advance and best regards, Eduardo de Oliveira Horta [[alternative HTML version deleted]]

Re: [R] Tukey HSD not working

2010-12-13 Thread PGZC
TukeyHSD(fit) Error in UseMethod(TukeyHSD) : no applicable method for 'TukeyHSD' applied to an object of class lm TukeyHSD(anova) Error in UseMethod(TukeyHSD) : no applicable method for 'TukeyHSD' applied to an object of class function TukeyHSD(group) Error in UseMethod(TukeyHSD) : no

Re: [R] ggplot2 histograms

2010-12-13 Thread Hadley Wickham
Hi Sandy, The way I'd describe it is that you expected the width parameter of the position adjustment to be relative to the binwidth of the histogram - but it's actually absolute, and it has to be this way because there's currently no way for the position adjustment to know about the parameters

[R] How to ignore data

2010-12-13 Thread Steve Sidney
Dear list I have quite a small data set in which I need to have the following values ignored - not used when performing an analysis but they need to be included later in the report that I write. Can anyone help with a suggestion as to how this can be accomplished Values to be ignored 0 -

[R] Testing an interaction with a random effect in lmer

2010-12-13 Thread Reed, Nicola
Hi, I was hoping to get some advice regarding the testing of interactions, when one factor is modelled as a random effect... I have a model with binomial error structure where the response variable is the proportion of time spent at the main sett (animals were tracked for 28 consecutive days

Re: [R] Integration with LaTex and LyX

2010-12-13 Thread Jonathan P Daily
?Sweave LyX is a bit harder, although you can probably export LyX docs to a *.tex and Sweave those fairly painlessly. -- Jonathan P. Daily Technician - USGS Leetown Science Center 11649 Leetown Road Kearneysville WV, 25430 (304) 724-4480 Is the room still a

[R] odot symbol as a subscript in axes labels

2010-12-13 Thread pilchat
Dear R users, do you know how to print the latex \odot symbol subscripted to axes labels? Thank you in advance Gaetano __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Re : Re : descriptive statistics

2010-12-13 Thread effeesse
what am I supposed to put into function(x)? The indicator for extracting the subgroups? data is the df. cluster={1,...,14}. This is how I was compiling: for (i in 1:14) { my.summary-data$cluster==i c(mean(?),var(?)) summary(var_A~cluster, fun=my.summary,data=data) summary(var_B~cluster,

Re: [R] Integration with LaTex and LyX

2010-12-13 Thread Liviu Andronic
On Mon, Dec 13, 2010 at 4:55 PM, Jonathan P Daily jda...@usgs.gov wrote: ?Sweave LyX is a bit harder, although you can probably export LyX docs to a *.tex and Sweave those fairly painlessly. LyX can play very nicely with Sweave and R. For the 1.6.x series, you could get started here [1]. If

Re: [R] Tukey HSD not working

2010-12-13 Thread Peter Ehlers
On 2010-12-13 07:29, PGZC wrote: TukeyHSD(fit) Error in UseMethod(TukeyHSD) : no applicable method for 'TukeyHSD' applied to an object of class lm TukeyHSD(anova) Error in UseMethod(TukeyHSD) : no applicable method for 'TukeyHSD' applied to an object of class function TukeyHSD(group)

Re: [R] How does R compute sums of squares?

2010-12-13 Thread Douglas Bates
On Mon, Dec 13, 2010 at 8:20 AM, Ethan Arenson ethan.a.aren...@gmail.com wrote: Consider the following missing data problem:  y = c(1, 2, 2, 2, 3) a = factor(c(1, 1, 1, 2, 2)) b = factor(c(1, 2, 3, 1, 2)) fit = lm(y ~ a + b) anova(fit)  Analysis of Variance Table Response: y          

Re: [R] How to ignore data

2010-12-13 Thread Ben Bolker
Steve Sidney sbsidney at mweb.co.za writes: Dear list I have quite a small data set in which I need to have the following values ignored - not used when performing an analysis but they need to be included later in the report that I write. Can anyone help with a suggestion as to how

[R] Data permutation

2010-12-13 Thread matteop
Hello R User, I am new in R and trying to migrate from SAS. I have to convert a table that look like this YEARFIRMID_NAME VALUE 1994Microsoft John Doe5 1994Microsoft Mark Smith 3 1994Microsoft David Ring

[R] library(Defaults) throws evaluation nested too deeply

2010-12-13 Thread nbest937
Can anybody tell me why this is happening? library( Defaults) ls() [1] class foomlct [4] mlctTheme overlayFunctionpct [7] priprimaryFractionreclassFractions [10] reclassMatrix reclassNames

Re: [R] Re : Re : descriptive statistics

2010-12-13 Thread William Revelle
An alternative way of getting summary statistics by a grouping variable is to use describe.by in the psych package: using Jim Lemon's example: library(psych) testmat-data.frame(sample(1:14,50,TRUE),rnorm(50),runif(50)) #make up the data describe.by(test.mat,testmat[1]#get descriptive

Re: [R] Re : Re : descriptive statistics

2010-12-13 Thread Ivan Calandra
Do it with aggregate(), something like this should do: aggregate(.~cluster, FUN=summary, data=data) Now if you don't want to run summary(), replace it with the function you'd like. HTH, Ivan Le 12/13/2010 17:17, effeesse a écrit : what am I supposed to put into function(x)? The indicator

[R] 15' lag of an irregular Time Series

2010-12-13 Thread Vassilios Dimitrakas
Hi everyone, I am new to R and I have a beginner's question on Time Series: I have an irregular time series that goes like this: TIMESTAMP PRICE 2010-11-29 12:29:28 25.255 2010-11-29 12:30:47 25.255 2010-11-29 12:36:58 25.230 2010-11-29 12:43:14 25.235

Re: [R] How to ignore data

2010-12-13 Thread Bert Gunter
Values to be ignored 0 - zero and 1 this is in addition to NA (null) The reason is that I need to use the log10 of the values when performing the calculation. Currently I hand massage the data set, about a 100 values, of which less than 5 to 10 are in this category. This is probably a

Re: [R] Data permutation

2010-12-13 Thread Gabor Grothendieck
On Mon, Dec 13, 2010 at 11:37 AM, matteop mpr...@iese.edu wrote: Hello R User, I am new in R and trying to migrate from SAS. I have to convert a table that look like this YEAR    FIRM            ID_NAME         VALUE 1994    Microsoft       John Doe                5 1994    Microsoft      

[R] Qs re writing/reading arrays, dataframes

2010-12-13 Thread Roy Shimizu
Hi! I'm just getting started with R (and with the analysis of large datasets in general). I have several beginner-level questions whose answers I have not been able to find, and was hoping one of you would be kind enough to throw me a cluebrick or two. I have a 6-dimensional numeric array

Re: [R] Rapache on windows

2010-12-13 Thread Santosh Srinivas
Thanks to R ... I just got myself a new ubuntu setup just an hour back! It feels good! :-) On Mon, Dec 13, 2010 at 8:46 PM, Bos, Roger roger@rothschild.comwrote: I asked Jeff Horner that question a while back and he said it was Linux only. He doesn't have time to create a windows version.

Re: [R] How to ignore data

2010-12-13 Thread Steve Sidney
Thanks for the questions. 1) The data represents micro-organism counts and a count of zero in this case is highly unlikely given the info we have; including the other participants. 2) The data is submitted in duplicate and then a standardised sum and difference is established and is used to

Re: [R] Projecting data on a world map using long/lat

2010-12-13 Thread mathijsdevaan
Hi Patrick, Thanks! That worked perfectly! M -- View this message in context: http://r.789695.n4.nabble.com/Projecting-data-on-a-world-map-using-long-lat-tp3081298p3085834.html Sent from the R help mailing list archive at Nabble.com. __

[R] How to change leaf color by group in hclust plot or how to install A2R package in windows?

2010-12-13 Thread Soyeon Kim
I want to change leaf color by group in hclust plot. I've seen several answers about A2R package but I cannot install A2R and Rtools in windows. Do you know how to install A2R package in windows or how to change leaf color by group in hclust plot? Thank you in advance, Soyeon

[R] Complicated nls formula giving singular gradient message

2010-12-13 Thread Jared Blashka
I'm attempting to calculate a regression in R that I normally use Prism for, because the formula isn't pretty by any means. Prism presents the formula (which is in the Prism equation library as Heterologous competition with depletion, if anyone is curious) in these segments: KdCPM =

Re: [R] How to ignore data

2010-12-13 Thread Steve Sidney
Thanks for the comments Please see my reply to Stavros - the counts represent organisms and btw both mean and the median are virtually unaffected by the removal of these valuse. Furthermore, experience rather than statistics indicates that these values are in fact gross errors and as you of

Re: [R] How to change leaf color by group in hclust plot or how to install A2R package in windows?

2010-12-13 Thread Tal Galili
What error do you get when using: install.packages(A2R) ? 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

Re: [R] How to print colorful R output??

2010-12-13 Thread Greg Snow
One possibility, though not as simple as what you ask for, is to use etxtStart and friends from the TeachingDemos package. Other possibilities include using gui interfaces to R, possibilities (though they may do more than you ask, and color might be different) include emacs/ess; vim; jgr; and

[R] predict.lm[e] with formula passed as a variable

2010-12-13 Thread Thaler, Thorn, LAUSANNE, Applied Mathematics
Dear all, In a function I paste a string and convert it to a formula which I pass to lm[e]. The idea is to write a function which takes the name of the response variable and the explanatory variable and the data frame as an argument and calculates an lm[e]. (see example below) This works fine,

Re: [R] How to ignore data

2010-12-13 Thread Bert Gunter
Inline below. -- Bert On Mon, Dec 13, 2010 at 9:42 AM, Steve Sidney sbsid...@mweb.co.za wrote: Thanks for the questions. 1) The data represents micro-organism counts and a count of zero in  this case is highly unlikely given the info we have; including the other participants. ?? Censoring

Re: [R] How to change leaf color by group in hclust plot or how to install A2R package in windows?

2010-12-13 Thread Bryan Hanson
An example is described here that you can adapt: http://r.789695.n4.nabble.com/coloring-leaves-in-a-hclust-or-dendrogram-plot -tt795496.html#a795497 HTH. Bryan * Bryan Hanson Professor of Chemistry Biochemistry DePauw University, Greencastle IN USA On 12/13/10 12:54 PM, Soyeon

Re: [R] Qs re writing/reading arrays, dataframes

2010-12-13 Thread David Winsemius
On Dec 13, 2010, at 12:20 PM, Roy Shimizu wrote: Hi! I'm just getting started with R (and with the analysis of large datasets in general). I have several beginner-level questions whose answers I have not been able to find, and was hoping one of you would be kind enough to throw me a

[R] curve

2010-12-13 Thread Val
Hi All, I generated 5000 samples using the following script test- rnorm(5000,1000,100) test1 - subset(test, subset=(test 1100)) d - density(test) plot(d, main=Density of production) abline(v=mean(test1) I wanted to do the following but faced

Re: [R] Browsing through a dataframe page by page (like with shell command more)

2010-12-13 Thread Greg Snow
For data frames the best is probably the View function (note capitol V) which opens the data frame in a spreadsheet like window that you can scroll through. For more complicated, list or list-like objects, look at TkListView in the TeachingDemos package. For more general investigation of data

Re: [R] survival package - calculating probability to survive a given time

2010-12-13 Thread Terry Therneau
- --- included message i try to calculate the probabilty to survive a given time by using the estimated survival curve by kaplan meier. What is the right way to do that? as far as is see i cannot use the predict-methods from the survival package? end inclusion The

Re: [R] curve

2010-12-13 Thread Sarah Goslee
Here's one way to do what I think you want: test- rnorm(5000,1000,100) test1 - subset(test, subset=(test 1100)) d - density(test) plot(d, main=Density of production, xlab=) lines(d$x[d$x 1100], d$y[d$x 1100], col=blue, lwd=2) curveheight - d$y[abs((d$x -

Re: [R] predict.lm[e] with formula passed as a variable

2010-12-13 Thread David L Lorenz
Thorn, Here's how I do it: retval - list(as.name('lm'), formula=as.formula(paste(Response, ~, Explan, sep='')), data=as.name(Data)) #... optionally add other arguments retval - eval(as.call(retval)) Dave From: Thaler, Thorn, LAUSANNE, Applied

Re: [R] curve

2010-12-13 Thread Ashta
Thanks Sarah, 1. to shade or color (blue) the curve using the criterion that any values greater than 11,000 I think I was not clear in the above point. I want shade not the line but the area under the curve, and Your last line of code, segments(x0=mean(test1), y0=0, y1=curveheight) gave me

Re: [R] Re : Re : descriptive statistics

2010-12-13 Thread Matthieu Lesnoff
You could also use aggstat() of package tdisplay (available at http://forums.cirad.fr/logiciel-R/viewtopic.php?t=3367). See the help page. mydata - data.frame( + y1 = c(NA, rnorm(n = 8, mean = 10, sd = 5), NA), + y2 = c(rep(NA, 2), rnorm(n = 6, mean = 10, sd = 5), rep(NA, 2)), + y3 =

[R] Post positive reviews

2010-12-13 Thread Eugene Zola
Google’s Huge Change and How it affects you. •Anyone can now post bad reviews and kill your rank. •We post good reviews and improve your rank. •We post good reviews to keep others from killing your rank. Google: Judge, Jury and Online Shopping Executioner Google rank is based on

[R] stepAIC: plot predicted versus observed

2010-12-13 Thread sergeetienneparent
Hi, stepAIC generic plot function creates useful graphics for the diagnosis of multiple regressions. To create predicted versus observed plots, I use to look for the coefficients, copy them by hand, calculate R², then plot. Is there a more automated way to plot predicted versus observed with its

Re: [R] How to ignore data

2010-12-13 Thread Steve Sidney
Oh dear oh dear!!! another arrogant statistician/scientist One asks for help and instead one gets an ear full!!! So much for the much vaunted helpful R community. But thanks anyway, I guess you were trying Steve On 2010/12/13 08:17 PM, Bert Gunter wrote: Inline below. -- Bert On

Re: [R] curve

2010-12-13 Thread Sarah Goslee
On Mon, Dec 13, 2010 at 2:12 PM, Ashta sewa...@gmail.com wrote: Thanks Sarah, 1. to shade or color (blue) the curve using the criterion that any values greater than 11,000 I think I was not clear in the above point. I want shade not the line but the area under the curve, Here's an example

Re: [R] How to ignore data

2010-12-13 Thread Ben Bolker
Steve Sidney sbsidney at mweb.co.za writes: Oh dear oh dear!!! another arrogant statistician/scientist One asks for help and instead one gets an ear full!!! So much for the much vaunted helpful R community. But thanks anyway, I guess you were trying Steve I know I

Re: [R] How to ignore data

2010-12-13 Thread Barry Rowlingson
On Mon, Dec 13, 2010 at 7:36 PM, Steve Sidney sbsid...@mweb.co.za wrote: Oh dear oh dear!!! another arrogant statistician/scientist One asks for help and instead one gets an ear full!!! So much for the much vaunted helpful R community. But thanks anyway, I guess you were trying

[R] Does a formula object have a left hand side

2010-12-13 Thread Erik Iverson
Hello, Does anyone know of a function that will determine whether or not a formula object has a left hand side? I.e., can differentiate between y ~ x + z and ~ x + z Perhaps I'm overlooking the obvious... Thanks! __ R-help@r-project.org mailing

Re: [R] Does a formula object have a left hand side

2010-12-13 Thread William Dunlap
attr(terms(formula), response) is 1 if the formula has a left hand side and 0 otherwise. At a lower level, you can look at length(formula): 2 means there is no LHS, 3 means there is (any other value indicates that someone made a call object that the parser would not make). Bill Dunlap Spotfire,

Re: [R] Does a formula object have a left hand side

2010-12-13 Thread Phil Spector
Erik - Perhaps the response attribute of the terms() function? formula1 = formula(y ~ x + z) formula2 = formula(~x + z) attr(terms(formula1),'response') [1] 1 attr(terms(formula2),'response') [1] 0 Although there may be more direct ways. - Phil

Re: [R] How to ignore data

2010-12-13 Thread Peter Ehlers
On 2010-12-13 11:36, Steve Sidney wrote: Oh dear oh dear!!! another arrogant statistician/scientist One asks for help and instead one gets an ear full!!! So much for the much vaunted helpful R community. But thanks anyway, I guess you were trying Steve Ouch!! I didn't offer advice

Re: [R] Does a formula object have a left hand side

2010-12-13 Thread Erik Iverson
Excellent, thank you all! William Dunlap wrote: attr(terms(formula), response) is 1 if the formula has a left hand side and 0 otherwise. At a lower level, you can look at length(formula): 2 means there is no LHS, 3 means there is (any other value indicates that someone made a call object that

Re: [R] Integration with LaTex and LyX

2010-12-13 Thread Yihui Xie
I tried hard to write an automagic script to configure LyX so that you don't need to go to the instructions on CRAN (http://cran.r-project.org/contrib/extra/lyx/): http://yihui.name/en/2010/10/how-to-start-using-pgfsweave-in-lyx-in-one-minute/ This works for LyX 1.6.x and major OS'es with

[R] Date variable error

2010-12-13 Thread Kurt_Helf
Greetings In attempting to create a date variable based on month (e.g., February, April, etc.) and year (e.g., 2006) data, wherein I converted Month to a factor with Jan=1...Dec=12, I used the following command: data$Date-mdy.date(month=data$Month,day=15,year=data$Year) however, I get a

[R] OFFTOPIC: BAD SCIENCE by Ben Goldacre

2010-12-13 Thread Bert Gunter
Folks: This is off topic, but I believe many R-Help participants would be interested in this. My apologies to my British colleagues, who probably already know about this, and to others for whom this is a waste of their time. Dr. Ben Goldacre, a British Physician and science columnist, has

  1   2   >