search for: maternity

Hi all, I tried to produce some kriged surfaces with geoR (latest version). The size of the grid should be around 900 x 650 cells (what I find is not a very big grid), and the number of points is around 2500. The command krige.conv stopped after arround 5 min saying it can not allocate a vector with around 1.5 billion units. Sounds reasonable. Is there a workaround? How would I partition the

Hi all, >str(df) `data.frame': 31837 obs. of 3 variables: $ x : num 410683 410700 410720 410740 410324 ... $ y : num 43136 43126 43123 43125 42709 ... $ wz: num -101.1 -94.9 -93.3 -94.5 30.8 ... >library(gstat) >g<-gstat(id="rv",form=wz~1,loc=~x+y,data=df,model=mat,nmax=500,set=list(average=1)) >str(g) List of 3 $ data :List of 1 ..$ rv:List of 10 ..

Hi, As part of the analysis of a longitudinal study, I have a list, called wg, on R 1.3.1 for windows created as follows - wg <- by(xx, xx$mrn, function(xx) lm(xx$maternal.wt ~ xx$gestation, data=xx)) xx is a data.frame, and so wg contains the results of 293 linear regressions, as lm objects. Now what I want is to extract the 293 intercepts and regression coefficients as a vector, and

Hi, R-folks: I have been tryin many combination of parameter to make Matern variogram to work, but I can't find the available one. I'm near to be crazy. I tiped: A?o2003Selg.lf<-likfit(A?o2003Selg,cov.model="matern",ini.cov.pars=c(1.5,14),kappa=2.5,fix.kappa=FALSE,nugget=0.08,lambda=0.008,fix.lambda=FALSE,hessian=TRUE) the hessian shows: $hessian [,1]

Dear All, I am analysing a dataset on levels of herbivory in seedlings in an experimental setup in a rainforest. I have seven classes/categories of seedling damage/herbivory that I want to analyse, modelling each separately. There are twenty maternal trees, with eight groups of seedlings around each. Each tree has a TreeID, which I use as the random effect (blocking factor). There are two

i'm trying to build a little summary table for the contents of a data frame. t<-sapply(macro, data.class) c<-sapply(macro, length) m<-sapply(macro, mean, na.rm=T, digits=2) cbind(type=t, n=c , mean=m) I want to drop the variables that are factors so I can include -max- and -min- in my table. -macro- contacts the data--how do I drop the variables according to their data.class

Hi, I would be grateful for help with a problem which is irritating me. I am quite sure that I am doing something stupid, but I can't see what it is. I am running R 1.7 on Windows 2000. The graphics device is the PC screen. The graphics from the nlme demonstration in Bates an Pinheiro's manual work just as advertised. The CO2 data and the Orthodont data dsiplay

Hi, I'm a long-time STATA user and a R newbie. I'm doing ok, but I'm addicted to STATA macro variables. Is there something like a macro variable in R? Specifically, I'd like to be able to do something like for (i in 1:3) { ..... x`i' <- ... } where R would resolve x`i' to the objects named x1, x2 and x3 as I move through the loop. I guess I could create these in

ComRades: I am geeting the error message Error in ldots[[which(MET)]] : attempt to select less than one element when I try to fit the geostatistical model with the likfit() function of geoR. I have tried with old data for which likfit() successfully maximised the likelihood in previous versions of geoR, and yet the current version fails. I have tried in Windows Vista and Windows XP (I haven't

I am trying to generate two dimensional random coordinates. For randomly distributed data I have simply used >xy<-cbind(runif(100),runif(100)) However I also want to generate coordinates that are more uniformly distributed, and coordinates that are more contagiously distributed than the above. Can anyone make any suggestions Thanks. Dr Terry Beutel Rangeland Scientist Animal

Hi, I'm moving (slowly) to R from STATA. I often have need to move through a set of tasks across a series of years. In this case, you can see that I'm mimicking -reshape- in STATA, but I'm less interested in the task than in programming R. library(foreign) mydata<-read.dta("z:\example.dta") for (y in 2000:2002) {

-- Dear R users, I have data on around 2000 birds from 3 generations for which I know an individual's pedigree (i.e. the relationship it shares with other individuals e.g brother, uncle, mother) and also a pedigree based on foster-families, because half broods were removed from their nest of origin and placed in a foster parent's nest. From this I want to model two types of random

-------- Original Message -------- Subject: [AGDG-LIST:428] Summer Course in Guelph Date: Fri, 22 Dec 2006 09:12:24 -0500 From: Larry Schaeffer <lrs at uoguelph.ca> Reply-To: lrs at uoguelph.ca To: Animal Geneticist's Discussion <agdg-list at colostate.edu> The Centre for Genetic Improvement of Livestock at the University of Guelph is pleased to announce a one week summer course

I'm trying to get my graphics so that I can use them in LaTeX to create (via ) a pdf presentation. I've tried controlling inner and outer margins and figure size using par(), to no avail. The ps output keeps appearing as a portrait page with a centered figure. Nothing I have been able to do so far has changed that. Below is the code: par( bg="yellow", lab=c(10,5,7)

Hi I'm using geoR package to perform linear spatial interpolation(OK). The function likfit() fails to compute REML. The error meassage is : Error in solve.default(v$varcov, xmat); How I can find out that likfit() is failed to process and retrieving the error message ? Thank you so much for your help.

Nir Simionovich is about to become a father. He will be joining our conference at 12 Noon EDT today from the Maternity Ward to talk about Amazon EC2 cloud computing with Asterisk. Nir gave a very good presentation on this at AMOOCON a few weeks ago (see http://www.amoocon.de for more on that). The advantage here though is that he'll be live with us for your questions. All the details on how to join us are here...

I have a problem with the textmatrix() function of the LSA package whenever I specify 'removeNumbers=TRUE'. The data for the function are stored in a directory LSAwork which consists of a series of files that houses the text in column form. As long as removeNumbers = FALSE or it is not present the textmatrix function works just fine. The error message I get seems to suggest it is

Hello everyone, at the moment I'm using the tgp package for modelling a nonstationary data set on a two dimensional area D and I'm interested in the prediction and the estimated covariance matrix. For this purpose I'm using the function btgp. As far as I understand, btgp uses a MCMC algorithm to split up D along lines parallel to the coordinate axes and estimates independent

*Hello All,* I tried figuring out the problem, I was trying to use laplacedot to predict the long term, which however would not do a good job. Then, I tried to do a point by point prediction and building the model again, everytime. It shows me better results. I tried writing my kernel function (matern covariance function), and attached is the result of that. Red lines show the fit and blue lines

> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch [SMTP:r-help-bounces at stat.math.ethz.ch] On Behalf Of ernesto > Sent: Friday, January 13, 2006 9:25 AM > To: Mailing List R > Subject: [R] help with gepRglm::likfit.glsm > > Hi, > > I'm exploring likfit.glsm and I need some help. I have to say that I'm > not an MCMC expert ... >