similar to: Error Message from Variogram.lme Example

HI Dave, My comment was based? on: " >The main question with this test was if the interaction term is significant (i.e. growth rate). However, my question is could I also look at the p-values of the main effects to ?>say if body mass increase significant with body mass?" Here, the result shown were from the summary of the linear model.?? We report the p-values of the main

Dear R-help, The following nlme example, from help(Variogram.gls), does not work on my system. > data(BodyWeight) > fm1 <- gls(weight ~ Time * Diet, BodyWeight) > Variogram(fm1, form = ~ Time | Rat)[1:10,] Error in FUN(X[[1]], ...) : unused argument(s) (method ...) It looks like I might not have loaded a necessary package. I did load nls and lattice (the latter not actually being

Dear R-help, In nlme package, anova () can be used to test the difference between two coefficients as shown on page 225 of "Mixed Effects Models in S and S-Plus": anova(fm2BW.lme, L = c("Time:Diet2" = 1, "Time:Diet3" = -1)) Now my question is instead of test the difference between two coefficients, can I use anova to test the significance of the sum of two

Les consulto por lo siguiente: Se me ocurrió usar datos de MCMCglmm en otros procedimientos de R, comparar los resultados y el tiempo como facilidad de uso, para lo cual use LME4, los resultados fueron parecidos, entendiendo que internamente los algoritmos tienen diferencias, es lógico que los resultados no coincidan exactamente, pero aproximadamente tienen que ser compatibles, en este caso

  Hola,   Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?   http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5   Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.   Muchas gracias.   Eva [[alternative HTML

Dear all, This question may be too basic quesition for this list, but if someone has time to answer I will be happy. I have tried to find out, but haven't found a consice answer. As an example I use "Pinheiro, J. C. & Bates, D. M. 2000. Mixed-effects models in S and S-PLUS. Springer, New York." page 225, where rats are fed by 3 different diets over time, which body mass has

Dear Dr. Bates, and R-help, I've tried the anova function to test the difference between two coefficients, as shown on page 225 of your book "Mixed Effects Models in S and S-Plus (Statistics and Computing)". When I type: anova( fm2BW.lme, L = c(TimeDiet2 = 1, TimeDiet3 = -1) ) I got the following error message: Error: unexpected '=' in "anova( fm2BW.lme, L =

Hello, I am new to R and have a couple of questions. My data set contains the variables "Bwt" and "Hwt", which are bodyweight and heartweight, respectively, of a group of cats. With the following code, I am making two plots, both to be viewed in the same plot window in R: library(MASS) maleData <- subset(cats, Sex == "M") linreg0 <- lm(maleData$Hwt ~

Dear Members, I am trying to fit a variogram model using fit.variogram function from the gstat package. The figure showing my experimental variogram can be seen here: https://i.stack.imgur.com/UZXw4.png My code line for this operation is: > c2.vgm.fit<-fit.variogram(c2.vgm.exp,vgm(nugget=0, psill=400,model="Exp",range =40000),fit.method = 7) The system throws following error

Dear All, Is there anybody has the experience in using variogram(gstat) ? Please kindly give me some hints about the results. I used variogram() to build a semivariogram plot as: tr.var=variogram(Incr~1,loc=~X+Y,data=TRI2TU,width=5) then fir the variogram to get the parameters as: v.fit = fit.variogram(tr.var,vgm(0.5,"Exp",300,1)) v.fit model psill range 1 Nug 1.484879

Dear all can you help me? i have this problem: i have a dataset in a text file in a matrix of 3 columns: x, y, z where x and y are the coordinates and Z are the mesurements. How can i do a variogram with R?

Hello, I have data at 10 locations, in each location there are time series (T=56). Question is: when I?m fitting variogram what happens with those measures in each location? Are they taken as repeated measures? It's very important for my to know this Thanks a lot

Hi: I need some help with the legend. I got 14 samples(Muestreo) and I am trying to plot a smooth line for each sample. I am able to accomplish that but the problem is that the legend only displays every other sample. How can I force the legend to show all of my Muestreos? Thanks in advance. fish_ByMuestreo <- structure(list(data = structure(list(SampleDate = structure(c(3L, 3L, 3L, 3L,

I wonder what is the default value for the argument 'cutoff' when not specified in the variogram.formula function of gstat. Computing variogram envelops within gstat, I am comparing the results obtained with variog in geoR and variogram in gstat, and it took me a while before understanding that the cutoff default value is not the maximum distance. Can Edzer tell us about it? All the

Hello, I'm french and I have some difficulties in carry out the semiones under R with an aim to carry out an interpolation by krigeage. My goal is to obtain a chart of the distribution of precipitations/temperatures in Europe starting from 73 different stations (and, of course, distributed irregularly on the chart, where use of the krigeage). Here, I carried out this to test to obtain

Hi Im attempting to build a sample variogram for 300 obersvations of longitudinal data. So what I need to do is compute the half squared differences between pairs of residuals (for instance if a subject has 4 obersvations, this is 4 choose 2 paird differences) for each subject. Also, then I need the corresponding time differences within each individual. So the end result will be a 300 by 2

hello I need basic explanation of variogram and how to use in R. I will appreciate if you supply them. regards Ahmet Temiz ______________________________________ XamimeLT - installed on mailserver for domain @deprem.gov.tr Queries to: postmaster at deprem.gov.tr ______________________________________ The views and opinions expressed in this e-mail message are ...{{dropped}}

Hi I'm interested in doing a anisotropy analysis of some spatial data and I want to do a contour plot of a covariance or variogram surface to visualize the directions of maximal and minimal anisotropy. I know how to do a contour plot :-) but how can I obtain the variogram surface having a geoeas format dataset (X, Y, Z) ? Thanks EJ

Hello everybody, I want to run the following R-code, but it does'nt work. > res.vgm = variogram(ertrag ~ rep, # !!! + loc = ~ xpos + ypos, # !!! + width = start.range / 10, + data = d) Error in vector("double", length) : element 2 is empty; the part of the args list of '.Internal' being evaluated was:

Hello, I am attempting to get the standard deviation in multiple distance bins in my spatial data. It appears as though the 'variog' command in the geoR package will do the trick, as one of the outputs from 'variog' is 'variog$sd', which, according to the manual, is the "standard deviation of the values in each bin". However, when I run this command, the