search for: topy

Dear R-lings, I did not know which list to post to, because it is a studentship so not really a job, so it did not fit the r-sig-jobs list .... and it is about devloping an extension package interfaced with R .... I hope I did not upset anyone. If so apologies. The Centre For Complex systems Analysis at the University of York (YCCSA) in UK in collaboration with Stockholm Environment

I am using nls to fit a non linear function to some data. The non linear function is: y= 1- exp(-(k0+k1*p1+ .... + kn*pn)) I have chosen algorithm "port", with lower boundary is 0 for all of the ki parameters, and I have tried many start values for the parameters ki (including generating them at random). If I fit the non linear function to the same data using an external

I'm trying to diagnose a bizarre problem in which the help files invoked from R are partially defective. help(), for instance, is missing the Description, Arguments and See Also sections; as are help files for other commands (see below). I've built 2.8.1, 2.10.1 and HEAD with the same results; do I need to set an environment or configuration variable? help

Dear friends, what is this error message in gam???? I cannot understand what it means .... is it a bug? gam_bray_scot24_pc_0505<gam(bray~s(PC1,PC2,PC3,PC4,PC5, PC1.1,PC2.1,PC3.1,PC4.1,PC5.1),data=dist_scot24_vector_with_climate) Error in if (length(data) != vl) { : missing value where TRUE/FALSE needed Calls: gam ... smooth.construct -> smooth.construct.tp.smooth.spec -> array In

...or(rep(c("Placebo", "Not Placebo"), each=2500), levels=c("Placebo", "Not Placebo")) xdens <- tapply(x, trtp.f, density, na.rm=TRUE) ydens <- tapply(y, trtp.f, density, na.rm=TRUE) topx <- max(unlist(lapply(xdens, function(x) max(x$y)))) topy <- max(unlist(lapply(ydens, function(x) max(x$y)))) xrange <- range(x, na.rm=TRUE) yrange <- range(y, na.rm=TRUE) nf <- layout(matrix(c(2,0,1,3),2,2,byrow=TRUE), c(3,1), c(1,3), TRUE) par(mar=c(4,4,0,0), oma=c(0,0,0,0)) getcolors <- rainbow(length(levels(trtp.f)), al...

Dear list, I have a strange requirement .... I have a vector, for example v<- c(0,0,0,0,1,2,4,6,8,8,8,8). I have a value,for example x<- 4.8. I would like to understand in which sub interval of v is x. In this case, v would be in the sub interval [4,6] that is in the subinterval starting from element j=7 to the element j+1=8. Can we do that with an R command? Regards -- Corrado Topi

Dear All, The Decision and Policy Analysis (DAPA) program of the International Center for Tropical Agriculture (CIAT), a not-for-profit agriculture, biodiversity and ecosystems research center, supported by the Consultative Group of International Agricultural Research (CGIAR), with headquarters in Cali, Colombia is looking for a visiting researcher to work on various modeling projects, primarily

Dear friends, Do you know how to calculate the CONCAVE hull of a set of points (2- dimensional or n-dimensional)? is that possible in R? (With a "smoothing" parameter of course). Best, -- Corrado Topi Global Climate Change & Biodiversity Indicators Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904 328645, E-mail: ct529 at york.ac.uk

Dear friends, I am testing glm as at page 514/515 of THE R BOOK by M.Crawley, that is on proportion data. I use glm(y~x1+,family=binomial) y is a proportion in (0,1), and x is a real number. I get the error: In eval(expr, envir, enclos) : non-integer #successes in a binomial glm! But that is exactly what was suggested in the book, where there is no mention of a similar warning. Where am I

Dear R-friends, How do you test the goodness of prediction of a model, when you predict on a set of data DIFFERENT from the training set? I explain myself: you train your model M (e.g. glm,gam,regression tree, brt) on a set of data A with a response variable Y. You then predict the value of that same response variable Y on a different set of data B (e.g. predict.glm, predict.gam and so on).

Hi all, I started to ask this in the other post, but it is off topis...so here it is again. I have a data.frame (created with the helpof this mail list) that looks like this : 'data.frame': 22801 obs. of 15 variables: $ V1 : chr "HUMUS" "SLABO" "MALO" "SLABO" ... $ V2 : chr "IN" "GRANULIRAN"

Dear friends, I have to use a very large matrix. Something of the sort of matrix(80000,80000,n) .... where n is something numeric of the sort 0.xxxxxx I have not found a way of doing it. I keep getting the error Error in matrix(nrow = 80000, ncol = 80000, 0.2) : too many elements specified Any suggestions? I have searched the mailing list, but to no avail. Best, -- Corrado Topi Global

Dear R list, I just obtained a negative AIC for two models (-221.7E+4 and -230.2E+4). Is that normal? Regards -- Corrado Topi Global Climate Change & Biodiversity Indicators Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904 328645, E-mail: ct529 at york.ac.uk

Hi, I've been planning to develop a support for tunneling between "local_tcp => server_AF_UNIX". This way, every user of server machine, can have: 1. personal address space (if socket is located on personal directory). Currently one must check assigned local port every time starting a server (e.g. vncserver), and redirect a local port to "random" remote port. 2. Added

How do people feel about eliminating the list serv in favor of the forum? Or limiting the list serv to news? I have one vendor I work with that has a knowledgebase, a forum, and two mailing lists - with much overlap of topics. It is truly enough to drive you totally crazy. I'm not suggesting anything. Just bringing it up. Geoff

Dear R gurus, I have some climatic data for a region of the world. They are monthly averages 1950 -2000 of precipitation (12 months), minimum temperature (12 months), maximum temperature (12 months). I have scaled them to 2 km x 2km cells, and I have around 75,000 cells. I need to feed them into a statistical model as co-variates, to use them to predict a response variable. The climatic

Dear R-experts, I have a question on the formulas used in the gam function of the mgcv package. I am trying to understand the relationships between: y~s(x1)+s(x2)+s(x3)+s(x4) and y~s(x1,x2,x3,x4) Does the latter contain the former? what about the smoothers of all interaction terms? I have (tried to) read the manual pages of gam, formula.gam, smooth.terms, linear.functional.terms but

Dear R-experts, I have a question on the formulas used in the gam function of the mgcv package. I am trying to understand the relationships between: y~s(x1)+s(x2)+s(x3)+s(x4) and y~s(x1,x2,x3,x4) Does the latter contain the former? what about the smoothers of all interaction terms? I have (tried to) read the manual pages of gam, formula.gam, smooth.terms, linear.functional.terms but

Being very old school, I think mailing lists rock :-) Maybe we could have a web interface to the mailing list instead of a forum? :-D Kelv Eric Robinson wrote: > i wouldnt say eliminate the list serv > however i would limit it to newsletters > or other specific topics which you would NOT need to > look for in the new forums > > I would also say that the majority of issues

Hi, Just want to ask if ocfs2-2.6.9-89.0.11.ELsmp-1.2.9-1.el4.x86_64.rpm compatible on kernel release 2.6.9.-89ELsmp? Regards, Francis Delos Reyes Oracle Apps Engineer PS & MIS Division MSI-ECS Phils. Inc. Topy II Bldg., #3 Economia St., Libis, Quezon City 1110 Website : www.msi-ecs.com.ph Trunkline : +63 (2) 688-3333 / 439-7777 Direct Line : +63 (2) 688-3505 Mobile No : +63 (916) 713-3135 Email : fdelosreyes at msi-ecs.com.ph -------------- next part -------------- An HTML...