Re: [R] x-axis tick marks on log scale plot

Very, very nice. Thanks for sharing. On 5/20/2016 4:21 AM, Martin Maechler wrote: >> Brian Smith >> on Thu, 19 May 2016 11:04:55 -0400 writes: > > Thanks all !! On Thu, May 19, 2016 at 9:55 AM, Ivan > > Calandra

Re: [R] Change sum of squares type for ANOVA

Dear Michael and Michael, I took a quick look at the Statistica website, and their definition of "types" of sums of squares doesn't appear to be the same as the common definition, which originated with SAS. As I read their materials, what they call "Type VI" sums of squares would correspond to

Re: [R] Change sum of squares type for ANOVA

Dear Michael F What you say is of course true but Michael E (or possibly Simone) seems to be using a function glht from an unnamed package. Since we do not have the data or the offending output it is hard to be sure though. Michael D On 20/05/2016 13:49, Michael Friendly wrote: Use

[R] Change sum of squares type for ANOVA

Dear R-list members, I compared my statistics with my supervisor yesterday. He is using STATISTICA, I am using R. We both loaded the same data-file and did a two-way ANOVA with treatment and trial as factors; treatment means were then compared to the untreated control using Dunnett's test.

[R] rjags crashes RStudio every time it tries to load

Hello, Every time I try to require(rjags) in Studio, Studio crashes and says that a fatal error has occurred. This seems to happen when it's trying to require coda, which I have installed. All of my software is updated. Anyone else encountered this before? The package does seem to load

Re: [R] x-axis tick marks on log scale plot

> Brian Smith > on Thu, 19 May 2016 11:04:55 -0400 writes: > Thanks all !! On Thu, May 19, 2016 at 9:55 AM, Ivan > Calandra wrote: >> Hi, >> >> You can do it by first plotting your values without the >>

[R] Multistate survival models

A quick (I think) question. Does the survreg() function in the survival package handle Surv() objects of type mstate? It's not clear from the documentation that it doesn't. The Vignette only discusses the Cox model and the Fine-Gray model. I tried using survreg() on an mstate Surv()

Re: [R] rjags crashes RStudio every time it tries to load

On 20/05/2016 2:28 PM, Ryan Utz wrote: Hello, Every time I try to require(rjags) in Studio, Studio crashes and says that a fatal error has occurred. This seems to happen when it's trying to require coda, which I have installed. All of my software is updated. Anyone else encountered this before?

Re: [R] Plots for lmrob function

> On May 20, 2016, at 8:28 AM, varin sacha via R-help > wrote: > > Dear R-helpers, > > I have fitted a robust regression using lmrob function from robustbase > package. I try to get the different plots for diagnostics of residuals and > others. I can't get them, a

Re: [R] stopifnot() doesnt work as I expect it to. Are my expectations correct?

On 20/05/2016 4:44 AM, Vasanth Mohan wrote: Hi, *stopifnot(FALSE, someOtherExpression)* For the above code I expect stopifnot() to always say that 'FALSE is not TRUE' regardless of what someOtherExpression is(It may evaluate to TRUE or FALSE or throw an error). Is my expectation correct? Is

Re: [R] stopifnot() doesnt work as I expect it to. Are my expectations correct?

The following usage of stopifnot seems reasonable to me and it would be nice if the 2nd call caused the message 'is.data.frame(df) is not TRUE'. f <- function(df) { stopifnot(is.data.frame(df), is.integer(df$ID)) range(df$ID) } f(data.frame(ID=4:7)) # [1] 4 7 f(4:7) # Error in df$ID : $

Re: [R] Choosing subset of data

> On May 20, 2016, at 11:46 AM, oslo via R-help wrote: > > Hi all; > I have a big data set (a small part is given below) and V1 column has > repeated info in it. That is rs941873, rs12307687... are repeating many > times. I need choose only one SNP (in first column named

Re: [R] (no subject)

> On May 20, 2016, at 11:43 AM, Sibasish Saha wrote: > > What is the code in R to get the goodness of fit in nls?? Can you "back up" a bit and present the _theory_ you have in mind that lets you assess the GOF? > Please sir... > please help me You should also

[R] (no subject)

What is the code in R to get the goodness of fit in nls?? Please sir... please help me [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Choosing subset of data

Hi all; I have a big data set (a small part is given below) and V1 column has repeated info in it. That is rs941873, rs12307687... are repeating many times. I need choose only one SNP (in first column named rs) which has the smallest  Pvalue withing V1 column. That is I need choose only one SNP

[R] Inverse normal transformation of ranked data

I am using ddply on a data set that contains 2+ million rows; trying to rank the values of a variable within groups, and then transform the ranks to (approximate) z-scores --- i.e generate quantiles on the normal scale. Here is some sample data for one group:x <- NA 0.3640951 0.1175880 0.3453916

Re: [R] stopifnot() doesnt work as I expect it to. Are my expectations correct?

Dear Bill, assertthat has such capabilities f <- function(df){ require(assertthat) assert_that(is.data.frame(df)) assert_that(is.integer(df$ID)) range(df$ID) } f(data.frame(ID=4:7)) # [1] 4 7 f(4:7) # Error: df is not a data frame f(data.frame(ID=letters)) # Error: df$ID is not an

Re: [R] stopifnot() doesnt work as I expect it to. Are my expectations correct?

On 20/05/2016 10:33 AM, William Dunlap wrote: The following usage of stopifnot seems reasonable to me and it would be nice if the 2nd call caused the message 'is.data.frame(df) is not TRUE'. f <- function(df) { stopifnot(is.data.frame(df), is.integer(df$ID)) range(df$ID) } Yes, that's

[R] Instruments for endogenous interaction term in Simultaneous Equation Model

Hi I want to estimate the following SEM model: Y1 = a0 + b0*D + c0* Y3 + d0 * Y3 * D + control variables Y2 = a1 + b1*D + c1 *Y3 + d1 * Y3 * D + control variables Y3 = a2 + b2*D + c2 *Y1 + d2 * Y1 * D + e2*Y2 + f2 * Y2* D + control variables where, D is dummy variable, Y1, Y2 and Y3 are

[R] Plots for lmrob function

Dear R-helpers, I have fitted a robust regression using lmrob function from robustbase package. I try to get the different plots for diagnostics of residuals and others. I can't get them, a window opens but nothing appears on it (the window remains white, no graph appears) and I get this error

[R] stopifnot() doesnt work as I expect it to. Are my expectations correct?

Hi, *stopifnot(FALSE, someOtherExpression)* For the above code I expect stopifnot() to always say that 'FALSE is not TRUE' regardless of what someOtherExpression is(It may evaluate to TRUE or FALSE or throw an error). Is my expectation correct? Is that how stopifnot() is supposed to work? The

Re: [R] error in data.farme--duplicate row.names error

Don, Thank you for your helpful response. At this point, I do believe it is a package error and I have contacted the developer. Thanks, Lisa -Original Message- From: MacQueen, Don [mailto:macque...@llnl.gov] Sent: Thursday, May 19, 2016 4:27 PM To: Rees, Lisa Marie (MU-Student);

Re: [R] how to perform multiple comparison？

thanks for your help. 发自 网易邮箱大师 On 2016-05-20 09:48 , David Winsemius Wrote: > On May 19, 2016, at 5:19 PM, Jim Lemon wrote: > > Hi laomeng_3, > Have a look at the padjust function (stats). > > Jim > > > On Fri, May 20, 2016 at 1:56 AM, laomeng_3

Re: [R] Change sum of squares type for ANOVA

Use car::Anova() for type II and type III sums of squares and F-tests. The sequential type I tests computed by anova() are rarely sensible On 5/20/2016 3:11 AM, michael.eisenr...@agroscope.admin.ch wrote: Dear R-list members, I compared my statistics with my supervisor yesterday. He is using

Re: [R] Linear Regressions with non-negativity constraint

Hi all, I hope you are doing well? I'm currently using lm() to estimate a linear multi-factor (5 factors without intercept) model as follows ... factor.lm <- lm(y~x1+x2+x3+x4+x5-1, data = data.frame.rbind) Using nnls(A,b) I estimated the same model, extended by a non-negativity constraint on

Re: [R] how to perform multiple comparison？

ok,many thanks. 发自 网易邮箱大师 On 2016-05-20 08:19 , Jim Lemon Wrote: Hi laomeng_3, Have a look at the padjust function (stats). Jim On Fri, May 20, 2016 at 1:56 AM, laomeng_3 wrote: > Hi all: > As to the anova, we can perform multiple comparison via TukeyHSD. > But as to