similar to: AR(1) with NLME

Hello! I'm having an echo problem with a TDM card. The TDM card is being fed by a channel bank just 12 or so feet away. When you put an analog handset on the line, both the RX and TX volume seem to be just fine. However, when I use the TDM card, I have to have an rxgain of 13.5, and even then, the audio is relatively quiet. I'm also getting echo on these lines, so I have turned

Hi there, I am using gls from the nlme library to fit an AR(1) regression model. I am wondering if (and how) I can separate the auto-correlated and random components of the residuals? Id like to be able to plot the fitted values + the autocorrelated error (i.e. phi * resid(t-1)), to compare with the observed values. I am also wondering how I might go about calculating confidence (or

Hi, I have a binary variable and corresponding predicted probability (using logistic regression on some explanatoey variables); I want to check that the model is well-calibrated using a calibration plot. how can I have the calibration plot for my data? thanks. [[alternative HTML version deleted]]

Hi R-sters, Just wondering what I might be doing wrong. I'm trying to fit a multiple linear regression model, and being ever mindful about the possibilities of autocorrelation in the errors (it's a time series), the errors appear to follow an AR1 process (ar(ts(glsfit$residuals)) selected order 1). So, when I go back and try to do the simultaneous regression and error fit with gls,

Hello, I am attending a course in Computational Statistics at ETH and in one of the assignments I am asked to prove that a time series is not autocorrelated using the R function "acf". I tried out the acf function with the given data, according to what I found here: http://landshape.org/enm/options-for-acf-in-r/ this test data does not look IID but rather shows some trends so how can I

Dear R users: In Prof. Harrell's library rms, calibrate.rms plot the Bias-corrected Probability and Apparent Probability. The latter one can be retrieved from class calibrate.default. But how to retrieve the former one. BW *Yao Zhu* *Department of Urology Fudan University Shanghai Cancer Center Shanghai, China* [[alternative HTML version deleted]]

Dear R users, someone knows how to remove auto-correlation from a frequencies time series? I've tried by differencing (lag 1) the cumulative series (in order to have only positive numbers) , but I can't remove all auto-correlation. If it's useful I can send my db. x <- # autocorrelated series new1<-cumsum(x) new2<-diff(new1,lag=1,differences = 1) acf(new2) #

Hello, currently I'm working on a model based on Monte-Carlo-Simulations. I observed that a generated normal distributed times series using rnorm(100,mean=0,sd=1) is far away from being not autocorrelated. Is there any other gerenator implemented in R, which might solve my problem? -- View this message in context:

Hola a todos, Estoy iniciandome en R y la verdad es que aun estoy muy muy verde.... Estoy intentando clasificar unos datos con support vector machines, pero me da fallo al usar la funcion predict. El código que uso es el siguiente: calibrate<-read.table("calibration.txt", header=TRUE) calibrate$calibration<-as.factor(calibrate$calibration)

Hi all, I have a question on the package « survey” I have some difficulties to use the function ‘calibrate’. Although it works well with one single factor variable, I cannot use it for 2 and get the message “Erreur dans regcalibrate.survey.design2(design, formula, population, aggregate.stage = aggregate.stage, : Population and sample totals are not the same length.” Here is the format

Hi, I am trying to learn time series, and I am attending a colleague's course on Econometrics. However, he uses e-views, and I use R. I am trying to reproduce his examples in R, but I am having problems specifying a AR(1) model. Would anyone help me with my code? Thanks in advance! Reproducible code follows: download.file("https://sites.google.com/a/proxima.adm.br/main/ex_32.csv

Hola de nuevo, Esta es la información de mi sesion: R version 2.10.1 (2009-12-14) i386-pc-mingw32 locale: [1] LC_COLLATE=Spanish_Spain.1252 LC_CTYPE=Spanish_Spain.1252 [3] LC_MONETARY=Spanish_Spain.1252 LC_NUMERIC=C [5] LC_TIME=Spanish_Spain.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base Lo que yo prentendo es

Hello, I'm a newby in R. I have created a data.frame holding panel data, with the following columns: "id","time","y", say: periods = 100 numcases = 100 df = data.frame( id = rep(1:numcases,periods), time = rep(1:periods, each = numcases) ) df = transform(df,y=c(rnorm(numcases*periods)+id) I want to check whether "y" is autocorrelated. I came across

Hi, I would like to specify a spherical correlation structure for spatially autocorrelated residuals in a model based upon the logistic function of a response that is a proportion (0 to 1) (so usual binary logistic regression is not an option). There is no need for a g-side random effect with grouping in this model. Am I correct that nlme requires this (meaning a correlated error structure only

Good afternoon! As mentioned in the subject, my question regards more the methodological part that accompanies survey design and the statistical part that is involved. So, I have the following data: a<-data.frame (id_hh=c(1:5), strata=c(1,1,2,2,1), Nhstrata=c(100,100,200,200,100), Nrmemb=c(2,4,2,5,4)) a$ocmemb1<-c("wk","jl","st","jl","st")

Hi R users, I'm here trying to understand correlated residuals in nonlinear estimation. I'm reading/studying the book Bates, D. M. and D. G. Watts, (1988), /Nonlinear regression analysis and its applications/, Wiley, NY. pages 92-94, trying to reproduce the figures and to find out the code in R to perform the necessary calculations. I also consulted Pinheiro and Bates, but without

Hello, I am using the val.prob function in the Design package. I understand how the Brier quadratic error score is calculated, but I do not know how the Brier score computed on the calibrated rather than raw predicted probabilities (B cal) is calculated. My question is: how are the calibrated probabilities calculated? Any explanation of this, or references to explanations of this, would be

I'm trying to do a simple (?) analysis of a 1D spatial data set, allowing for spatial autocorrelation. (Actually, I'm comparing expected vs. observed for a spatial model of a 1D spatial data set.) I'm using models like gls(obs~exp,correlation=corExp(form=~pos),data=data) or gls(obs~exp,correlation=corLin(form=~pos),data=data) This form is supposed to fit a linear model of

Hi, I have spatially autocorrelated data (with a binary response variable and continuous predictor variables). I believe I need to do an autologistic model, does anyone know a method for doing this in R? Many thanks C Bell

Hi! How can I exclude the legends from calibration plots  generated by calibrate.cph regards, Salvo [[alternative HTML version deleted]]