[R] Summer student internship placement at University of York / YCCSA / SEI (paid)

reference (not just contact details). Your application will not be accepted without this reference (max 1 page). Best, -- Corrado Topi Stockholm Environment Institute Mob: +44 (0) 7769 601784 Tel: +44 (0) 1904 322893 Skype: corrado-eeos Website: sei-international.org University of York York

Re: [R] Error singular gradient matrix at initial parameter estimates in nls

in an interpolating function rather than doing statistics, and I can understand the reluctance of statisticians to countenance such a use of nls. And Bert's comment on overparametrization is almost certainly valid also. JN -- Corrado Topi PhD Researcher Global Climate Change and Biodiversity Area

[R] Error singular gradient matrix at initial parameter estimates in nls

is 4, then the nls works perfectly (but that excludes all the k5 kn). Can anyone help me with suggestions? Thanks in advance. Alternatively, what do you suggest I should do? Shall I abandon nls in favour of optim? Regards -- Corrado Topi PhD Researcher Global Climate Change

[R] From THE R BOOK - Warning: In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!

)x -0.3603 0.4480 Degrees of Freedom: 510554 Total (i.e. Null); 510553 Residual Null Deviance: 24420 Residual Deviance: 23240AIC: 700700 Warning message: In eval(expr, envir, enclos) : non-integer #successes in a binomial glm! Regards -- Corrado Topi PhD

Re: [R] Error singular gradient matrix at initial parameter

Hi Gabor, same problem even using nls2 with method=brute-force to calculate the initial parameters. Best, Gabor Grothendieck wrote: You could try method=brute-force in the nls2 package to find starting values. On Tue, Mar 30, 2010 at 7:03 AM, Corrado ct...@york.ac.uk wrote: I am using

Re: [R] Error singular gradient matrix at initial parameter

Yes, of course. The problem still stays. Gabor Grothendieck wrote: Sorry, its algorithm=brute-force On Tue, Mar 30, 2010 at 10:29 AM, Corrado ct...@york.ac.uk wrote: Hi Gabor, same problem even using nls2 with method=brute-force to calculate the initial parameters. Best, Gabor

Re: [R] From THE R BOOK - Warning: In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!

I am certain I am not doing anything wrong. Regards -- Corrado Topi PhD Researcher Global Climate Change and Biodiversity Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904 328645, E-mail: ct...@york.ac.uk __ R

Re: [R] From THE R BOOK - Warning: In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!

at page 514/515. Rubén Roa wrote: -Mensaje original- De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] En nombre de Corrado Enviado el: martes, 30 de marzo de 2010 16:52 Para: r-help@r-project.org Asunto: [R] From THE R BOOK - Warning: In eval(expr, envir, enclos) : non

[R] Changing content of column in data.frame + efficient join extraction between 2 data.frames

) where pcs$str_1 is the same as East$str_1. I have even tried regexp, but cannot modify the content of pcs$str_1 to remove the inveretd commas and change the case to lowercase. How do I do that? Regards -- Corrado Topi PhD Researcher Global Climate Change and Biodiversity Area 18,Department

Re: [R] S4: Multiple inheritance

://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Corrado Topi PhD Researcher Global Climate Change and Biodiversity Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904 328645, E-mail: ct

Re: [R] Factor variables with GAM models

the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Corrado Topi PhD Researcher Global Climate Change and Biodiversity Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904 328645, E

Re: [R] Constrained non linear regression using ML

/R2WinBUGS, JAGS/r2jags, etc.) On Wed, Mar 17, 2010 at 10:03 AM, Corrado ct...@york.ac.uk wrote: Dear Arne, Gabor, I solved the problem with betareg (downloaded the package). I run it on my data, and unfortunately the constraint is definitively active, if I remove the active variables, I

Re: [R] Constrained non linear regression using ML

Dear Gabor, dear R users, I had already read the betareg documentation. As far as I can understand from the help, it does not allow for constrained regression. Regards Gabor Grothendieck wrote: Check out the betareg package. On Tue, Mar 16, 2010 at 2:58 PM, Corrado ct...@york.ac.uk wrote

Re: [R] Constrained non linear regression using ML

is strange because, because the function dimnames is not called anywhere. Regards Gabor Grothendieck wrote: Try it anyways -- maybe none of your constraints are active. On Wed, Mar 17, 2010 at 6:01 AM, Corrado ct...@york.ac.uk wrote: Dear Gabor, dear R users, I had already read

Re: [R] Constrained non linear regression using ML

distribution, it does not mean that your error term also follows a beta distribution. And it the distribution of the error term which is crucial for specifying the likelihood function. /Arne -- Corrado Topi PhD Researcher Global Climate Change and Biodiversity Area 18,Department of Biology University

[R] Constrained non linear regression using ML

. The nearest thing was nlrob, but even that does not allow for constraints, as far as I can understand. Any suggestion? Regards -- Corrado Topi PhD Researcher Global Climate Change and Biodiversity Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904 328645, E-mail

[R] Distance between sets of points in transformed environmental space

worth considering? Kind 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: ct...@york.ac.uk __ R-help@r-project.org mailing

Re: [R] Distance between sets of points in transformed environmental space

-- 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: ct...@york.ac.uk __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Concave hull

) ) contour(dens2, add=TRUE) # You can pick a single contour if you like: contour(dens2, level=0.05, col=red, add=TRUE) contour(dens2, level=0.10, col=blue, add=TRUE) contour(dens2, level=0.005, col=red, add=TRUE) (More bat-like.) -- Corrado Topi Global Climate Change Biodiversity Indicators Area

[R] Concave hull

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

[R] IRLS or other iteratively re weighted optimization algorithms with constraints in R

Dear list, is there an iterative re weighted least square based algorithm or any or other iteratively re weighted optimisation algorithms for non linear (and possibly non parametric) optimisation problems with constraints available in R? Regards -- Corrado

[R] Problem with dist (bug?)

? Package dist 0.4-3 on R 2.9.2 on Kubuntu 904 64 bit. 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: ct...@york.ac.uk __ R-help@r

Re: [R] Problem with dist (bug?)

of not posting it to the list We need to see the data and the script that produced the error. On Fri, Oct 2, 2009 at 5:06 AM, Corrado ct...@york.ac.uk wrote: Dear list, using package proxy. In one situation, the dissimilarity between two vectors based on method=correlation returns a value

Re: [R] Fetch large sized file from SQL

and provide commented, minimal, self-contained, reproducible code. -- 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: ct...@york.ac.uk __ R

Re: [R] Fetch large sized file from SQL

, therefore odbc can sqlQuery() can not handel it. On Fri, Oct 2, 2009 at 2:14 PM, Corrado ct...@york.ac.uk wrote: You can try using RODBC, it allows you to connect to databases using the ODBC driver. I had some difficulties using it with the postgresSQL driver in the past, because of some

[R] A point in a vector?

element j=7 to the element j+1=8. Can we do that with an R command? 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: ct...@york.ac.uk

Re: [R] Something wrong with my function Please Help

) + } tt-matrix(data=1:4,nrow=2,ncol=2) OR.RR(tt) Error in OR.RR(tt) : subscript out of bounds Many Thanks Tu -- 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: ct

Re: [R] SAS user now converting to R - Help with Transpose

/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Corrado Topi Global Climate Change Biodiversity Indicators Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone

Re: [R] define a new family (and a new link function) for gam in gam package

16:06:31 Simon Wood wrote: I am using gam in gam package (not in mgcv) it is possible to force gam in mgcv to behave like gam in gam package? -- not *exactly*, no. But what do you want to do? (i.e. what feature of `gam' do you need?) -- Corrado Topi Global Climate Change

Re: [R] define a new family (and a new link function) for gam in gam package

Dear David, I am using gam in gam package (not in mgcv) it is possible to force gam in mgcv to behave like gam in gam package? On Thursday 17 September 2009 23:00:17 David Winsemius wrote: On Sep 17, 2009, at 1:39 PM, Topi, Corrado wrote: Dear R list, is it possible to define a new

[R] define a new family (and a new link function) for gam in gam package

Dear R list, is it possible to define a new family (and a new link function) for gam in gam package? How? I read the help for gam, family, gam.model, make.link but I did not find a solution. Regards -- Corrado Topi Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone

[R] Help with error on function: Error in .... attempt to apply non-function

what I am doing wrong? kubuntu 904 64 bit, R 2.9.2 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: ct...@york.ac.uk __ R

Re: [R] Help with error on function: Error in .... attempt to apply non-function

Dear Duncan, this is a reproducible example: it is the function copied straight from my Eclipse. I found the mistake (thanks to Peter) On Tuesday 15 September 2009 11:15:29 Duncan Murdoch wrote: Corrado wrote: Dear R gurrus, I wrote this function http://scsys.co.uk:8002/33852

Re: [R] Help with error on function: Error in .... attempt to apply non-function (Solution)

Dear friends, the problem with the error has been solved (thanks to Peter). The line 41 in http://scsys.co.uk:8002/33852 should be rewritten as M-k*((x-t[i])*m0+(t[i+k]-x)*m1)/((k-1)*(t[i+k]-t[i])) On Tuesday 15 September 2009 11:32:53 Corrado wrote: -- Corrado Topi Global Climate Change

[R] identical(length(x), 1) returns FALSE, but print(length(x)) is 1, length(x)==1 is TRUE, and is.integer(lenght(x)) is TRUE????

Dear R, the condition: identical(length(x),1) returns FALSE but print(length(x)) returns 1 and: is.vector(x) is TRUE. is.integer(length(x)) is TRUE length(x) ==1 is TRUE I am puzzled. Regards -- Corrado Topi Global Climate Change Biodiversity Indicators Area 18,Department of Biology

[R] Function returns different results on the vector as a whole vs. on the values of the vector

the inner for iteration with the expression: y-mspline(i=h,x=X,k=3,t=Q) 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: ct...@york.ac.uk

Re: [R] identical(length(x), 1) returns FALSE, but print(length(x)) is 1, length(x)==1 is TRUE, and is.integer(lenght(x)) is TRUE????

any vector of length 1, for example x-1. Plus all the command that where in my previous email What is the logic behind identical(length(x),1) being false? Regards -- Corrado Topi Global Climate Change Biodiversity Indicators Area 18,Department of Biology University of York, York, YO10

[R] Intercept=0 in gam from gam package

Dear R list, is it possible to force the intercept to assume the value of 0 (that is no intercept) in gam from gam package? 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

[R] Negative AIC

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: ct

Re: [R] Negative AIC

My worry is: can I compare negative AIC with positive AIC? does the comparison still hold? On Thursday 10 September 2009 15:57:01 Ben Bolker wrote: Corrado-5 wrote: Dear R list, I just obtained a negative AIC for two models (-221.7E+4 and -230.2E+4). Is that normal? It's

Re: [R] Negative AIC

that something was wrong, if you have some AIC values that are on the order of -20,000 (as below) and others that are +20,000 ... Ben Bolker Corrado wrote: My worry is: can I compare negative AIC with positive AIC? does the comparison still hold? On Thursday 10 September 2009 15:57:01

[R] AIC and goodness of prediction - was: Re: goodness of prediction using a model (lm, glm, gam, brt,

) R^2, but I am unhappy with them. 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: ct...@york.ac.uk __ R

[R] Matrix regression

). Any suggestion? I will be more than happy to share the results of my quest with the list or developers. 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: ct

Re: [R] Strange error returned or bug in gam in mgcv????

-0.629 -0.629 -0.629 ... $ PC10 : num 0.26 0.26 0.26 0.26 0.26 ... It looks ok to me. What do you think? On Tuesday 01 September 2009 18:43:24 Gavin Simpson wrote: On Tue, 2009-09-01 at 17:55 +0100, Corrado wrote: Dear Simon, I have stored all information at the link: http://scsys.co.uk

[R] goodness of prediction using a model (lm, glm, gam, brt, regression tree .... )

=A for training Y~T,data=B for predicting I have devised a couple of method based around 1) standard deviation 2) R^2, but I am unhappy with them. Regards -- Corrado Topi Global Climate Change Biodiversity Indicators Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904

Re: [R] Google's R Style Guide

prescriptive which makes our research group life easier? 2) What class type should I use? From what you two say, I should use S3 because is easier to use what are the disadvantages? Is there an advantages / disadvantages table for S3 and S4 classes? Thanks -- Corrado Topi Global Climate

[R] Strange error returned or bug in gam in mgcv????

/FALSE needed Calls: gam ... smooth.construct - smooth.construct.tp.smooth.spec - array In addition: Warning message: In array(0, n * k) : NAs introduced by coercion Execution halted Thanks in advance, Best regards -- Corrado Topi Global Climate Change Biodiversity Indicators Area 18,Department

Re: [R] Strange error returned or bug in gam in mgcv????

Nope Of course, it was just a copy and paste problem On Tuesday 01 September 2009 15:00:34 David Winsemius wrote: On Sep 1, 2009, at 9:51 AM, Corrado wrote: Dear friends, what is this error message in gam I cannot understand what it means is it a bug

Re: [R] Strange error returned or bug in gam in mgcv???? - additional information

Here I pasted the code from when I opened the R shell, so that it possible to see what is going on: http://scsys.co.uk:8002/33309?hl=onsubmit=Format+it! Thanks in advance -- Corrado Topi Global Climate Change Biodiversity Indicators Area 18,Department of Biology University of York, York

Re: [R] Strange error returned or bug in gam in mgcv???? - yet more additional information

I am using mgcv 1.4-1.1 on Fedora 9 64 bit on an Opteron server with 8Gb of RAM. On Tuesday 01 September 2009 15:19:28 Corrado wrote: Here I pasted the code from when I opened the R shell, so that it possible to see what is going on: http://scsys.co.uk:8002/33309?hl=onsubmit=Format

Re: [R] Strange error returned or bug in gam in mgcv????

trying to fit, and I'll try and track down exactly where it's failing, and put in a more informative message. best, Simon On Tuesday 01 September 2009 14:51, Corrado wrote: Dear friends, what is this error message in gam I cannot understand what it means is it a bug

Re: [R] Google's R Style Guide

all the hints you need. Hadley -- 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: ct...@york.ac.uk __ R-help@r-project.org mailing

Re: [R] Best R text editors?

interface (e.g. Tinn-R, Komodo+SciViews). Thanks! --j -- 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: ct...@york.ac.uk __ R-help@r

Re: [R] Best R text editors?

I am using 3.4.2 not 3.5, I would not know. But it is worth visiting the STATET mailing list archive and subscribe to the mailing list. On Friday 28 August 2009 09:22:14 [Ricardo Rodriguez] Your XEN ICT Team wrote: Hi! Corrado wrote: Eclipse + StatET (the R plugin) both on Linux

Re: [R] [Rd] Formulas in gam function of mgcv package

ANOVA. best, Simon -- 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: ct...@york.ac.uk __ R-help@r-project.org mailing list https

Re: [R] [Rd] Formulas in gam function of mgcv package

this should have been directed - redirected there now] On Mon, 2009-08-24 at 17:02 +0100, Corrado wrote: 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

Re: [R] [Rd] Formulas in gam function of mgcv package

: [Note R-Devel is the wrong list for such questions. R-Help is where this should have been directed - redirected there now] On Mon, 2009-08-24 at 17:02 +0100, Corrado wrote: Dear R-experts, I have a question on the formulas used in the gam function of the mgcv package. I am trying

[R] Distance between clusters

(dissimilarityMatrix) hctr-hclust(distPredTurn,ward) cutree(hctr,k=48) 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: ct...@york.ac.uk

Re: [R] Using very large matrix

. Best, On Thursday 26 February 2009 15:43:11 Jay Emerson wrote: Corrado, Package bigmemory has undergone a major re-engineering and will be available soon (available now in Beta version upon request). The version currently on CRAN is probably of limited use unless you're in Linux

[R] Using very large matrix

Any suggestions? I have searched the mailing list, but to no avail. 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: ct...@york.ac.uk

[R] Similarity between clusters generated by hclust + cutree

, for example as (1-distance or dissimilarity) between the centroids. Any clue on how I can do that? 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: ct...@york.ac.uk

Re: [R] Principal Component Analysis - Selecting compo nents? + right choice?

-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of S Ellison Sent: Thursday, December 11, 2008 9:37 AM To: r-help@r-project.org; Corrado Subject: Re: [R] Principal Component Analysis - Selecting components? + right choice? If you're intending to create a model using PCs

[R] Principal Component Analysis - Selecting components? + right choice?

would have been better off by using something different from PCA? 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: [EMAIL PROTECTED

[R] *** buffer overflow detected ***: /usr/lib64/R/bin/exec/R terminated on R 2.6.2 to 2.8.0: logging a bug?

/LC_MAborted [EMAIL PROTECTED]:~$ OS: Mandriva 2008.1 x86_64 Postgresql: 8.3.1 (PostGIS enabled) R: from 2.6.2 from repository to 2.8.0 repackaged Is it my doing, or R's doing? Best, -- Corrado Topi Global Climate Change Biodiversity Indicators Area 18

[R] R base 2.7.2 packaged for Mandriva 2008.1 x86_64: anyone interested?

I have packaged R base 2.7.2 for Mandriva 2008.1 x86_64 who should I send it to so that it can be made available to everybody? It is my first attempt and it works well on my computer, but it will need some testing. Best, -- Corrado Topi Global Climate Change Biodiversity Indicators

[R] installing source package RGtk2 on mac os x 10.4.11

/ RGtk2/libs/i386/*: No such file or directory ERROR: compilation failed for package 'RGtk2' ** Removing '/Library/Frameworks/R.framework/Versions/2.7/Resources/ library/RGtk2' CAN ANYONE HELP IN RESOLVING THIS ERROR? Thank you very much for your support. Corrado Giannasca