similar to: Residuals in gls

Dear all, I am trying to create a vector name, for example tmax.195012 from tmax., 1950 and 12. Obviously I don't wish to simply type it because the 3 name components are changing in each iteration within a loop. Is there any way of concatenating those 3 components (which are a mixture of numbers and letters)? Thanks for reading, Panos

Dear R colleagues, I am trying to plot some geophysical data as a filled contour on a continent map and so far the guidance from the R-help archives has been invaluable. The only bit that still eludes me is the colour key (legend) coming with filled.contour: I prefer to generate my own colour palette, mainly based on the quantiles of tenths of the data in order to capture the whole range (of

Dear R people, I notice that the confidence intervals of a very small sample (e.g. n=6) derived from the one-sample wilcox.test are just the maximum and minimum values of the sample. This only occurs when the required confidence level is higher than 0.93. Example: > sample <- c(1.22, 0.89, 1.14, 0.98, 1.37, 1.06) > summary(sample) Min. 1st Qu. Median Mean 3rd Qu. Max.

Hi to everybody, i tried to find an asnwer before posting this, but most astcc billing issues i searched refer to the case when no billing occurs at all. In my case i get only initial charges and any subsequent minute does not count for billing. In my iax.conf i entered the "notransfer = yes" but nothing changed. My test card and test calls are the following TEST-CARD en N/A

Dear Rers, I have threes series, x, y, z and I want to fit a model z ~ x + y. First of all, I fit a lm. I found the residuals are correlated, by looking at the acf() and pacf(). Then I tried to fit a gls model allowing residuals to be correlated (correlation = corARMA(p=5, q=1)): y.na <- as.data.frame(y[complete.cases(y),]) y.gls <- gls(z ~ x + y, data = y.na, correlation=corARMA(p=5,

Under the description of the 'type' argument, ?residuals.gls says 'Defaults to ?"pearson"?.' But residuals.gls starts residuals.gls <- function(object, type = c("response", "pearson", "normalized"), ...) { type <- match.arg(type) ... which sure looks to me like it defaults to "response", not "pearson" (and it

Hello all I wanted to forward all incoming requests in port 80 to a server in my LAN, and by using DNAT lines, it actually works. However, it is unstable, in the sense that in the beggining of each connection (one or two seconds) it is extremely fast, then it sometimes pauses and waits 30 seconds or so, then it starts again and so on. The line i used is : DNAT net loc:192.168.0.210 tcp 80 DNAT

Dear all, I am a beginner user of R and I tried to fit a gls model with explanatory variables and an AR(1) correlation component using the function "gls" with: correlation = corAR1 (form = ~ 1) It should mean that the residual follows an AR(1) process, isn't it? The problem is that, if I use the funcion "predict" I noticed that the predicted values are the same as if I

Dear colleagues, hi. I am using the function gls in R-package nlme and I want to compute either the leverage values or the studentized residuals. Can any body tell me how can do that? With regards ___________________________________________________________ 围蟻畏蟽喂渭慰蟺慰喂蔚委蟍[elided Yahoo spam]] 螔伪蟻蔚胃萎魏伪蟿蔚 蟿伪 蔚谓慰蠂位畏蟿喂魏维 渭畏 未喂伪胃苇蟿蔚喂 蟿畏谓 魏伪位蠉蟿蔚蟻畏 未蠀谓伪蟿萎 蟺蟻慰蟽蟿伪蟽委伪 魏伪蟿维 蟿蠅谓 蔚谓慰蠂位畏蟿喂魏蠋谓 渭畏谓蠀渭维蟿蠅谓

I am not sure if anyone has the answer. Here is the setup: I am using redhat linux 9 with the latest rpm from www.samba .org (3.0.6). I am connecting to the linux server with a G4 running 10.3.5. Once the share is mounted on the G4 using the appropriate smb://.etc command. You get the little global and everything is good. Opening the global will give you 2 folders. Disk1 and Disk2. Both these

I can't seem to get rsync to restart where it left off when I am syncing a large file (> 5GB). Below is some info on what I have been doing. If someone has the energy to barrel through my comments that would be great. Out of curiosity is there an alternative to rsync for large files? I believe I looked all over the web for answers and I found "partial" answers and just

Hello I got about 20 different people using 20 different PC''s in my Lan. Five of them should have unrestricted access to the internet any time of the day. Yet the rest of them should only have access at certain timeframes during the day. Two of them should have no access at all. I am well experienced in writing rules (for the rules file) that make all this possible, yet i was

I am using the gls function to fit a two-stage least squares model with first order autoregressive error terms. Since there is no automated adjustment for the use of two-stage least squares in this package, I am trying to manually replicate standard errors of the coefficient estimates in order to adjust for a first stage OLS estimate of endogenous variables. However, thus far I have been unable to

Dear all, I am facing a problem with running a script that i have written in R and i hope you can help me spot which exactly is the problem and how i can solve it. The error I get is the following: ---- Error in .jcheck(silent = FALSE) : Java Exception <no description because toString() failed>.jcall(row[[ir]], "Lorg/apache/poi/ss/usermodel/Cell;", "createCell",

Hello to all, I am trying to build speex for the Windows Mobile platform. However I am not able to build it even on Windows x86 because a source file (medfilter.c) while declared in the project settings is missing. So linker cannot find medfilter.obj and I am not able to build it. Please could anyone tell me if this is an error that I should fix manually (for example ignore the file and fix

I am trying to analyse a data with gls/lm using the following set of models prcn.0.lm <- lm( log10(Y)~(cond-1)+(cond-1):t ,prcn) prcn.1.gls <- gls( log10(Y)~(cond-1)+(cond-1):t ,prcn,cor=corAR1()) prcn.0.gls <- gls( log10(Y)~(cond-1)+(cond-1):t ,prcn) prcn.1m.gls <- gls( log10(Y)~(cond-1)+(cond-1):t ,prcn,cor=corAR1(),method="ML") I get the following AICs for these models:

Hello listmembers, I've been thinking of using gls in the nlme package to test for serial correlation in my data set. I've simulated a sample data set and have found a large discrepancy in the results I get when using the default method REML vs. ML. The data set involves a response that is measured twice a day (once for each level of a treatment factor). In my simulated data set, I

Dear all, I have been using the function gls in the package nlme in R to fit some spatial regressions (as described in Dormann et al.). However, I have been struggling trying to find a way to calculate a measure of effect size from these models, so I wanted to know if any of you had an idea on how to do this. More precisely, I am producing a multiple model with an exponential correlation

Hi, I have fit a model using the rms package with the Gls() function. Is there a way to get the model estimates, std errors, and p-values (i.e. what you get with print(fit)) into latex format? I have tried: f <- Gls(...) latex(f, file='') ... but I get the following error Error in replace.substring.wild(s, old, new, test = test, front = front, : does not handle > 1 * in

Hi I have a huge array with series of data. For each cell in the array I fit a linear model, either using lm() or gls() with lm() there is no problem, but with gls() I get an error: Error in glsEstimate(glsSt, control = glsEstControl) : computed gls fit is singular, rank 2 as soon as there are data like this: > y1 <- c(0,0,0,0) > x1 <- c(0,1,1.3,0) > gls(y1~x1)