similar to: anova p value extraction

Hello everyone. I am trying to build up understanding in R by trying to develop just some simple scenarios. I would like to explain you what I am trying to do and what I did so far. I would like to put inside a RasterLayer (raster package) a Gaussian field (for given covariance) using grf function (geoR package) 1. First I created a Raster Layer object r <- raster() # Default values are

Hello, I am a graduate student of University of Florida. I am trying to run a process variation model and I am facing the following error. ------------------------------------------------------------- Analysis of geostatistical data For an Introduction to geoR go to http://www.leg.ufpr.br/geoR geoR version 1.6-27 (built on 2009-10-15) is now loaded

Dear friends. I have trobpe with path. Below I have tried to paste a path and file name but it will not work > path <- "D:/rw0641/own/own/procedures/GFR_TIME" > path [1] "D:/rw0641/own/own/procedures/GFR_TIME" > path1 <- paste(path,"gfr.txt",sep="/") > path1 [1] "D:/rw0641/own/own/procedures/GFR_TIME/gfr.txt" -- this is good

Hello all, I am trying to create a grid of large number of points 4096*4096, on processing the data I am writing it into a file. phi <- 0.5 N <- 4096 mu <- 90 sim<-grf(N*N,grid="reg",cov.model="spherical",cov.pars=c(1,phi),method="RF") sim$data <- (sim$data - mean(sim$data)) / sd(sim$data) local <- array(sim$data, dim=c(N,N)) sysvar <- local *

I have asked the same question on stackoverflow but did not get a satisfying answer. I am trying to simulate a lognormal spatial random field but I need the simulated value in a certain range. So I need some easy to use functions to generate a truncated Gaussian field to start with. To be specific, I need a function like GaussRF from the RandomFields package or grf from the geoR package to

Hi, I want to automatically fit variograms to a large number of different sample data sets, and call the funtion "likfit" (in package geoR) from within a for-loop. "likfit" does again call "optim". After ssuccessfully fitting variograms to some of the data sets, the procedure crashes and I get the error message: Error in optim(par = ini, fn = negloglik.GRF,

Been on linux (ubuntu 9.10) for about a month now. Starting to get the hang of it. The only thing from windows I am missing is being able to play this game. I've read a lot of stuff on the net and tried different things. Some say i need to trick the ip but I am unable to figure out the server IP address. A lot of private servers have certain things in place to prevent hacking, so I think

Hi all, I would like to test the relevance of a vector field (i.e. if the vectors are organized or not). To do so, I would like to use an autocorrelation test, so that I have two questions: - is the Watson test applicable to that perpuse ? - is the kuiper test applicable to that purpuse ? Regards ------------------------------------------------ Emmanuel Poizot Cnam/Intechmer B.P. 324 50103

Dear R People: I know that this is COMPLETELY off topic. Does anyone know how to put colors into a LaTex document, please? This is LaTex via WinEdt. I can get colors via PcTex but not in WinEdt. Thanks in advance! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: hodgess at gator.uhd.edu

dear list I am student M.S. statistics in department statistics . I am working in the function "nls" in the [R 2.3.1] with 246 data and want to fit the "exp" model to vectors( v and u ) but I have a problem to use it u 5.000000e-13 2.179057e+03 6.537171e+03 1.089529e+04 1.525340e+04 1.961151e+04 2.396963e+04 2.832774e+04 3.268586e+04 3.704397e+04 4.140209e+04

Hiya, I've posted this 'enigma' on the newsgroup : http://groups-beta.google.com/group/comp.emulators.ms-windows.wine/browse_thread/thread/57dcb7531ebfb529/65018c77a5e58651#65018c77a5e58651 The answers were pretty rare so I post the problem here, too. " I'm getting connection problems with an application run under Wine which isn't acting like that under a native Windows

Dear R developers A question about the class 'listof', defined in package 'stats'. Other than its definition and use in the code for 'anova', we can't see that the class 'listof' is used for anything else (in recommended packages, or elsewhere). In the spatstat package we have been using a 'listof' to represent a list of spatial objects of the same

I'm trying to write a function to automate doing a variance analysis, part of which involves doing some further calculations. The method I've been using isn't very robust, if variable names change then it stops working. For this dummy data > dput(assayvar,"") structure(list(Run = structure(c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L, 4L), .rk.invalid.fields =

Hi All , Can anybody tell me if there's any way to get the summarized anova values.Now i will explain what i mean , when i say "*summarized*". Below you can see the anova table of recmeanC1 with rest* all* i.e from recmeanC2 to i15(predictors),as shown in table. Df Sum Squares Mean Square F value Significance [Pr(>F)] recmeanC2 1 89.272 89.272

Hi, all I've written some R script to calculate the linear regression of a matrix. Here below is my script: >x<-matrix(scan("h:/data/xxx.dat",0),nrow=46,ncol=561,byrow=TRUE) >year <- NULL >year <- cbind(year,as.matrix(x[,1])) >lm.sol<-lm(x~year) >xtrend<-coef(lm.sol)[2,] # get the matrix of regression coefficient >t.test<- ?

Hello, I have a list of data with multiple elements, and each element in the list has multiple variables in it. Here's an example: ### Make the fake data dv <- c(1,3,4,2,2,3,2,5,6,3,4,4,3,5,6) subject <- factor(c("s1","s1","s1","s2","s2","s2","s3","s3","s3",

Dear r-users, Suppose I have three datasets: Dataset-1: Date x y Jan-1,2005 120 230 Jan-2,2005 123 -125 Jan-3,2005 -110 300 Jan-4,2005 114 -21 Jan-7,2005 11 299 Mar-5,2005 200 311 Dataset-2: Date x y Jan-2,2005 123 -125 Jan-3,2005 -110 300 Jan-4,2005 114 -21

The PDF file on the web, which is an appendix on nonlinear regression associated with the CAR book, is very nice. When I ran through the code presented there, I found something odd. The code does a certain model in 3 ways: Vanilla NLS (using numerical differentation), Analytical derivatives (where the user supplies the derivatives) and analytical derivatives (using automatic differentiation). The

Dear all, I'm trying to access the Sums of Squares resulting from a summary(aov(....)) so I can use them in a function. Is there an easy way to access these? The structure of a summary.aov object is something like this: > str(summary(fhc.aov)) List of 5 $ Error: PPNR :List of 1 ..$ :Classes anova and `data.frame': 1 obs. of 5 variables: .. ..$ Df : num 70 .. ..$

----=_NextPart_ST_21_29_48_Monday_February_23_2004_11097 Content-Type: text/plain; Your message could not be sent. A transcript of the attempts to send the message follows. The number of attempts made: 1 Addressed To: john@spaar.com Mon, 23 Feb 2004 21:29:47 -0600 (2) Failed to send to identified host, [10.0.0.2], 550 5.1.1 User unknown [SMTP server error] --- Message non-deliverable.