[R] rpart.object help

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

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

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 m

Re: [R] descriptive statistics

?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 wi

Re: [R] package sampling

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

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

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 t

Re: [R] descriptive statistics

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 couldn

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

R-11.1 on both 32-bit and 64-bit on Windows with the windows() device. Best Marcin 2010/12/13 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,

[R] Re : descriptive statistics

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 À : effeesse Cc : r-help@r-project.org Envoyé le : Lun 13 décembre 2010, 11h 23min 15s Objet : Re: [R] descriptive

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

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 once

[R] How to compare two square matrices

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

Another way is the remix function of the remix package. On Monday, December 13, 2010, justin bem 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 > > > > > > De : Ji

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

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 func

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

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

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 'xpinc

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

> 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 wrote: > > On 2010-12-11 03:12, Francesco N

Re: [R] rpart.object help

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,

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

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 (don'

[R] post-hoc test for ANCOVA method

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()

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 lm.fit()

Re: [R] simple plotting question

to add to Michael's response: http://www.statmethods.net/advgraphs/parameters.html On Mon, Dec 13, 2010 at 2:23 AM, Michael Bedward wrote: > Hello Erin, > > Try this... > > plot(x, y, type="b", pch=16) > > Michael > > On 13 December 2010 18:11, Erin Hodgess wrote: >> Dear R People: >> >> When

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

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

Re: [R] Need help on nnet

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 > -4

Re: [R] Help..Neural Network

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 do

[R] help on SAS Macro in R

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

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)

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)? Se

Re: [R] help on SAS Macro in R

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 ev

Re: [R] help on SAS Macro in R

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 world

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

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 go

Re: [R] Re : descriptive statistics

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 se

[R] How does R compute sums of squares?

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

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

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

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, Kattankulathur

[R] Re : Re : descriptive statistics

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 variabl

[R] check for item in vector

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 i

Re: [R] check for item in vector

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 i

[R] Reg : null values in kmeans

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: http://r.789695.n4

Re: [R] check for item in vector

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

Re: [R] check for item in vector

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." To: R-help Date: 12/13/2010 08:50 AM Subject: [R] check for item in vector Sent by: r-help-boun...

[R] booklet on using R for biomedical statistics

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" http://a-little-book-of-r-for-biomedi

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

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 batchfil

[R] How to plot Ellipsoid like function

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 usi

Re: [R] Rapache on windows

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 8:3

Re: [R] How to plot Ellipsoid like function

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

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

[R] Integration with LaTex and LyX

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

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("TukeyHS

Re: [R] ggplot2 histograms

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 of

[R] How to ignore data

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 - ze

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

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 in

Re: [R] Integration with LaTex and LyX

?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 ro

[R] odot symbol as a subscript in axes labels

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 http://www.R

Re: [R] Re : Re : descriptive statistics

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, fun=m

Re: [R] Integration with LaTex and LyX

On Mon, Dec 13, 2010 at 4:55 PM, Jonathan P Daily 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 that's not e

Re: [R] Tukey HSD not working

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

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

On Mon, Dec 13, 2010 at 8:20 AM, Ethan Arenson 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 >          Df  Sum Sq Mean

Re: [R] How to ignore data

Steve Sidney 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 this

[R] Data permutation

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"

Can anybody tell me why this is happening? > library( Defaults) > ls() [1] "class" "foo""mlct" [4] "mlctTheme" "overlayFunction""pct" [7] "pri""primaryFraction""reclassFractions" [10] "reclassMatrix"

Re: [R] Re : Re : descriptive statistics

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 st

Re: [R] Re : Re : descriptive statistics

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 for

[R] 15' lag of an irregular Time Series

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 2010-11-

Re: [R] How to ignore data

>> >> 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. >>

Re: [R] Data permutation

On Mon, Dec 13, 2010 at 11:37 AM, matteop 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       Mark

[R] Qs re writing/reading arrays, dataframes

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 (which

Re: [R] Rapache on windows

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 wrote: > 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 Me

Re: [R] How to ignore data

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 c

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

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-help@r-pro

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

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

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 = KdnM*SpAct*Vo

Re: [R] How to ignore data

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?

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 (En

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

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 o

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

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, bu

Re: [R] How to ignore data

Inline below. -- Bert On Mon, Dec 13, 2010 at 9:42 AM, Steve Sidney 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 or an experiment

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

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

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 cluebrick

[R] curve

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 diffic

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

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 o

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

- --- 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 sur

Re: [R] curve

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$

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

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, Appl

Re: [R] curve

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 th

Re: [R] Re : Re : descriptive statistics

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

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 rev

[R] stepAIC: plot predicted versus observed

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

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 Mon

Re: [R] curve

On Mon, Dec 13, 2010 at 2:12 PM, Ashta 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 of how to

Re: [R] How to ignore data

Steve Sidney 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 shoul

Re: [R] How to ignore data

On Mon, Dec 13, 2010 at 7:36 PM, 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, we're

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

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 l

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

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"

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. - Phi

Re: [R] How to ignore data

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 e

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

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

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 probabi

[R] Date variable error

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 mes

[R] OFFTOPIC: BAD SCIENCE by Ben Goldacre

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 written

  1   2   >