similar to: modified R/S

Dear R-users, I have some problems with the geweke.diag-Function of the coda-package. I try to obtain the Geweke-diagnostic by using the following, simple code: library(coda) input.geweke = as.mcmc(input.matrix) #input.matrix is a 25.000 x 35 matrix with the 25.000 saved draws of the 35 parameters of interest output = geweke.diag(input.geweke) However, I get the following

Hi all: The CODA package provides Geweke Diagnostic method for convergence checking. The geweke.diag in CODA returns Z-score value but not give a conclustion that it is convergence or not. So I'd like to know how small/big the magnitude of Z-score is corresponding to the convergence of a chain. That is, Doese Z-score smaller or more than *threshold *determine the convergence? If so, how big

Intermittently I'm getting this error from the geweke.diag function of the coda package. Would anyone be kind enough to enlighten me as to the possible source of such an error, or how to debug/locate it? Error in { : task 22 failed - "inner loop 1; cannot correct step size" Thanks a bunch. J -- Dr. Jim Maas University of East Anglia

Hi R users, Does anybody knows for the following erro after running geweke.diag(MCMC.sampled, frac1=0.1, frac2=0.5) Erro em glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : la??o interno 1; n??o ?? poss??vel corrigir o tamanho do passo Al??m disso: Warning messages: 1: algoritmo n??o convergiu in: glm.fit(x = X, y = Y, weights = weights, start = start,

Hi samba lists, we have a samba-ldap domain running on a debian squeeze (samba 3.5.6)server (pdc and bdc). I try to configure a domain member server on an other debian squeeze that will serve as cifs and nfs server. My Debian server member use winbind (on ldap) for mapping the users windows sid to the unix uid. The users mapping is write in the the ldap directory : ou=idmap,dc=exemple,dc=com

Dear R users: I am using hierchical modeling for a response varible (ordinal) given by patients for each doctor. My MCMC run with MCMMCglmm ( thin=20,nitt=208000,burnin=24000, family='ordinal') still can't have all the parameters pass the GEWEKE diagnosis test. I am trying to find the reason for the slow convergence. Can the reason be 20 out of 50 doctors having 5 or less response

Si vous ne visualisez pas correctement l???e-mail, cliquez ici <http://url.mp32.ch/visu-90BD7B44-1B18-42B2-BB27-BF4B02C4DE10-203221958-356307-08042014.html> D??couvrez les nouveaux programmes immobiliers du Grand Ouest <http://url.mp32.ch/url-203221958-1612521-08042014.html> <http://url.mp32.ch/url-203221958-1612521-08042014.html> a s??lectionn?? pour vous les meilleurs

Si vous ne visualisez pas correctement l???e-mail, cliquez ici <http://url.mp32.ch/visu-C2035EE0-46EE-4FCF-9C8F-E48C8E5E008C-203221958-355072-03042014.html> RETROUVEZ-NOUS AU SALON NATIONAL DE L'IMMOBILIER ?? PARIS <http://url.mp32.ch/url-203221958-1607886-03042014.html> <http://url.mp32.ch/url-203221958-1607887-03042014.html> a s??lectionn?? pour vous les meilleurs

Si vous ne visualisez pas correctement l???e-mail, cliquez ici <http://url.snd10.ch/visu-E9D92C62-F56A-4768-8594-94E7BC64B1C3-203221958-423876-31072014.html> HABITER OU INVESTIR SUR LE LITTORAL ATLANTIQUE <http://url.snd10.ch/url-203221958-1804885-31072014.html> <http://url.snd10.ch/url-203221958-1804885-31072014.html> a s??lectionn?? pour vous les meilleurs programmes

Hi, I try to configure a samba 4.1.17 server member of a samba3.5.6 (ldap backend ) domain. i can mount the smb share on a windows client but i cant modify the acl. this is the samba server logs: create_canon_ace_lists: unable to map SID S-1-5-21-856890099-1868262392-538272213-2012 to uid or gid when i try to manually map SID to UID with wbinfo: wbinfo -S

Hi: I have a question regarding printer support/services. Does Samba support shared printers directly attached to a Windoz or Linux machine? (i.e. attached via a parallel port or USB port to either the Samba machine or one of the Linux or Windoz boxes). If so, how does one configure Samba and/or the other network machines to configure the printers to make them reachable to the rest of the

Dear List: the CODA and BOA packages for the analysis of MCMC output yield different results on two dignostic test of convergence: 1) Geweke's convergence diagnostic; 2) Heidelberger and Welch's convergence diagnostic. Does that imply that the CODA and BOA packages implement different ``flavors'' of the same test? I paste below an example. Geweke's test

Hello, I'm using the fArma package to estimate Hurst exponent by R/S method. I've some measurements in the following format: Call_number Call_duration I'm using the following steps. Am new to R, so it would help if someone could please confirm if my steps are correct. Further, this method seems to give a value H ~= 0.8 (> 0.5). Is the use of R/S method to estimate H OK? >

Hello, I'm using fArma package to estimate the value of Hurst exponent using R/S method. However, for a certain set of data I get H ~ 1.8. How do I interpret this? Following are the output that I get for this set: > mean(data[,2]) [1] 400.5433 > sd(data[,2]) [1] 1139.786 > > rsFit(data[,2], levels = 64) Title: Hurst Exponent from R/S Method Call: rsFit(x = data[, 2], levels

Dear R users, The function lmacfPlot( ) in fBasics, returns a list including the Hurst Exponent. I get sometimes a value for Hurst Exponent using lmacfPlot bigger than 1 which is theoretically incorrect. What could be reason ? From the related documents, I could not find how the Hurst Exponent in lmacfPlot is estimated. Could you please say how it is estimated here? Thank you so

Hi, i agree with Stefan. I have the same problem with Win2K Client and SAMBA Server. But on my situation i have one PC that i can make a connect and all the other not. Alle workstation are on the same group, are all in the same SubNet; but i can?t see my SAMBA _ Server so well like on this machine that is working. So is there any other problems on SAMBA ? Kind regards Matthias

Dear R Users, I am using code from the following links to do Hurst exponent estimation. link: http://r-forge.r-project.org/plugins/scmsvn/viewcvs.php?rev=1&root=rmetrics&view=rev file: LongRangeDependence.R Take a look at the following run: > x<-(cos((1:200)/10)) > rsFit(diff(x,15))@hurst$H [1] 1.027420 > aggvarFit(diff(x,15))@hurst$H [1] 0.02301331 First of all, can anyone

Dear All, I am a novice when it comes to time-series analysis and at the moment I am actually interested in calculating the Hurst exponent of a time series. This question has already been asked quite some time ago http://bit.ly/98dZsi and I trust some progress has been made ever since. I was able to find some functions in the packages http://cran.r-project.org/web/packages/Rwave/index.html

Hi, Can u please tell me which all packages do i need to install to estimate the hurst parameter in R. I have tried installing all the possible options but still it doesnt work. basically i want to use 9 functions to estimate hurst parameter like aggvarfit, rsfit, etc. i will be very thankful if u could be of some help. -- Regards, Deepak Jadhav. [[alternative HTML version

Dear All; I have used fbmSim to simulate a fbm sequence, however, when I tried to estimate the Hurst effect, none of the nine procedures gave me an answer close enough to the real value, which is 0.5 (n=1000). So, would you please advice, 1. which is the best method to estimate the H among the 9 mehods, R/S, higuchi or Whittle? 2. how to choose the levels (default=50), minnpts, cutoff values or