Re: [R] Binary data on JAGS

2021-06-01 Thread Massimiliano Tripoli
Thanks, for suggestions (To define p before). I thought it and I put this line of code inside the loop: p[i] <- y[i] / n[i] but I received a RUNTIME ERROR and: (Attempt to redefine node p[1]). Thanks Massimiliano Il giorno mar 1 giu 2021 alle ore 17:22 Bert Gunter ha scritto: > Wher

[R] Binary data on JAGS

2021-06-01 Thread Massimiliano Tripoli
Dear R users, I'm trying to reproduce the example 6.5.1 (Dobson (1983)) in BUGS book in linux using JAGS. Below the code as https://www.mrc-bsu.cam.ac.uk/software/bugs/the-bugs-project-the-bugs-book/bugs-book-examples/the-bugs-book-examples-chapter-6-6-5-1/ # By R code: library('rjags') jags <

[R] Wilcoxon Test

2017-04-21 Thread Tripoli Massimiliano
Dear R users, Why the result of Wilcoxon sum rank test by R is different from sas https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_npar1way_sect022.htm The code is next: sampleA <- c(1.94, 1.94, 2.92, 2.92, 2.92, 2.92, 3.27, 3.27, 3.27, 3.27, 3.7

[R] Freq in dataframe

2016-08-05 Thread Tripoli Massimiliano
) Dose 1.A Dose 2.A Dose 3.A Dose 1.B Dose 2.B Dose 3.B Dose 1.C Dose 2.C Dose 3.C Week 32878767 117 8 Week 52 10 13 1599 10 13 16 16 Someone could help me ? T

[R] par function for barchart

2015-01-14 Thread Massimiliano Tripoli
;, "X1", "X2"), class = "data.frame", row.names = c(NA, -19L)) require(lattice) op <- par(mfrow=c(1,2)) barchart(REGIONE ~ X1,data=mydata,col="red") barchart(REGIONE ~ X2,data=mydata,col="green") par(op) Can anyone help me, please? Thanks,

Re: [R] overlapping coefficient bidimensional distribution

2015-01-13 Thread Meli Massimiliano
Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, Jan 12, 2015 at 5:09 AM, Meli Massimiliano wrote: Hello All, I have two set of data (A and B) a

[R] overlapping coefficient bidimensional distribution

2015-01-12 Thread Meli Massimiliano
Hello All, I have two set of data (A and B) and i would calculate the percentage of superimposition of these two set. Or the number of point that belong to the intersection A ∩ B of these two sets. How i can measure it ? thanks a lot, max __ R-help

Re: [R] Converting list to character

2014-11-26 Thread Massimiliano Tripoli
Thanks David, that's I was looking for. Thanks to Chel too. Massimiliano - Messaggio originale - Da: "David L Carlson" A: "Chel Hee Lee" , "Massimiliano Tripoli" , r-help@r-project.org Inviato: Martedì, 25 novembre 2014 19:40:51 Oggetto: RE: [R] Co

[R] Converting list to character

2014-11-25 Thread Massimiliano Tripoli
chr "GR.3.1" "GR.3.8" ..$ 9 : chr "GR.3.8" ..$ 10: chr "GR.3.7" ..$ 11: chr "GR.3.1" ..$ 12: chr "GR.3.8" ..$ 13: chr "GR.3.1" "GR.3.8" ..$ 15: chr "GR.3.8" ..$ 16: chr "GR.3.1" ..

Re: [R] Replacing NA's in the data set

2014-04-29 Thread Massimiliano Tripoli
Hi, Try this: for (i in 1:10) Ts[is.na(Ts[,i]),i] <- Ts$V11[is.na(Ts[,i])] Massimiliano - Messaggio originale - Da: "dila radi" A: r-help@r-project.org Inviato: Martedì, 29 aprile 2014 10:46:32 Oggetto: [R] Replacing NA's in the data set Hi all, I have this d

Re: [R] problem of tapply function

2014-03-30 Thread Massimiliano Tripoli
Thanks Bert - Messaggio originale - Da: "Bert Gunter" A: "Massimiliano Tripoli" Cc: r-help@r-project.org Inviato: Sabato, 29 marzo 2014 15:52:59 Oggetto: Re: [R] problem of tapply function The argument to tapply is FUN, not fun. R is case sensitive. -- Bert B

[R] problem of tapply function

2014-03-29 Thread Massimiliano Tripoli
1 7 9 5 5 2 4 6 3 16146 16147 16150 16151 16153 16154 16155 16157 16158 5 3 1 2 1 1 4 1 1 > # By looking tapply function help why the table's output is different from tapply's output? # Anyone coul

[R] Problem of "!" operator

2014-03-28 Thread Massimiliano Tripoli
haracter(esn) Error: object 'esn' not found > for (i in 2:length(esn)) if(esn[i]==esn[i-1]) + idi[i] <- idi[i-1] else idi[i] <- idi[i-1] + 1 Error: object 'esn' not found > idi Error: object 'idi' not found > } Error: unexpected '}' i

Re: [R] sub function problem

2014-03-03 Thread Massimiliano Tripoli
;PIKAL Petr" A: "Massimiliano Tripoli" , r-help@r-project.org Inviato: Lunedì, 3 marzo 2014 10:36:02 Oggetto: RE: [R] sub function problem Hi Everything works for me as expected. Nothing happens R starts as usual. You can shorten your function corregge2 <- function(stringa){ nc <

[R] sub function problem

2014-02-28 Thread Massimiliano Tripoli
fhsfh") [1] 8 > sub("^ ", "", " fghaghf") [1] "fghaghf" > And restarting R I had this message: Error:10061 R is not in server mode. Please, start R and/or run the startSocket Server function available in svSocket package! Anyone could hel

[R] fitting

2012-11-16 Thread Meli Massimiliano
Hello All, i would fit my data with a function like this : y = a0 + a1 * exp(-x/a2) + a3 * exp(-x/a4) + a5 * exp(-x/a6) + a7 * exp(-x/a8) + a9 * exp(-x/a10) plus i have to impose that a1 + a3 + a5 + a7 + a9 = 1 a1 , a3 , a5 , a7 , a9 >= 0 i try with nls, but i do not know how to impose the

[R] R graphs differ from exported one

2011-06-12 Thread Massimiliano
Hello everybody! This is my first mail so I'll write a couple of lines of self-introduction. My name is Massimiliano, I'm from Italy and I'm studying Mathematical Engineering. I started using R in my Statistics course and have to use it to make a project which I'll discuss at t

[R] paste function by INDICES

2011-04-27 Thread Massimiliano Tripoli
ccording to ID variable. For example for ID=2 "b" should be 13, for ID=3 "b" should be 114, and for ID=4 should be 1241. How can I do that ? Thank you for any advices, Regards, Massimiliano Tripoli [[alternative HTML version deleted]]

Re: [R] Include C++ DLL, error in ...: C symbol name not in load table

2011-04-21 Thread Massimiliano
Maybe you have to add #include #include to your C++ source code. Massimiliano Tripoli Il giorno mer, 20/04/2011 alle 22.06 +0200, Sascha Vieweg ha scritto: > #include > using namespace std; > void foo (double* x, double* y, double* out) > { > out[