search for: michela

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060627/84d9bd54/attachment.pl

My name is Michela Marignani and I'm an ecologist trying to solve a problem linked to knight' s tour algorithm. I need a program to create random matrices with presence/absence (i.e. 1,0 values), with defined colums and rows sums, to create null models for statistical comparison of species distribution phe...

Hi, can you tel me how can I make a Forrest Plot with R? It is possible and easy or are there a more practical free software available? Than you Mic [[alternative HTML version deleted]]

Dear R-Users, I have a model with a latent variable for a spatio-temporal process. I would like to use EM algorithm to estimate the parameters. Does anybody know how to implement the algorithm in R? Thank you very much in advance, Michela

..., see the modifications Diaconis and Sturmfels make (the "Metropolis" step). Reid Huntsinger -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Huntsinger, Reid Sent: Thursday, January 27, 2005 10:50 AM To: 'Michela Marignani'; r-help at stat.math.ethz.ch Subject: RE: [R] Request for help Persi Diaconis and Bernd Sturmfels have an article on generating random contingency tables uniformly distributed subject to having fixed marginals for the same purpose (null distribution of conditional test) and they us...

....000 You can see that using the master it takes 44.168 seconds for computing the function on matrix XX while it takes 71.562 seconds (more time!!!) with the cluster. Can you give us some advice in order to understand why the cluster is slower than the master? Thank you very much in advance, bye Michela and Marco Ps: we have a gigabit ethernet between the master and the nodes

...th x values = 0, 10, 20, 30 while I want a graph with values 0, 6, 12, 18, 24 in the x axis. How can I do this? In R 2.1.1 version there was "time.inc" in survplot, but in version R 2.2.0 there isn't it! I am sorry for my english and I hope that you understand my problem. Thank you Michela ************************************** Dr.ssa Michela Ballardini Unità di Biostatistica e Sperimentazioni Cliniche c/o Osp. Morgagni-Pierantoni - Pad. Valsalva Via Forlanini, 34 47100 Forlì Tel 0543-731836 Tel/Fax 0543-731612 ************************************** [[alternative HTML version del...

...I does not show any problem, e.g.: a<-list.files() for (x in a) { u<-read.csv(x, header=T, sep=",", check.names=FALSE) } h<-u[,2] I would appreciate any suggestion or pointer to solve the problem or to do the same thing in a different way. Thanks for your consideration -- Michela Leonardi

...ime,family=poisson,fit.method="mgcv") H<-predict.gam(mod_spline,type="lpmatrix") H<-H[,2:16] I presume that the vectors obtained using the predict.gam() are different from the ones obtained using GAM.setup(). Thank for your help. Looking forward to hearing from you. Michela [[alternative HTML version deleted]]

...uot; "utils" "datasets" [13] "base" #searching for function "abind" belonging to library "abind" > abind Error: object "abind" not found Is it possible to know which libraries each node has loaded? Thank you very much best regards Michela and Marco

...st ############################## sumMatrices <- function(matrices){ if (length(matrices) > 2) matrices[[1]] + Recall(matrices[-1]) else matrices[[1]] + matrices[[2]] } ############################## result = sumMatrices(e1) * Thank you in advance for all your help, best Michela [[alternative HTML version deleted]]

...a 2 y a 3 z a 4 x b 5 y b 6 x c 7 y c 8 z c >From the above data I'd like to obtain the following matrix: a b c x 1 1 1 y 1 1 1 z 1 0 1 Do you have any advice? Can you help me please? Thank you in advance, Michela

...would like to add a new "day" variable to group data by day, like this: id; words; day 2005-07-07T09:59:56Z; 35; 1 2005-07-07T10:01:39Z; 13; 1 2005-07-08T10:02:22Z; 1; 2 2005-07-09T10:03:16Z; 23; 3 2005-07-10T10:04:23Z; 39; 4 2005-07-10T10:04:39Z; 15; 4 How can I do that? Many thanks! Michela

I'm fitting a Poisson gam model, say model<-gam(a65tm~as.factor(day.week )+as.factor(week)+offset(log(pop65))+s(time,k=10,bs="cr",fx=FALSE,by=NA,m=1),sp=c( 0.001),data=dati1,family=poisson) Currently I've difficulties in obtaining right predictions by using gam.predict function with MGCV package in R version 2.2.1 (see below my syntax).

...itted to IEEE Xplore and for EI indexing. SPECIAL JOURNAL ISSUE: The best papers of Cluster 2015 will be included in a Special Issue on advances in topics related to cluster computing of the Elsevier International Journal of Parallel Computing (PARCO), edited by Pavan Balaji, Satoshi Matsuoka, and Michela Taufer. This special issue is dedicated for the papers accepted in the Cluster 2015 conference. The submission to this special issue is by invitation only. IMPORTANT DATES September 27, 2014 .... Submissions open for Workshops January 1, 2015 ........... Submissions open for Papers, Posters, and T...

...itted to IEEE Xplore and for EI indexing. SPECIAL JOURNAL ISSUE: The best papers of Cluster 2015 will be included in a Special Issue on advances in topics related to cluster computing of the Elsevier International Journal of Parallel Computing (PARCO), edited by Pavan Balaji, Satoshi Matsuoka, and Michela Taufer. This special issue is dedicated for the papers accepted in the Cluster 2015 conference. The submission to this special issue is by invitation only. IMPORTANT DATES September 27, 2014 .... Submissions open for Workshops January 1, 2015 ........... Submissions open for Papers, Posters, and T...

...nd at each time time the observed variable y has dimension equal to d=3. Moreover I want to use white noise errors. But with the command SS(y) where y is my nXd data matrix (d=observations at each time, n=numbers of time points) I obtain a different model. Can you help me please? Thanks in advance Michela

...ch for OpenMPI but we are wondering if in the meanwhile the package has been compiled also for OpenMPI. Does anybody use RScalapack and has any suggestions? Sorry for writing here but nobody answered in the R-sig- hpc (High performance computing with R) mailing list. Thank you very much. My best, Michela [[alternative HTML version deleted]]

Hello, someone knows which function has replaced the GamSetup function in mgcv? With R 1.6 version we have GamSetup but this function doesn't exist with R 2.7.2 version [[alternative HTML version deleted]]

...y but when we try to install the package we get ** testing if installed package can be loaded Error in library.dynam(lib, package, package.lib) : shared library 'mylib' not found we cannot understand what's wrong and where R searches for the mylib.so files. Any suggestion? Thanks, Michela and Francesco