[R] how to read association of variable in multiple outcomes using MCMCMGLMM

2018-11-15 Thread Victor Chikwapulo
Dear all, I am  using package MCMCglmm and I would like to request for an assistant on  what to look   in the output which can tell me  whether there is significant  association among the three antibody titers(logiga,logigm,logigg)  for example association between antibody titers and exposure I

Re: [R] Corrupting files while copying (was Re: saveRDS() and readRDS() Why? [solved, pretty much anyway])

2018-11-15 Thread Ista Zahn
Hi Patrick, I think it would help to start from the beginning and give complete (but concise!) replication instructions, including telling us what host and gest operating systems you are using (including the versions), the version of virtualbox you used, and exactly what steps are needed to

Re: [R] glmutli package assistance please

2018-11-15 Thread Bert Gunter
OK. Then post here but *not* on mixed models list. One or the other, exclusive or. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Nov 15, 2018 at

[R] create a heatmap for findAssocs results based on time

2018-11-15 Thread Elahe chalabi via R-help
Hi all, I have the following data for which I create a document term matrix first and then I add the time available to the dtm. In order to see the correlations to the term "updat" in the different years, I would like to have a heat-map for findassoc in a way that x-axis shows the time.

Re: [R] How to create gridded data

2018-11-15 Thread David L Carlson
It would depend on the format of the gridded data. Assuming it is a data frame like DF2 in my earlier answer, you just reverse the steps: > DF2 110.5 111 111.5 112 46 6.1 4.5 7.8 5.5 45.5 3.2 5.0 1.8 2.0 > DF3 <- data.frame(as.table(as.matrix(DF2))) Var1 Var2 Freq 1 46 110.5

Re: [R] glmutli package assistance please

2018-11-15 Thread Michael Dewey
Dear Bert Since glmulti operates on glm/lm models I think, although I agree about not cross-posting, that it was OK here. Perhaps I do not understand the full significance of mixed models though. Michael On 15/11/2018 15:43, Bert Gunter wrote: Please do not cross post (see te posting

Re: [R] very slow "memoise"

2018-11-15 Thread Jeff Newmiller
Yes, you have misunderstood what memoise is for. First, it is for when you call your function with the same inputs frequently as part of your calling-level algorithm. For your iterative calculation you would have a stuck (cycling) process if the same value of current were to be revisited...

Re: [R] glmutli package assistance please

2018-11-15 Thread Bert Gunter
Please do not cross post (see te posting guide). This should go only to the mixed models list. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Nov

Re: [R] help with grouping data and calculating the means

2018-11-15 Thread Bert Gunter
On Thu, Nov 15, 2018 at 10:40 AM Boris Steipe wrote: > > Use round() with the appropriate "digits" argument. Then use unique() to > define your groups. No. > round(c(.124,.126),2) [1] 0.12 0.13 As I understand it, the OP said he wanted the last decimal to be ignored. The OP also did not

[R] glmutli package assistance please

2018-11-15 Thread Bill Poling
Hello Michael, thank you for your responses. I will have to look through the myriad of iterations of attempts at all this on this project but it seems to me that I tried that? Or tried to convert the integers to numeric and it errored out? If not I certainly will do so and let know, thank you

Re: [R] help with grouping data and calculating the means

2018-11-15 Thread Boris Steipe
Use round() with the appropriate "digits" argument. Then use unique() to define your groups. HTH, B. > On 2018-11-15, at 11:48, sasa kosanic wrote: > > Dear All, > > I would very much appreciate the help with following: > I need to calculate the mean of different lat/long points that

[R] glmutli package assistance please - error Solved

2018-11-15 Thread Bill Poling
Hello Michael and all. 1. I replaced INT with Numeric that did not solve the error problem. 2. However I discovered that my matrix was the problem, I had it set to 2,5 and when I set it to 6,6 I no longer get the error, which was due to not having the correct or more than necessary fields in

Re: [R] help with grouping data and calculating the means

2018-11-15 Thread Bert Gunter
On further thought -- and subject to my prior interpretation -- I think a foolproof way of truncating to 4 decimal digits is to treat them as character strings rather than numerics and use regex operations: > with(df,tapply(TK.QUADRANT, >

Re: [R] help with grouping data and calculating the means

2018-11-15 Thread Anthoni, Peter (IMK)
Hi Sasa, Those latitude look equidistant with a separation of 0.05. I guess you want to calculate the zonal mean along the latitude, right? #estimate the lower and upper latitude for the cut lat.dist=0.05 #equidistant spacing along latitude lat.min=min(df$LAT,na.rm=T)-lat.dist/2

[R] glmutli package assistance please

2018-11-15 Thread Bill Poling
Hi, I have removed the pdf which was causing my e-mail to be blocked by moderators, my apologies. https://www.jstatsoft.org/article/view/v034i12/v34i12.pdf Original post: Hello. I am still trying to get some of the examples in this glmulti pdf to work with my data. I have sent e-mails to

Re: [R] glmutli package assistance please

2018-11-15 Thread Michael Dewey
Dear Bill I am not sure what is going on here but I notice that 2 of your covariates are numeric and 3 integer. What happens if you make them all numeric? Michael On 15/11/2018 11:46, Bill Poling wrote: Hi, I have removed the pdf which was causing my e-mail to be blocked by moderators, my

[R] very slow "memoise"

2018-11-15 Thread Martin Møller Skarbiniks Pedersen
Hi, I want to compute a lot of values and I have tried to use memoise::memoise to speed-up the computation. However it is much slower using the memoised version. I guess I have misunderstood how to use the package memoise or the purpose of the package. The code takes more than 2 minutes to

Re: [R-es] 55 Reunión Grupo Usuarios de R de Madrid...

2018-11-15 Thread Carlos Ortega
Buenas a todos, Ya está disponible el video y las referencias de la reunión que celebramos ayer: http://madrid.r-es.org/55-miercoles-15-de-noviembre-2018/ (Lamentablemente no pudimos grabar la primera charla). Gracias, Carlos Ortega www.qualityexcellence.es El lun., 12 nov. 2018 a las 11:47,

Re: [R-es] package scorecard

2018-11-15 Thread Carlos Ortega
Hola Dayana, No es realmente un "error" es un aviso ("warning"). Te está indicando que en el proceso para transformar a binario tu variable target, una de las variables de tu conjunto, como es constante, no la va a usar. Saludos, Carlos Ortega www.qualityexcellence.es El mié., 14 nov. 2018 a