[R] Mediation using 'mediate' in R with an interaction

Hi everyone, I have a question about mediation when there is an interaction between 2 categorical variables, one with two levels and one with three levels. Using the ‘mediation’ package by Tingley and colleagues in R, I’m trying to conduct a bootstrapped mediation analysis for an interacti

Re: [R] Lattice Histogram with Normal Curve - Y axis as percentages

Just an afterthought if any one really needs to do it again. The crux of the matter is the different limits in the 2 plot x and y scales. It may be easier to accomplish this with a prepanel function to get the limits eg panel.loess Duncan -Original Message- From: r-help-boun...@r-proj

Re: [R] adding rows

Dear Rui and Murphy, Thanks for your help. Eliza > Date: Fri, 9 May 2014 22:55:27 +0100 > From: ruipbarra...@sapo.pt > To: eliza_bo...@hotmail.com; r-help@r-project.org > Subject: Re: [R] adding rows > > Hello, > > Try the following. > > sapply(1:(30 - 2), function(i) sum(el[i:(i+2), ])) > >

Re: [R] adding rows

Hello, Try the following. sapply(1:(30 - 2), function(i) sum(el[i:(i+2), ])) but with number of rows instead of 30. Hope this helps, Rui Barradas Em 09-05-2014 22:35, eliza botto escreveu: Dear useRs, I have a matrix, say "el" of 30 rows and 10 columns, as el<-matrix(sample(1:300),ncol=10)

Re: [R] as.character(quote(x$y) ) = "$" "x" "y" not "x$y"?

Dear Hadley: Thanks for that. Digits are not numbers. Numbers are not data. Data is not information. Information is not intelligence. Intelligence is not knowledge. Knowledge is not wisdom. And your "is." warnings are more useful than my trivia here. Spencer On 5/9/2014

Re: [R] transition from depends to imports

Duncan Murdoch writes: > On 09/05/2014, 5:39 PM, Andreas Leha wrote: >> Hi Duncan, >> >> Thank you for your follow-up and fast response! >> >> Duncan Murdoch writes: >> >>> On 09/05/2014, 5:08 PM, Andreas Leha wrote: Hi all, I am having trouble to do a particular transition from d

Re: [R] as.character(quote(x$y) ) = "$" "x" "y" not "x$y"?

Beware of the is.* functions: * is.object() does not test the usual definition of objects * is.vector() does not test the usual definition of vectors * is.numeric() does not work the same way as is.character() or is.integer() * is.Date() doesn't exist * is.nan() doesn't return TRUE for some NaNs

Re: [R] transition from depends to imports

On 09/05/2014, 5:39 PM, Andreas Leha wrote: Hi Duncan, Thank you for your follow-up and fast response! Duncan Murdoch writes: On 09/05/2014, 5:08 PM, Andreas Leha wrote: Hi all, I am having trouble to do a particular transition from depends to imports in one of my packages. This packages

Re: [R] transition from depends to imports

Hi Duncan, Thank you for your follow-up and fast response! Duncan Murdoch writes: > On 09/05/2014, 5:08 PM, Andreas Leha wrote: >> Hi all, >> >> I am having trouble to do a particular transition from depends to >> imports in one of my packages. >> >> This packages uses 'wilcoxsign_test' from th

[R] adding rows

Dear useRs, I have a matrix, say "el" of 30 rows and 10 columns, as el<-matrix(sample(1:300),ncol=10) I want to sum up various sets of three rows of each column in the following manner sum(el[c(1,2,3),]) ##adding row number 1, 2 and 3 of each column sum(el[c(2,3,4),])##adding row number 2, 3 and 4

Re: [R] as.character(quote(x$y) ) = "$" "x" "y" not "x$y"?

Hi, Duncan: Thanks very much. I used to think that everything in R was a object. Now I know that is.object(quote(x)) is FALSE. (A decade ago, S-Plus asked me if I wanted to save changes to history. I thought, "Wow! Do I get to change history?" Hadley's "Advanced R" book me

Re: [R] transition from depends to imports

On 09/05/2014, 5:08 PM, Andreas Leha wrote: Hi all, I am having trouble to do a particular transition from depends to imports in one of my packages. This packages uses 'wilcoxsign_test' from the 'coin' package. And this is the only function from the coin package that it uses (directly). (Ever

Re: [R] Help with passing in name of XTS file into write.table as a variable.

It is a simple change. If you don't save the results of getSymbol() into a variable, then it's not available for write.table() to write. Try mydata <- getSymbols(ticker,from='1990-01-01') write.table( mydata, {etc}) The help for write.table says that "x" (the first argument to write.table),

[R] transition from depends to imports

Hi all, I am having trouble to do a particular transition from depends to imports in one of my packages. This packages uses 'wilcoxsign_test' from the 'coin' package. And this is the only function from the coin package that it uses (directly). (Everything works fine, as long as my package depen

Re: [R] How to write a loop in which i has to be in POSIX format?

i is an element in zn so replace zn[i] with just i for (i in zn){ treat$su[treat$Zugnacht==as.POSIXct(i, "UTC")] <- min(treat$Vollzeit[treat$Zugnacht==as.POSIXct(zni, "UTC")]) treat$sa[treat$Zugnacht==as.POSIXct(zni, "UTC")] <- max(treat$Vollzeit[treat$Zugnacht==as.POSIXct(zni, "UTC")])

Re: [R] Pearson III distribution

On May 9, 2014, at 8:17 AM, NOELIA LEGAL wrote: > Please I need your help.I'm intrested to know if there is any R-package for > fit a Pearson III extreme value distribution to data.Thanks a lot You are requested to do some searching yourself before posting. If you have failed to find what yo

Re: [R] Error during working with wgcna and R

WGCNA maintainer here. When working with a large data set, you have a few options. 1. Without being snarky, the best option is to get (or get access to) a computer with large-enough RAM. Many universities, departments, and other research institutes have computer clusters with nodes with at least 6

[R] Pearson III distribution

Please I need your help.I'm intrested to know if there is any R-package for fit a Pearson III extreme value distribution to data.Thanks a lot Noelia [[alternative HTML version deleted]] __ R-help@r-project.o

Re: [R] Error during working with wgcna and R

I am also working on co-expression analysis. It seems like there is no way to use TOMsimilarityFromExpr for large datasets. The option 'maxBlockSize' exists for module detection but not for TOMsimilarity? The only solution seems to reduce the dataset On Sunday, July 8, 2012 1:02:54 PM UTC+2, de

Re: [R] Tukey "no factors"

HI, May be the below links help: http://r.789695.n4.nabble.com/TukeyHSD-troubles-td1570205.html http://r.789695.n4.nabble.com/TukeyHSD-error-td3051059.html A.K. Hi I'm trying to run a Tukey test on mortality data, where I want to test whether mortality is influenced by the amount of copper (

[R] Resquest

Dear Sir, I need Your help,I am working on the comparison of l moments and trimmed l moments .,Now i face a problem in R programming of trimmed l moments. for this purpose i need your help .for example let suppose i have the following data vector abchttps://stat.ethz.ch/mailman/listinfo/r-help PL

Re: [R] DESeq vs DESeq2 different DEGs results

You're on the wrong list. This is more appropriate on the bioconductor mailing list. On Mon, May 5, 2014 at 9:42 AM, Catalina Aguilar Hurtado wrote: > Hi, > > I want to compare DESeq vs DESeq2 and I am getting different number of DEGs > which I will expect to be normal. However, when I compare t

Re: [R] How to examine the parameter values during iteration (nls)

Read ?nls and note the "trace" argument. Does nls(...,trace = TRUE) not do what you want? By using ?capture.output you could then capture the trace I would think, but I haven't tried. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Inform

Re: [R] as.character(quote(x$y) ) = "$" "x" "y" not "x$y"?

Ahhh. Thanks Duncan. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." H. Gilbert Welch On Fri, May 9, 2014 at 2:41 AM, Duncan Murdoch wrote: > On 09/05/2014, 2:41 AM, Bert G

[R] How to examine the parameter values during iteration (nls)

Hi everyone, I am using nls() to run a simple sigmoid curve type of regression and would like to see how the parameter values change through iterations. How can I see those values? Or ideally, can I even extract those values? For example, nls(Response~E0+Emax*Conc/(EC50+Conc), data=data, start=

Re: [R] How can I make this nested loop faster?

Your code is not re-producable. Can you provide a working example using a standard dataset from R? But, you could first try to use compiler package, see ?enableJIT. Another option would be to use doMC/foreach packages if you can run your assignment in the nested loop in parallel, see %dopar%. On

Re: [R] problem in r-code

Wrong list. This is an R list not Bugs. You may want to consult Bugs materials: http://www2.mrc-bsu.cam.ac.uk/bugs/weblinks/webresource.shtml On 8 May 2014 11:36, thanoon younis wrote: > dear all members > > is there anyone explain to me the code below and how can i transfer this > code to winbug

Re: [R] as.character(quote(x$y) ) = "$" "x" "y" not "x$y"?

On 09/05/2014, 2:41 AM, Bert Gunter wrote: Spencer: Hmmm Well, I don't get what's going on here -- as.character.default is internal -- but could you method-ize a simple synonym: See ?InternalMethods: "For efficiency, internal dispatch only occurs on objects, that is those for which is.ob

Re: [R] group by and merge two dataframes

yes thanks, that's correct! here a slight variation inspired by your solution: a cartesian product restricted to non duplicated records to get the logical vector i to be used in the next natural join i<-!duplicated(merge(df1$id,df1$item, by=NULL)) merge(df1[i,],df2) thanks Il 08/05/2014 18: