[R] How can i add a color bar with base package

2006-06-09 Thread Yves Magliulo
with contributed package. it's sounds very easy but i can't get it! :-( thanks by advance. --- Yves Magliulo RD Engineer CLIMPACT : www.climpact.com tel : 01 44 27 34 31 [[alternative HTML version deleted]] __ R-help

Re: [R] Justifying R to anti open-source management

2006-05-17 Thread Yves Magliulo
hi, * http://www.r-project.org/ R-home * http://www.itp.phys.ethz.ch/econophysics/R/index.html R-financial * http://wiki.r-project.org/rwiki/doku.php#r_wiki_-_overview R-overview * http://addictedtor.free.fr/graphiques/index.php R-graphics should give you

Re: [R] Fitting model with varying number of predictors

2005-11-18 Thread Yves Magliulo
try : my_formula = as.formula(paste(y~,xxx)) lm(my_formula) note that you can play with substr in paste to change your formula paste(y~,substr(xxx,1,5)) [1] y~ x1+x2 paste(y~,substr(xxx,4,8)) [1] y~ x2+x3 furthermore, with 10 predictors x1 x2 ...x10 for instence, you can create 2 array of

[R] adding sequence for each value in a vector

2005-10-27 Thread Yves Magliulo
hi, i have a vector like : x-c(1,15,30,45,60,90,115) i know that step by step i have always more than 10 min(diff(x)) =11 i want to add for each value a sequence of value:value+9 result should be : 1 2 3 4 5 6 7 8 9 10 15 16 17 18 19 20 21 22 23 24 30 31 (...) 39 45 46 (...) 54 60 61

Re: [R] problème d'import de fichier

2005-10-21 Thread Yves Magliulo
hello!! si tu veux plus de réponse, you should speak in english here. mais ok je vais essayer de te répondre vu ton erreur, tu devrais regarder la ligne 15 qui doit ne pas avoir le même nombre de séparateur que tes autres lignes. de façon générale, quand tu utilises read.table, assure toi que,

Re: [R] x axis

2005-10-17 Thread Yves Magliulo
hi, first shut down automatic axis in plot with plot( ...,axes=F) then customize your axis with ?axis axis(side, at = NULL, labels = TRUE, tick = TRUE, line = NA, pos = NA, outer = FALSE, font = NA, vfont = NULL, lty = solid, lwd = 1, col = NULL, padj = NA, ...) you can

Re: [R] shell scripts in R

2005-10-14 Thread Yves Magliulo
Le ven 14/10/2005 à 00:39, Benedict P. Barszcz a écrit : Dnia piątek, 14 października 2005 00:13, Sundar Dorai-Raj napisał: Don't you mean system(ls)? See ?system. Arguments: command: the system command to be invoked, as a string. This is the kind of obstacles a newbie has to

Re: [R] p-level in packages mgcv and gam

2005-09-29 Thread Yves Magliulo
of your smooth, how your model is improved, if the parametrization look fine according of what you want/know of your predictor etc... but that's my point of view, maybe someone have Criterion more precise. Yves Le jeu 29/09/2005 à 09:55, Henric Nilsson a écrit : Yves Magliulo said the following

[R] Smooth terms significance in GAM models

2005-09-23 Thread Yves Magliulo
hi, i'm using gam() function from package mgcv with default option (edf estimated by GCV). G=gam(y ~ s(x0, k = 5) + s(x1) + s(x2, k = 3)) SG=summary(G) Formula: y ~ +s(x0, k = 5) + s(x1) + s(x2, k = 3) Parametric coefficients: Estimate std. err.t ratioPr(|t|) (Intercept)

[R] Estimate predictor contribution in GAM models

2005-09-20 Thread Yves Magliulo
hi, i'm using gam() function from package mgcv. if G is my gam object, then SG=summary(G) Formula: y ~ +s(x0, k = 5) + s(x1) + s(x2, k = 3) Parametric coefficients: Estimate std. err.t ratioPr(|t|) (Intercept) 3.462e+07 1.965e+05 176.2 2.22e-16 Approximate

Re: [R] How to change x axes' range

2004-12-08 Thread Yves Magliulo
matplot(.., axes=FALSE) then custozie by yourself ayour axes -axis(1,at=1:4,label=c(11 Vars,10 Vars,6Vars,4Vars)) ?axis and see label option for more details HTH, Yves MAGLIULO, PARIS Le mer 08/12/2004 à 11:36, Qin Liu a écrit : Hi, there: When plot ann predicted results I need

[R] install bug with specific JPEG library by exporting CPPFLAGS variable

2004-12-08 Thread Yves Magliulo
Hi there, I think I have found a small problem in the /my/path/R-2.0.1/src/modules/X11/MakeFile generation. During the configure step, I have specified a specific JPEG library by exporting CPPFLAGS variable. All compilation works well for individual files in the src/modules/X11/ directory, but

Re: [R] Gam() function in R

2004-12-06 Thread Yves Magliulo
. thanks! yves magliulo, Paris. Le lun 06/12/2004 à 09:09, Jari Oksanen a écrit : On 6 Dec 2004, at 7:36, Janice Tse wrote: Thanks for the email. I will check that out However when I was doing this :gam(y~s(x1)+s(x2,3), family=gaussian, data=mydata )it gives me the error

Re: [R] Gam() function in R

2004-12-06 Thread Yves Magliulo
so mgcv package is the one i need! indeed, i want integrated smoothness selection and smooth interactions rather than stepwise selection. i have a lot of predictor, and i use gam to select those who are efficient and exclude others. (using p-value) thanks a lot for those precious information.

Re: [R] help! a urgent question

2004-11-19 Thread Yves Magliulo
hi, you can download The R Reference Index contains all help files of the R standard and recommended packages in printable form at http://cran.r-project.org/doc/manuals/fullrefman.pdf you will find other Manuals at http://www.r-project.org/ have a nice week-end all. -- Yves Magliulo

Re: [R] Combining columns of different length

2004-10-26 Thread Yves Magliulo
hi, something like : toto is your data.frame toto= V1 V2 3 2 4 2 5 8 and tata= 3 3 so to do what you want make cbind(toto,c(tata,rep(NA,len=(length(toto$V1)-length(tata) maybe there is a easier way but it's work! -- -- Yves Magliulo [EMAIL PROTECTED] RD Engineer, CLIMPACT

Re: [R] A question in R

2004-10-12 Thread Yves Magliulo
hi, check R-data import/export document at this link cran.r-project.org/doc/manuals/R-data.pdf let me give you an advice for further questions : it's easiest to find answer to your question searching directly with google -- -- Yves Magliulo [EMAIL PROTECTED] RD Engineer, CLIMPACT Tel

Re: [R] Warning: number of items to replace is not a multiple of replacement length

2004-09-29 Thread Yves Magliulo
/posting-guide.html -- -- Yves Magliulo [EMAIL PROTECTED] RD Engineer, CLIMPACT Tel. : +33 (0) 1 44 27 34 31 Fax. : +33 (0) 1 44 27 49 96 Universite Pierre et Marie Curie Boite 101 - Tour 45 - 5eme etage - Couloir 45/46 4 place Jussieu, 75252 Paris CEDEX 05, France

Re: [R] Problems with png()

2004-09-06 Thread Yves Magliulo
hi, here some alternative solution. 1-use postscript function then convert it in png in shell with convert command (higher resolution) 2-use dev2bitmap(file, type = png256) many other type available ?dev2bitmap for more info -- -- Yves Magliulo [EMAIL PROTECTED] RD Engineer, CLIMPACT

Re: [R] obtaining exact p-values in mixed effects model

2004-09-01 Thread Yves Magliulo
://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- -- Yves Magliulo [EMAIL PROTECTED] RD Engineer, CLIMPACT Tel. : +33 (0) 1 44 27 34 31 Fax. : +33 (0) 1 44 27 49 96 Universite Pierre et Marie Curie Boite 101 - Tour 45

Re: [R] Sort a data frame

2004-07-20 Thread Yves Magliulo
: [EMAIL PROTECTED] Web:www.frs.fo __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- -- Yves Magliulo [EMAIL

Re: [R] extract columns from a dataframe

2004-07-05 Thread Yves Magliulo
__ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- -- Yves Magliulo [EMAIL PROTECTED] Climatology research department, CLIMPACT Tel. : +33 (0) 1 55 07 85 77 Fax. : +33

[R] gam

2004-06-16 Thread Yves Magliulo
)). Does someone know how? the purpose is to access to equation of smooths terms in order to have the equation of my additive model. best regards, -- Yves Magliulo, Climatology research departement [EMAIL PROTECTED] Climpact __ [EMAIL PROTECTED] mailing list