similar to: Interfacing between Optimization Software and R

Dear All I am running Sweave with xtable and want to put the caption placement on top. But this does not work. Any idea what is going wrong? Here is an example that runs properly with the exception of the caption placement in the pdf-file. \documentclass[11pt,a4paper]{article} \usepackage{Sweave} \begin{document} <<>>= x = runif(100, 1, 10) y = 2 + 3 * x + rnorm(100)

Hi I am trying to run a R-script in the windows command window: e.g. rscript example.R ... It plots some stuff, but using par(ask=T) doesn't let the plotting pause. Any idea how to do this? Thanks Renger __________________________________________________ Renger van Nieuwkoop Goldiwilstrasse 16 F / 3600 Thun Tel +41 33 221 53 05 / +41 79 818 53 73 mailto: renger at vannieuwkoop.ch

Hi When I have a plot under Windows and want to save it as a ps or jpg using right click, the window for choosing the directory opens. But if I want to choose a directory R freezes. This happens on my XP notebook as well on my P x64 machine. Renger _______________________ Renger van Nieuwkoop=20 Ecoplan - Thunstrasse 22=20 3005 Bern Switzerland=20 Tel. +41 31 356 61 61 Mail: renger at

Hi I imported a lot of text-files with R and saved them by using: save(list=ls(),file="Grunddaten.Rdata") After that I wanted to check my saved data and wrote source("Grunddaten.Rdata") This gives me an error: Error in parse(file, n = -1, NULL, "?") : Syntax error at 1: Can somebody tell me what I am doing wrong (I am using R.2.4.1)? Renger

Dear R-helpers, I have encountered the following unexpected behaviour of R-2.10.1, but not R-2.9.0, on both RHEL 4 and Ubuntu Karmic (precompiled via synaptic or built from source). I have a character vector from which I want to extract a certain pattern that is surrounded by junk as in: > nn <- sprintf("junk_%02d_junk", 1:2) > nn [1] "junk_01_junk"

At office I have been introduced by another company to new, complex energy forecasting models using gams as the basic software. I have been told by the company offering the models that gams is specialised in dealing with huge, hevy-weight linear and non-linear modelling (see an example in http://www.gams.com/modtype/index.htm) and they say it is almost the only option for doing it. I would

hi! I hope all of you can help me this question for example GAMs: ozonea<-gam(newozone~ pressure+maxtemp+s(avetemp,bs="cr")+s(ratio,bs="cr"),family=gaussian (link=log),groupA,methods=REML) formula(ozonea) newozone ~ pressure + maxtemp + s(avetemp, bs = "cr") + s(ratio,bs = "cr") #formula of gams coef(ozonea) # extract the coefficient of GAMs

I have 5 GAMs ( model1, model2, model3, model4 and model5) Before I use some data X(predictor -January to June data) to form a equation and calculate the expected value of Y (predictand -January to June). After variable selection, GAMs (Model 1)were bulit up! R-square :0.40 NOW, I want to use new X'( predictor -July - December data) and put into Model 1, then get the expected value of Y'

Hello, My name is Dontrece Smith. I am creating figures for my GAMs. I change my qualitative variables to 1 or 2 in my dataset, so I would be able to run my GAMs. However, R will only display plots for my quantitative variables and not my qualitative variables. Is there any way to fix this issue? I listed some of my code below: > library(mgcv) This is mgcv 1.7-13. For overview type

This is my first time running GAMs in R. My csv file has these column headings: "X" "Y" "Sound" "Atlantic" "Blacktip" "Bonnet" "Bull" "Finetooth" "Lemon" "Scalloped" "Sandbar" "Spinner" "Abundance" "Diversity"

Hi everybody, I am working with gams and I have found some questions when plotting gams models. I am using mgcv, and my model looks something like this: model<- gam(x ~ s(lat,long)) I can plot the output of the model using plot(model) or plot.gam(model) and I get a surface plot. That is ok, but what I want to do now is to extract the data used to perform the surface plot. Like that I

Hi I was wondering if anyone knew how to work out the number of knots that should be applied to each variable when using gams in the mgcv library? Any help or references would be much appreciated. Thanks Kathryn Baldwin

I'm about to write my thesis in economics and will need to setup and solve a system of non-linear equations. At our university we usually use GAMS for this, and though GAMS is a fine program, it bugs me a that I won't be able to run my code after I finish my thesis without buying a license for the program(about $3.500 :-(( ) So I've looked around for NL-stuff for R, but I

Hi, I have a variable (current speed direction) which is circular (0=360 degrees), and I'd like my GLM to include the variable as a circular variable. Can I do this? And what is the code? I'm actually doing a GLM-GEE using the 'geepack' package, so want to use it in that, but also interested in whether it can also be used in GLMs and GAMs (I use the 'mgcv' package for

Dear r-helpers, I just read in an article by Virtanen et al. (2006) where vegetation-environment relationships are studied by fitting smoothed surfaces on an NMDS ordination using GAMs (Wood 2000). The authors describe, that they used R? as goodness-of-fit statistic, which they compare to the R? of fitted vectors. Calculations were carried out using the package vegan (Oksanen). I know that I can

I would like to see R packages arranged by topic. CRAN Task Views are at http://lib.stat.cmu.edu/R/CRAN/src/contrib/Views/ , but I'd like something more detailed. For example, the IMSL Fortran library, version 4 is easy to navigate and has procedures arranged according to following topics: Basic Statistics Regression Correlation Analysis of Variance Categorical and Discrete

Hi, I've been using gamm4 to build GAMMs for exploring environmental influences on genetic ancestry. Things have gone well and I have 2 very straightforward questions: 1. I've used method=REML. Am I correct that this is an alternative method for estimating the smooth functions in GAMMs rather than GCV that is often used for GAMs? I've read up on REML and it makes sense, but I'm

Hi, I am using GAMs (package mgcv) to smooth event rates in a penalized regression setting and I was wondering if/how one can select the order of the derivative penalty. For my particular problem the order of the penalty (parameter "m" inside the "s" terms of the formula argument) appears to have a larger effect on the AIC/deviance of the estimated model than the

Dear R Users, I am looking for some guidance on setting up an optimisation in R with non-linear constraints. Here is my simple problem: - I have a function h(inputs) whose value I would like to maximise - the 'inputs' are subject to lower and upper bounds - however, I have some further constraints: I would like to constrain the values for two other separate function f(inputs) and

Greetings This is a long email. I'm struggling with a data set comprising 2,278 hydroacoustic estimates of fish biomass density made along line transects in two lakes (lakes Michigan and Huron, three years in each lake). The data represent lakewide surveys in each year and each data point represents the estimate for a horizontal interval 1 km in length. I'm interested in comparing