search for: ruwpa

...er functions/code that will do this. - you can do this using par(new=TRUE) between calls to persp() to get plots overlayed: persp.gam() in package mgcv does this, for example. ______________________________________________________________________ > Simon Wood snw at st-and.ac.uk http://www.ruwpa.st-and.ac.uk/simon.html > The Mathematical Institute, North Haugh, St. Andrews, Fife KY16 9SS UK > Direct telephone: (0)1334 463799 Indirect fax: (0)1334 463748 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci...

Would anyone be able to provide insight for the following question, please? Setting: estimation of prediction intervals for age-period-cohort models using GAMs (rate ~ s(age,period)) Method: bootstrap (Davison and Hinkley, 1997) Issue: standardisation of the residuals for resampling requires an adjustment using the diagonals of the hat matrix. Is there a simple way to get the hat values out of a

I hope this is the right place for my question, my apologise if it is not. I am a long time Mac/UNIX user and my local machine is a Mac. At present I run an orphan version of R which a worried out of (*****). On reflection I guess I'd better not say, but it changed my life. It is howerver a little out of date ( say 0.6ish) So the obvious question, do I move to MacLinux or is there likely

Hello Is there a simple way of fitting a 'smooth curve' to a time series with the additional constraint that the first derivative be not negative? I know that if I choose a "large enough" window for a moving average sort of filter on my data I end up with a non-decreasing curve anyway, but I'm wondering if there is a pre-built method that can include the constraint from

To R users: I am running an R program with .Fortran() calls, and discovered that the memory, as reported by gc(), continues to grow with new loop iterations. I hadn't been able to find information on how to 'destroy' objects in earlier loops to keep the program size from continuous growth; calling rm() to remove some what I thought to be large arrays didn't seem to help. Any

...s. First of all E. Paradis pointed out, that he would prefer to see R presented in a different way, please read his whole reply https:// www.stat.math.ethz.ch/pipermail/r-help/2003-December/041951.html. Packages which has been recommended are: WISP http://www.ruwpa.st-and.ac.uk/estimating.abundance/WiSP/ ape www.stat.uni-muenchen.de/~strimmer/ publications/ape2003.pdf phylogr http://www.stat.ucl.ac.be/ISdidactique/Rhelp/library/PHYLOGR/html/00Index.html genetics cran.r-project.org/doc/packages/genetics.pdf VEG...

Hello, I'm a final-year student in statistics and I deal with as subject of thesis "rare and elusive populations". I want make a simulation, in which I want generate a population and estimate its density with "line transect method" using R language. Is there someone that can send me the program that generate a population and estimate its density with "line transect

Hi Does R have a proceedure/software for capture recapture? Thank you. Lawrence Lessner Please respond to LLessner at nycap.rr.com

I noticed that there is no generalised additive model functions in R (1.1.1) ... is there a package that implements them? Thanks Prasad ***************************************************************** Mr. Anantha Prasad, Ecologist/GIS Specialist USDA Forest Service, 359 Main Rd. Delaware OHIO 43015 USA Ph: 740-368-0103 Email: aprasad at fs.fed.us Web:

Hi Allstaters, Does anybody know if the R package can fit Generalised Additive Models?. Thanks indeed, Aurelio. atobias at ole.com ___________________________________________________________________ Consigue tu e-mail gratuito TERRA.ES Haz click en http://www.terra.es/correo/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Dear R-users, I've a question regardiing multiple linear regression. Is it possible to regress a function of the type y=b0+b1*x1+b2*x2 constraining the partial derivatives dy/dx1 and dy/dx2 to be greater than 0 ??? Thank you very much for any possible suggestion. Best regards Marco -- Marco Pagani - PhD Student Dipartimento di Scienze della Terra - Sezione Geofisica Universita' degli

Hi, I need help regarding estimating a linear model where restrictions are imposed on the coefficients. An example is as follows: Y_{t+2}=a1Y_{t+1} + a2 Y_t + b x_t + e_t restriction a1+ a2 =1 Is there a function or a package that can estimate the coefficient of a model like this? I want to estimate the coefficients rather than test them. Thank you for your help Ahmad Abu Hammour --------------

Hello R-user, sorry for this very off-topic question. But I shall present R to my dept. (pro's and con's and what it can do). The pro's and con's are easy but not what R can do (additional to the "normal" statistics). I looked through the packages, but the enormous amount of packages makes it very difficult for me to decide which one is worth mentioning. I used

Hello everybody, I think my question has been asked before but I am posing it once again since I need it. Is there any way to find the rank of a matrix in R or Splus? ______________________ Indrajit SenGupta Department Of Statistics St. Xavier's College Calcutta University mailto:indra_calisto at yahoo.com mailto:indrajitsg at vsnl.net

I'm trying to figure out a good way to smooth binary data. The ideal approach appears to be the "sm.logit" function in library "sm", but I haven't had success with it. Below, is some code illustrating what I've come up with so far, but I'm hoping there is a better approach. I'm using R 1.2 (development) under Windows 98. -Bill library(MASS) data(birthwt)

Hola! I have some problems with gam in mgcv. Firts a detail: it would be nice igf gam would accept an na.action argument, but that not the main point. I want to have a smooth term for time over a year, the same pattern repeating in succesive years. It would be natural then to impose the condition s(0)=s(12). Is this possible within mgcv? I tried to obtain this with trigonometric terms, aca:

Hi, Is there any difference in the gam() function in R and Splus? I just tried to do an Ace Plot (plotting the result of the gam() ) on a dataset. R gave me four (the data set has one response and four explanatory variables) ace plot with perfectly straight lines. I kind of felt that it is highly unlikely, so I tried it in Splus. Using the exact commands Splus gave me four completely

Hi helpteam, may you send me the exact syntax for getting the estimated degrees of freedom for a gam? I tried to do it by name <- edf(gam_name) but this didn't work. Thanks a lot, greetings from Munich, Thomas Nittner -- Thomas Nittner Tel: ++49 89 2180 3197 Akademiestra?e 1 Fax: ++49 89 2180 5042 D - 80799 M?nchen http://www.stat.uni-muenchen.de

Has anyone got a function for smooth monotonic interpolation of a univariate function? I'm after something like the NAG function PCHIM which does monotonic Hermite interpolation. Alternatively, montononic cubic spline interpolation. Please reply directly. Rob Hyndman ___________________________________________________ Rob J Hyndman Associate Professor & Director of Consulting

Hi, Is it possible to compare two GAM objects created with the gam() function from the mgcv package. I use a slightly modified version of anova.glm() named anova.gam(), modified from John Fox (2002). It often gives me some aberant responses, especially with "F" test. I use a quasibinomial model and scale (dispersion) is calculated and used in the calculation of the F value. Does someone