search for: maternal

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 do things with them - like examine their distributions. >wg[[1]] gives Ca...

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 fixed effects: DISTANCE - distance to maternal tree; two levels 'CLOSE' or 'AWAY' (four groups of seedlings each per tree), and PLATEAU -...

...tors 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 thanks, michael foster ________________________________________________ E. Michael Foster (W) 814-865-1923 (Fax) 814-863-0000 After 7/1 Work: Professor of Maternal and Child Health School of Public Health University of North Carolina, Chapel Hill Rosenau Hall, CB# 7445 Chapel Hill, NC 27599-7445 Home: 309 Old Larkspur Way Chapel Hill, NC 27516 UNC School of Public Health--voted #2 SPH by . US News . Association of Hand-to-Fin Carp Fisher...

...of the vertical extent of the graphics device. My graphics are correct horizontally, but are confined to about one-fifth of the vertical extent of the graphics device. As a result they look odd, and are almost usless for analysis. They do seem, however, to be the correct graphics (In this case maternal weight against gestational age for two groups of women). I have read the manual for the lattice commands, and the nlme command, and indeed for the plain plot command, but failed to get any ideas. Roughly, I need to tell lattice to use the full vertical extent of the display device. If it helps,...

...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 advance of the loop and fill them in, but I'd rather not. Is there a way to use an index of a loop in this manner? thanks, michael E. Michael Foster Professor of Maternal and Child Health School of Public Health University of North Carolina

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

...might guess, what I want to do at the end is rbind the little files, and the lack of consistent column names causes the program to choke. Suggestions? /m -- ________________________________________________ E. Michael Foster download vcard: www.unc.edu/~emfoster/Foster.vcf Professor of Maternal and Child Health of Health Policy and Administration Associate Chair for Faculty Advancement (W) 919-966-3773 (F) 919-966-0458 Office: 407C Rosenau Hall << use for FEDEX School of Public Health University of North Carolina, Chapel Hill Rosenau Hall, Campus Box# 7445 Chapel Hill, NC 27...

...tive-definite and will look something like the following: 1 0 0 0.5 0 0 1 0 0.5 0 0 0 1 0 0 0.5 0.5 0 1 0.25 0 0 0 0.25 1 The elements basically correspond to the proportion of genes shared by any two individuals. The second matrix will model additive maternal effects (Vm) and the variance-covariance matrix associated with these effects will usually not be positive definite as shown below. 1 1 0 0 0.5 1 1 0 0 0 0 0 1 1 0 0 0 1 1 0 0.5 0 0 0 1 The elements here correspond to the proportion of genes s...

...models Packages and graphics functions Distributions and simulations Developing your own functions Accessing the operating system The R API Database connectivity Programming in R using Fortran and C subroutines Animal breeding applications - Inbreeding routines - Animal models - Iteration on data - Maternal genetic effects - Multiple traits - Random regression - Threshold models - Genomic applications On the application form you will be able to indicate particular problems that you would like to see covered, and we will try to work them into the course. -- Lep pozdrav / With regards, Gregor Go...

...add=F, type= "l", lwd=2, xlab="Ocassion of Measurement", ylab="y", ) I had the same issue with a trellis plot: obs.plot<- xyplot(msen~Age, data=obs.msen, groups=ID2, type="l", lwd=2, col=1, xlab="Age", ylab="Maternal Sensitivity", ylim=c(-.5,4.5),xlim=c(0,5), panel=panel.superpose, par.settings=list(background=list(col="white")) ) obs.plot ----- ------------------------------- Project Scientist University of California, Irvine -- View this message in context: http://www.nabble.com/...

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

...get seems to suggest it is finding the files empty after filtering. However, all of the files are primarily words with only a few numbers mixed in. Any help appreciated. The data I am using is the MEDLINE data set and the first file in the data set med.000001 looks like this: correlation between maternal and fetal plasma levels of glucose and free fatty acids . correlation coefficients have been determined between the the command I am using looks like this, with the resulting error below: > > dtm <- textmatrix(LSAwork, stemming=TRUE, stopwords=StopListm, minGlobFreq=1, minWordLength=2, r...

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 ... >