similar to: Request for help

I referred in my reply to a paper by Diaconis and Sturmfels. The exact reference is: Diaconis and Sturmfels, Algebraic algorithms for sampling from conditional distributions, Ann. Stat 26 (1998) 363-397. They cite the following: Besag and Clifford, Generalized Monte Carlo significance tests, Biometrika 76 (1989) 633-42. which actually contains your problem (section 3, Testing the Rasch model)

----- Messaggio originale ----- > Da: "Rowland Penny" <rpenny at samba.org> > A: samba at lists.samba.org > Cc: "Daniele Bernazzi" <daniele at ao-siena.toscana.it> > Inviato: Mercoledì, 24 maggio 2017 13:12:07 > Oggetto: Re: [Samba] samba 4 in AD 2008R2 without winbind > > On Wed, 24 May 2017 12:54:48 +0200 (CEST) > Daniele Bernazzi

> Da: "Rowland Penny" <rpenny at samba.org> > A: samba at lists.samba.org > Cc: "Daniele Bernazzi" <daniele at ao-siena.toscana.it> > Inviato: Mercoledì, 24 maggio 2017 12:45:56 > Oggetto: Re: [Samba] samba 4 in AD 2008R2 without winbind > > On Wed, 24 May 2017 12:08:09 +0200 (CEST) > Daniele Bernazzi <daniele at ao-siena.toscana.it>

Hello, I write to know how can I modify the x axis : when I plot a survival object, R plots a graph with 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

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

Dear list, i'm a new R user, so I apologize if the topic is already being addressed by some other user. I'm trying to determine if the reproductive success of a species of bird is related to a list of covariates. These are the covariates: ? elev: elevation of nest (meters) ? seadist: distance from the sea (meters) ? meanterranova: records of temperature ? minpengS1: records

Dear R-Help Members, I am trying to read and analyse a set of 100 csv files. I need work only on some columns in each of those, so I decided to use a for loop, isolate the column in each file to work on, but then an error mesage appear: "undefined columns selected" Here is my code: setwd("F:/Data/") a<-list.files() for (x in a) { ? u<-read.csv(x, header=T,

Dear R users, we are trying to do some parallel computing using library(snow). In particular we have a cluster with 3 nodes >cl <- makeCluster(3, type = "MPI") 3 slaves are spawned successfully. 0 failed. and we want to compute the function op_mat (see below) first with the master and then with the cluster using system.time for checking the computational performance.

Dear Simon, Thank you for your quick reply! I used to perform the GAMsetup in the following manner: GAMsetup sintax: x.summer: vector used for construct the spline knots<-14 N<-length(x.summer) x<-array(x.summer,dim=c(1,N)) G<-list(m=1,n=N,nsdf=0,df=knots+1,dim=1,s.type=0,by=0,by.exists=FALSE,p.order=0,x=x,n.knots=knots,fit.method="mgcv") H<-GAMsetup(G) with the

Dear list members, First of all thank you for your helpful advices. After your answeres to my firt mail I studied a lot (R-News n?5) and I tried to perform my analysis: First, to fit a GLM with a nested design I decided to use the function "lmer" in package "lme4" as suggested by Spencer Graves and Filippo Piro. I remember to you that my data were: land use classes, 3 levels

Dear R-users, we are using library(snow) for computation on a linux cluster with RMPI. We have a problem with clusterEvalQ: after launching clusterEvalQ it seems loading the required library on each node but if we type a function belonging to the loaded package R doesn't find it. > library(snow) # making cluster with 3 nodes > cl <- makeCluster(3, type = "MPI") Loading

Dear R-users, I have a little problem that I can't solve by myself. I have a data frame with 2 factors and 8 observations (see the following code): y <- c(1,1,1,2,2,3,3,3) y <- factor(y) levels(y) <- c("a","b","c") x <- c(1,2,3,1,2,1,2,3) x <- factor(x) levels(x) <- c("x","y","z") X <-

Hi everybody, I'm a beginner in R and I'm having a hard time grouping my data by day. The data are in this format: id; words 2005-07-07T09:59:56Z; 35 2005-07-07T10:01:39Z; 13 2005-07-08T10:02:22Z; 1 2005-07-09T10:03:16Z; 23 2005-07-10T10:04:23Z; 39 2005-07-10T10:04:39Z; 15 I've transformed the date strings in dates with the function: london$id <- transform(london$id,

Dear R user, I'm trying to find a solution for optimizing my code. I have to run a 50.000 iteration long simulation and it is absolutely necessary to have an optimized code. I have to do this operation *sum_t ( t(X_t) %*% A %*% X_t )* where X_t is a (d*k) matrix which changes in time and A is a constant in time (d*d) matrix. I have put all my X_t in a three dimensional array X of dimension

? 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

Hi all, after upgrading from 13.7 to 14.2, asterisk cli (asterisk -r) don't show what's happens. I've trying setting debug and verbose to 100 but nothing, no show. All commands works as expected but i can't what's happens on my asterisk server. asterisk*CLI> core show settings PBX Core settings ----------------- Version: 14.2.0 Build Options:

$ sudo yum install devtoolset-7 Loaded plugins: fastestmirror, langpacks Loading mirror speeds from cached hostfile * base: mirror.vcu.edu * epel: mirror.siena.edu * extras: mirror.vcu.edu * updates: mirror.vcu.edu No package devtoolset-7 available. Error: Nothing to do On Wed, Sep 11, 2019 at 12:17 PM Gordon Messmer <gordon.messmer at gmail.com> wrote: > > On 9/11/19 9:05 AM,

Dear R-Users, I'm using SSPIR package for a spatio-temporal application. Is it possible to modify the structure of the involved matrixes (Fmat, Gmat, Vmat,Wmat)? I want to create a model like this #y(t)=k*theta(t)+epsilon(t) #theta(t)=h*theta(t-1)+eta(t) #epsilon(t) N(0,V) V=sigma2*I #eta(t) N(0,W) W=sigma2_eta where the state variable theta has dimension 1(p=1) and at

Hi all, I'm using RScalapack library for parallelizing some heavy matrix operations required by MCMC methods for spatio-temporal models. The package reference manuals (dated 2005) states that the library needs LamMPI to work but we have a Linux Cluster with OpenMPI. We have found (http://cvs.fedoraproject.org/viewvc/devel/R-RScaLAPACK/) a patch for OpenMPI but we are wondering if in the

Dear R users, we are trying to build a R package that includes a precompiled shared library, let's say mylib.so. We created the skeleton of the package and we moved the mylib.so file into the libs folder that we created at the same level of the folders man and R. Moreover we created the file NAMESPACE and we added the line useDynLib(mylib, .registration=TRUE). The building step seems to work