similar to: Conditional Sequential Gaussian Simulation

Hello, I'm developing an ordination using metaMDS (package vegan). The analysis identifies 3 distinct groups that I'd like to define by either adding circles or ellipses to help identify the groups. The data set is a spatial temporal data base depicting change in each of 4 areas over 3 time periods. I can add lines that would link each of the time - x - area trajectories, but

Hello everyone, I'm working with R 2.4.1 on a PC running with XP. Trying to run isoMDS as follows: Gquad.mat <- Gquads[4:10] # extracts only the metric data variables Gquad.dist <- dist(Gquad.mat) Gquad.mds <- isoMDS(Gquad.dist) plot(Gquad.mds$points, type = "n") text(Gquad.mds$points, labels = as.character(Gquads$Quadname), cex = .5) The data.frame

*Hello all ** I am a newbie to R plotting maps. I am trying to plot over a world map a layer of Biogeochemical provinces (BGCP) by A.R. Longhurst. Each ocean region unfortunately are quite irregular in shape (not perfect squares). In GIS this layer of ocean provinces would be a layer of polygons, which I am assuming it cannot be plotted with R. I was wondering If anybody has encounter this

Dear R-SIG-Debian Team, r-devel sent me to you. Situation: Ubuntu 16.04, R-3.3.1, and a file in a folder mounted with sftp via gnome virtual file system. for(i in 1:10){cat(i,file = "/run/user/1001/gvfs/sftp:host=server/path on server/test2",append = TRUE)} gives a file that reads |10| and not |12345678910| Further checking: * A mapped network drive on Windows 7 and

Dear all, I have a question regarding the output of the coxph function. What does the 'max possible' exactly mean in the output below? Many thanks. coef exp(coef) se(coef) robust se z Pr(>|z|) smocc_zyban -0.4384 0.6451 0.8667 0.9473 -0.463 0.644 self 1.1857 3.2728 0.1405 0.1443 8.216 2.22e-16 *** smocc_zyban:self

Hello, One of our admins called me some time ago and told me, that my samba server (server-A) is connecting his server (server-B), although it shouldn't. Authentication of my samba shares is working via winbind, configured to use "server-DC1" : winbind separator = / winbind uid = 10000-20000 winbind gid = 10000-20000 winbind enum users = yes

Dear R users, we are pleased to announce the release of the 'bigleaf' R package. This package is useful for the calculation of physical (e.g. aerodynamic conductance, surface temperature), and physiological (e.g. canopy conductance, intrinsic water-use efficiency) ecosystem properties from eddy covariance data. All calculations are based on a 'big-leaf' representation of the

Hello, I'm using tinc for a few years now and that this is my first post to this list tells a lot about the quality of tinc, I think. It's just great, thanks a lot for your fine work! Recently I've discovered these messages in my syslog: Dec 28 16:34:06 vdr tinc.grue[9731]: Possible intruder ix (192.168.178.25 port 60882): wrong keylength Dec 28 16:34:06 vdr tinc.grue[9731]: Error

[Just a gentle reminder that the deadline is only one week away] ======================================================================== Call for Papers 7th IEEE International Workshop on Storage Network Architecture and Parallel I/Os (SNAPI 2011) http://snapi2011.cis.fiu.edu

[Just a gentle reminder that the deadline is only one week away] ======================================================================== Call for Papers 7th IEEE International Workshop on Storage Network Architecture and Parallel I/Os (SNAPI 2011) http://snapi2011.cis.fiu.edu

Dear Knowledgeable R Community Members, Please excuse my ignorance, I apologize in advance if this is an easy question, but I am a bit stumped and could use a little guidance. I have a finite mixture modeling problem -- for example, a 2-component gaussian mixture -- where the components have a large overlap, and I am trying to use the "mclust" package to solve this problem. I need

(Yes, today is the "John has an idea day") So, in my now fairly extensive coding of extension priority lists, I'm getting very bored of re-numbering my priority lists every time I add something at the top of the list. If I have a 7 step priority list, and I need to add something in at priority 2, then I have to re-number five other priorities, and if my Dial statement is

Dear R Helpers, I have searched for any R package or code for simulating multivariate fractional Brownian motion (mFBM) or multivariate fractional Gaussian noise (mFGN) when a covariance matrix are given. Unfortunately, I could not find such a package or code. Can you suggest any solution for multivariate FBM and FGN simulation? Thank you for your help. Best Regards, Ryan ----- Wonsang You

[This email is either empty or too large to be displayed at this time]

Hi All, I am using em() function to estimate a poisson-gaussian process from a univariate one dimension time series, but not sure how to do. In the help manual, it specify that in "pro" of the argument "parameter", if the model includes a Poisson term for noise, there should be one more mixing proportion than the number of Gaussian components. But in the example, the parameter

Hi All, I am trying to use A Gaussian quadrature over the interval (-infty,infty) with weighting function W(x)=exp(-(x-mu)^2/sigma) to estimate an integral. Is there a way to do it in R? Is there a function already implemented which uses such weighting function. I have been searching in the statmode package and I found the function "gauss.quad(100, kind="hermite")" which uses

Hello, I'm Emanuela, I'm implemented a survival analysis and I'm trying to use a frailty model with inverse gaussian distribution, but I'm not able to find the right code, because it seems to be only a gamma and a gaussian distribution. Is there also the inverse gaussian distribution? Thanks a lot Emanuela Rossi

Hi, As far as I know, there exist some programs via the function INLA, but I'm so curious if there is a specific function directly used to fit the log Gaussian Cox process model and predict the latent Gaussian field. That is, if I have a data points, then I input it in the function and don't need to revise the program. Thanks for your help. Joseph -- View this message in context:

Dear all, I am trying to fit a double gaussian to some data using the mixdist package: --- begin code --- ? library(mixdist) ? time <- seq(673,723) ??counts <- c(3,12,8,12,18,24,39,48,64,88,101,132,198,253,331,419,563,781,1134,1423,1842,2505,374,6099,9343,13009, ? ? ? ? ? ? ? ? ? ?15097,13712,9969,6785,4742,3626,3794,4737,5494,5656,4806,3474,2165,1290,799,431,213,137,66,57,41,35,27,27,27) ?

Hello, I am using the parameterized EM- Algorithm in clustering context. I use generated datasets, containing 5000 points in 20 dimensions with 5 clusters. The following poses a problem for me: emEst<-me(modelName="EEV",data=data[,1:numDimensions],z=unmap(x)) emEst<-me(modelName="VEV",data=data[,1:numDimensions],z=unmap(x)) coordProj(main="Model