Re: [R] Issue installing gamm4 for R v4.2.2

2022-11-22 Thread Nedelec, Pierre via R-help
Hi All, Thank you for your help. It was an issue with Mac packages indeed. I found that thread and it resolved it in my case also. https://stackoverflow.com/questions/72409270/getting-error-while-installing-lme4-package-in-r-in-macosx-monterey Best, Pierre From

[R] Issue installing gamm4 for R v4.2.2

2022-11-15 Thread Nedelec, Pierre via R-help
vate/var/folders/2q/qbb1m4q52md022yt2cqbd9mmgn/T/RtmpYkdEv8/downloaded_packages� Warning messages: 1: In install.packages("gamm4") : installation of package �minqa� had non-zero exit status 2: In install.packages("gamm4") : installation of package �lme4� had non-zer

[R] problems with greek language on windows

2020-01-28 Thread Pierre Ratinaud
ble to use greek language in R in this context ? Best regards Pierre Ratinaud -- Pierre Ratinaud Maître de conférences Département des Sciences de l'Education et de la Formation Laboratoire LERASS : http://www.lerass.com/ Université Toulouse - Jean Jaurès : http://www.univ-tlse2.fr/ tel :

[R] Fractional deviance for lambda in cv.glmnet for non-negative Lasso

2019-11-19 Thread Pierre Maho
Hi, I want to solve the following optimisation problem: [image: \hat{\beta} = \arg \min_{\beta \geq 0} \| y-A\beta \|_2^2 + \lambda \|\beta\|_1] For that, I am using glmnet package (cv.glmnet for finding 𝜆 and lower.limits = 0 to impose non-negativity). I would like to modify the fdev parameter

[R] Glmnet R - can't modify fdev parameter when lower = 0

2019-11-17 Thread Pierre Maho
Hi, I want to solve the following optimisation problem: [image: \hat{\beta} = \arg \min_{\beta \geq 0} \| y-A\beta \|_2^2 + \lambda \|\beta\|_1] For that, I am using glmnet package (cv.glmnet for finding 𝜆 and lower.limits = 0 to impose non-negativity). I would like to modify the fdev parameter

[R] How to extract values after using metabin from the package meta?

2017-05-04 Thread jan Pierre
Hello, I’m trying to do a meta-analysis with R. I tried to use the function metabin from the package meta : data <- data.frame(matrix(rnorm(40,25), nrow=17, ncol=8)) centres<-c("SVP","NANTES","STRASBOURG","GRENOBLE","ANGERS","TOULON","MARSEILLE","COLMAR","BORDEAUX","RENNES","VALENCE","CAEN","NAN

[R] Announcing the biggest French-speaking R Conference in North America - R à Québec 2017

2017-01-10 Thread Pierre Racine
We are posting this info here because most French-speaking R users in North America use this forum even if English is not their first language. Le comité scientifique du Colloque R à Québec 2017 (http://raquebec.ulaval.ca) recueille jusqu'au 15 janvier 2017 des propositions de conférences toucha

Re: [R] Rmpfr drives me *Rmpfr*

2016-08-16 Thread Jean-Pierre Mueller
Hello, Before running R, in a terminal, type:  sudo apt-get install libmpfr-dev -y HTH, Le mercredi 17 août 2016 à 07:46 +0200, Ferri Leberl a écrit : > Thank you for your answer. > The installation of Rmpfr ends with an error: > > checking for mpfr.h... no > configure: error: Header file mpf

[R] Package rms: c-statistic from lrm function with weights

2016-06-15 Thread Marie-Pierre Sylvestre
Dear list, I am using the lrm function from the rms package to estimate a logistic model with weights. The c-statistic (or area under the curve) is part of the lrm output. To understand how the weights enter the computation of the c-statistics, I looked at the script of lrm and lrm.fit but I am o

[R] barplot -issues with axis and labels not appearing

2015-07-21 Thread Pierre Micallef
or y axis appear on my graph (bars appear as if they are floating)(2) no axis labels are showing up (3) chart title is too high on graph and is being cut off from view Please can you help solve these issues? Thanks Pierre

[R] Frequencies for a list of vectors

2014-08-05 Thread Marie-Pierre Sylvestre
Dear R users, I have a list of vectors (list is called HTNlist). Each vector is of length 1 to 4 and takes only 0 and 1 as values. E.g. head(HTNlist) $`30008` [1] 1 0 1 0 $`60008` [1] 0 0 1 0 $`90008` [1] 0 0 1 0 $`17` [1] 1 $`130001` [1] 0 1 $`130007` [1] 1 0 1 0 I would like to obtain

Re: [R] R and external C library " cannot open shared object file" while LD_LIBRARY_PATH is set

2014-07-28 Thread Pierre Lindenbaum
ot able to load the library. Pierre (...) ( cross-posted on SO: http://stackoverflow.com/questions/24955829/ ) I'm building a C extension for R, this library also uses the HDF5 library. I compiled a dynamic library (gcc flags: -fPIC -shared -Wl,-soname,libmy.so the library seems to be

[R] R and external C library " cannot open shared object file" while LD_LIBRARY_PATH is set

2014-07-26 Thread Pierre Lindenbaum
h/to/my/src/libmy.so': /path/to/my/src/libmy.so: undefined symbol: H5T_C_S1_g Execution halted I also tried load("my.so") instead of "libmy.so". How can I load my dynamic library ? Thanks Pierre __ R-help@r-project.org

Re: [R] update.formula() to add interaction terms

2014-02-04 Thread Marie-Pierre Sylvestre
in the formula is a compound like x4:x5 and you > want x1:x4:x5 added you will need to do more work (look at the > 'factors' attribute of terms()'s output) - currently it adds x1:x5. > > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > > > -Origina

[R] update.formula() to add interaction terms

2014-02-03 Thread Marie-Pierre Sylvestre
Hi, I have a list of formulae that I need to modify. For each formula, I need to add an interaction term between the first independent variable and the last one. I want to write a function to achieve that because the list of formulae to modify is long. For example, if the formula is y ~ x1 + x2 +

[R] object not found in a cca analysis with vegan package

2013-09-24 Thread Pierre Paradis
I'm a beginner, trying to do a canonical correspondance analysis with the vegan package, but I come up with that : > library(vegan) > abondance.txt <- read.table("C:/Users/abondance.txt",header=T,sep="\t") > abondance <- cca(abondance[,-c(1,2)]~altitude, data=abondance, nhea) Error in model.fra

[R] [R-pkgs] Probabilistic neural network (PNN)

2013-05-18 Thread Pierre-Olivier Chasset
. It is delivered with four functions - learn, smooth, perf and guess - and a dataset. The functions are documented with examples and provided with unit tests. Continue reading at http://flow.chasset.net/r-pnn/ Have a good week-end, Pierre-Olivier Chasset

Re: [R] hierarchical clustering with pearson's coefficient

2013-03-29 Thread Pierre Antoine DuBoDeNa
Anyone for that question? 2013/3/28 Pierre Antoine DuBoDeNa > Hello, > > I want to use pearson's correlation as distance between observations and > then use any centroid based linkage distance (ex. Ward's distance) > > When linkage distances are formed as

[R] hierarchical clustering with pearson's coefficient

2013-03-28 Thread Pierre Antoine DuBoDeNa
re it? so that the wikipedia's statement "To apply a recursive algorithm under this objective function, the initial distance between individual objects must be (proportional to) squared Euclidean distance." is valid? Best, Pierre [[alternative HTML version deleted]] __

[R] Constrain slope in segmented package

2013-03-12 Thread Pierre Hainaut
the way, I don't understand the ".-.-x" code in the update function. Thanks a lot by anticipation ! Pierre -- Pierre Hainaut Université Catholique de Louvain Earth and Life Institute Croix du Sud, 2 - bât. De Serres, B337 - bte L7.05.14 1348 Louvain-la-Neuve - Belgique +3210479

[R] pgirmess install under linux mint 14 / a solution

2013-02-22 Thread Pierre-Henri Puech
n of libgdal1 dev lib & dependances # to sort out problems of rgdal * in R installation of rgdal & dependances # to prevent dependance problems of pgirmess *then* installation of pgirmess & dependances hope this may help ubuntu / gnu linux users best p --

[R] pearson's correlation and cross-correlation issue

2013-02-14 Thread Pierre Antoine DuBoDeNa
s any code for that. Best, Pierre [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

[R] evaluation clusters

2013-02-04 Thread Pierre Antoine DuBoDeNa
Hello, I was wondering if there is any package with several metrics for cluster estimation (estimating the k=2...10) (or evaluation..). Also if instead of euclidean distance in k-means I use some other kind of distance.. then for those evaluation metrics which takes into consideration the distanc

[R] R function for computing Simultaneous confidence intervals for multinomial proportions

2012-12-18 Thread Pierre THIRIET
e intervals and sample size determination for multinomial proportions. Journal of the American Statistical Association, 90:366-369. Paper available at http://tx.liberal.ntu.edu.tw/~purplewoo/Literature/!

[R] Fitting a multinomial model to a multi-way factorial design with repeated measures: help on package and syntax

2012-12-06 Thread Pierre THIRIET
he boxplot (run R code) as graphical representation and would be able to assign letters for pairwise comparisons of positions within a treatment, and also of treatments within a position level. I hope I have been clear and I thank you in advance, Best regards, Pierre T

Re: [R] twitteR location?

2012-08-15 Thread Jean-Pierre Müller
Have look at : http://stackoverflow.com/questions/11674842/how-to-extract-tweet-geocode-in-twitter-package-in-r/11678037#11678037 HTH, -- Jean-Pierre Müller SSP / Anthropole 4123 / UNIL / CH - 1015 Lausanne Voice:+41 21 692 3116 / Fax:+41 21 692 3115 Please avoid sending me Word or

[R] time-series classification with k-nn

2012-06-17 Thread Pierre Antoine DuBoDeNa
Hello, i am trying to apply k-nn classification for my time-series, however the euclidean distance is not the best choice as the features i use are not all normalized (others have values form 0-1 others are negative etc.) and also it doesn't do any feature evaluation and give different weights to

Re: [R] time-series statistics collection

2012-06-01 Thread Pierre Antoine DuBoDeNa
> Brockwell/Davis (2002): Introduction to Time Series and Forecasting. > Cowpertwait/Metcalfe (2009): Introductory Time Series with R. > Cryer/Chan (2008): Time Series Analysis: With Applications in R. > > for some general introductions of using time series in r. > > > > > >

[R] time-series statistics collection

2012-06-01 Thread Pierre Antoine DuBoDeNa
> > Hello, > > I am trying to collect several global measures or statistics for > time-series as well as packages of R that can compute them. I have found > several of them in papers and books, but the literature is so big i am sure > i am missing several of them. > > skewness > kurtosis > min > ma

[R] time-series statistics collection

2012-05-31 Thread Pierre Antoine DuBoDeNa
Hello, I am trying to collect several global measures or statistics for time-series as well as packages of R that can compute them. I have found several of them in papers and books, but the literature is so big i am sure i am missing several of them. skewness kurtosis min max mean SD trend season

[R] OpenStreetMap Library

2012-05-14 Thread Pierre-Olivier Chasset
t or on the command line Have you got any idea in order to get it works? Best regards, Pierre-Olivier [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

[R] r-square for non-linear regression

2012-04-25 Thread Pierre Grison
Hi, I saw you discussed about the meaning of the R squared in a nls regression. Do you have a source or a quotation please? kind regards, Pierre Grison Tel: 06 01 79 74 22 Mail: pgri...@hotmail.fr [[alternative HTML version deleted

[R] meta-analysis, outcome = OR associated with a continuous independent variable

2012-04-04 Thread Marie-Pierre Sylvestre
Hello everyone, I want to do a meta-analysis of case-control studies on which an OR was computed based on a continuous exposure. I have found several several packages (metafor, rmeta, meta) but unless I misunderstood their main functions, it seems to me that they focus on two-group comparisons (bi

[R] Strange R behavior for product of two sum of integers

2011-10-19 Thread Lafaye de Micheaux Pierre
Dear gentlemen, Can you explain me why the following happens (any OS I think, and even on 64 bits)? > sum(1000:1205)^2 [1] 51581223225 > sum(1000:1205)*sum(1000:1205) [1] NA Warning message: In sum(1000:1205) * sum(1000:1205) : NAs produced by integer overflow Best, Pierre -- Pierre

[R] optimization with Sparse matrices

2011-06-16 Thread Dube, Jean-Pierre
ing sparse matrices for some of the data inputs? Thanks, JP ____ Jean-Pierre H. Dubé Sigmund E. Edelstone Professor of Marketing The University of Chicago | Booth School of Business 5807 S. Woodlawn Avenue Chicago, IL 60637 Tel: (773)-834-5377 e-ma

Re: [R] Get contour of a map

2011-05-24 Thread Pierre Bruyer
} dev.off() } Here I use the solution of Duncan Murdoch, it is almost that I want, now I try to resolve the default of his methods, e.g. I don't want smooth the contour ( the border) of my raster. Thank for your help :-)) (Sorry for the spam, I have f

Re: [R] Get contour of a map

2011-05-24 Thread Pierre Bruyer
constrains, and I don't see how gPolygon can help me (in it doc, I have seen this function cut polygon in a set of subpolygon. Le 24 mai 2011 à 10:18, Barry Rowlingson a écrit : > On Mon, May 23, 2011 at 8:00 PM, Duncan Murdoch > wrote: >> On 11-05-23 5:55 AM, Pierre Bruyer w

Re: [R] Get contour of a map

2011-05-24 Thread Pierre Bruyer
24 mai 2011 à 10:52, David Winsemius a écrit : > > On May 23, 2011, at 5:55 AM, Pierre Bruyer wrote: > >> Hello everybody, >> >> I search a function which returns the contour of map with levels like >> contourLines, but I would like this function return the border

Re: [R] Get contour of a map

2011-05-24 Thread Pierre Bruyer
I think it is a very good idea, thank a lot, I will try this! Le 23 mai 2011 à 21:00, Duncan Murdoch a écrit : > On 11-05-23 5:55 AM, Pierre Bruyer wrote: >> Hello everybody, >> >> I search a function which returns the contour of map with levels like >> contourL

[R] Get contour of a map

2011-05-23 Thread Pierre Bruyer
. Thanks in advance Pierre Bruyer __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] recursive function

2011-05-19 Thread Tremblay, Pierre-Olivier
nyone could point me in the right direction as on how to code conditional events in R that would be much appreciated. many thanks in advance Pierre-Olivier [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://st

Re: [R] Smooth contour of a map

2011-05-18 Thread Pierre Bruyer
. Le 18 mai 2011 à 15:18, David Winsemius a écrit : > You may be looking for the par settings of xaxs="i", yaxs="i", which if you > add them to the plot call will prevent the "regular" behavior of adding 4% > padding to the axis widths. > > ?par >

Re: [R] Smooth contour of a map

2011-05-18 Thread Pierre Bruyer
found how delete it. Someone have a solution please? Pierre Bruyer ##smooth contour contours <- contourLines(V2b,levels=paliers) par(mar=c(0,0,0,0)) plot(1,col="white",main="polygon()", asp = 1, axes = FALSE, ann = FALSE,xlim=c

Re: [R] Smooth contour of a map

2011-05-17 Thread Pierre Bruyer
27;t get great results. A better way is to use a finer > grid, e.g. by fitting a smooth surface to your set of points and using > predictions from the model to interpolate. > > Duncan Murdoch > > > On 17/05/2011 9:35 AM, Pierre Bruyer wrote: >> I work with large dataset

Re: [R] Smooth contour of a map

2011-05-17 Thread Pierre Bruyer
evel[5,4] <- 102 level[6,4] <- 33 level[7,4] <- 59 level[8,4] <- 63 level[9,4] <- 14 level[10,4] <- 51 Le 17 mai 2011 à 15:17, Duncan Murdoch a écrit : > On 17/05/2011 8:24 AM, Pierre Bruyer wrote: >> Thank you for your answer, but the function spline() (and a lot of ot

Re: [R] Smooth contour of a map

2011-05-17 Thread Pierre Bruyer
;=length(level[,1])) { paliers <- c(paliers,level[i,1]) i <- i+1 } paliers <- c(paliers, 1.0E300) Le 17 mai 2011 à 13:05, Duncan Murdoch a écrit : > On 11-05-17 5:58 AM, Pierre Bruyer wrote: >> I'm a French develope

[R] Smooth contour of a map

2011-05-17 Thread Pierre Bruyer
et of colors, which is define by levels. The problem is the curve of my contour is so linear, and I would like a more continuous contour. I use the function fitted.contour to draw my map. Thanks, Pierre Bruyer __ R-help@r-project.org mailing list

Re: [R] Extracting the dimnames of an array with variable dimensions

2011-05-16 Thread Pierre Roudier
Thanks Bill, Your solution seems moreover more efficient, in terms of conputing cost. Cheers, Pierre 2011/5/16 : > Here is an alternative solution > >> foo <- array(data = rnorm(32), dim = c(4,4,2), > + dimnames=list(letters[1:4], LETTERS[1:4], letters[5:6])) >>

Re: [R] Extracting the dimnames of an array with variable dimensions

2011-05-15 Thread Pierre Roudier
Hi Dennis, Thanks for your answer, it works very well - clever way to sort the problem! Cheers, Pierre 2011/5/16 Dennis Murphy : > Hi: > > Does it have to be an array?  If all you're interested in is the > dimnames, how about this? > > library(plyr) > foo <- ar

[R] Extracting the dimnames of an array with variable dimensions

2011-05-15 Thread Pierre Roudier
[3], drop=FALSE] # But it does require to know the dimensions of foo I would like to do something like: ind <- which(foo > 0, arr.ind = TRUE) t(apply(ind, 1, function(x) unlist(dimnames(foo[x, drop=FALSE] but in that case the dimnames are dropped. Any suggestion? Cheers, Pierre

Re: [R] another question on shapefiles and geom_point in ggplot2

2011-04-03 Thread Pierre Roudier
depend on the package gpclib, which has a restricted licence. It is disabled by default; to enable gpclib, type gpclibPermit() Checking rgeos availability as gpclib substitute: TRUE > ?gpclibPermit Pierre [1] http://cran.r-project.org/web/packages/rgeos/ [2] http://trac.osgeo.org/geo

Re: [R] Problem with Snowball & RWeka

2011-03-30 Thread Jean-Pierre Müller
elate it back to how you invoked it. > This likely points to a missing object > in your call or maybe the installation lacked a dependency as this > occured during init, but hard to speculate with what you have provided. > You could try reinstalling and check for errors. > > >

[R] variance explained by each term in a GAM

2011-03-10 Thread Pierre Kleiber
d as in Simon's example, i.e. (deviance(b1)-deviance(b))/deviance(b0) ## prop explained by fv1 (deviance(b2)-deviance(b))/deviance(b0) ## prop explained by x1 etc. Am I more or less on the right track? Thanks, Pierre On Tue Oct 9 20:25:24 CEST 2007, Simon Wood wrote: > I think that

Re: [R] bug in qplot (library ggplot2)

2011-01-14 Thread Pierre-Olivier Chasset
The corrected code: library(ggplot2) for (i in 1:2) { png(file=paste('test ',i,'.png',sep='')) print(qplot(carat, data=diamonds, fill=color,geom='histogram')+scale_y_continuous(i)) dev.off() } Thanks Le 14/01/2011 02:17, Pierre-Olivier Chasset a é

Re: [R] bug in qplot (library ggplot2)

2011-01-14 Thread Pierre-Olivier Chasset
David, I haven't found any FAQs for ggplot. Do you mean the http://cran.r-project.org/faqs.html? In this one, no subject seems to be linked with my problem about qplot in a loop. Best regards, Pierre-Olivier Le 14/01/2011 05:43, David Winsemius a écrit : On Jan 13, 2011, at 8:17 PM, P

[R] bug in qplot (library ggplot2)

2011-01-13 Thread Pierre-Olivier Chasset
wing code doesn't make any file: /library(ggplot2) for (i in 1:2) { png(file=paste('test ',i,'.png',sep='')) qplot(carat, data=diamonds, fill=color,geom='histogram')+scale_y_continuous(i) dev.off() } / Have you got any idea? Have I made someth

[R] pairs and panel.smooth for two groups

2010-09-08 Thread Marie-Pierre Sylvestre
r each subplot. I don't know how to tell it to do a smooth for each of my group in each of the subplot. Any help would be appreciated! Best, Marie-Pierre __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] Error on random forest variable importance estimates

2010-08-08 Thread Pierre Dubath
Hello Andy, Thank you for your quick and helpful reply. I will try to follow your suggestions. Also, thank you for the R implementation of random forest. It is very useful for our work. Best, Pierre Liaw, Andy wrote: From: Pierre Dubath Hello, I am using the R randomForest package to

[R] Error on random forest variable importance estimates

2010-08-06 Thread Pierre Dubath
any idea how to estimate this? From the variation between runs or with the help of "importanceSD" somehow? Thanks very much in advance for any help. Pierre Dubath __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-he

[R] Multivariate Kalman Filter with inequality constraints

2010-07-08 Thread Pierre Lapointe
I'm trying to do a multivariate Kalman Filter with inequality constraints. Specifically, I'm trying to force all betas to be positive in a y = beta * x equation. Any ideas? Pierre Lapointe [[alternative HTML version deleted]] _

[R] Converting a data set from 'long' format to 'interval' format

2010-03-24 Thread Marie-Pierre Sylvestre
Hi, I have a data set in which the variable 'dose' is time-varying. Currently, the data set is in a long format, with 1 row for each time unit of follow-up for each individual "Id". It looks like this: orig.data <- cbind(Id = c(rep(1,4), rep(2,5)), time = c(1:4, 1:5), dose = c(1,1,1,0,1,0,1,1,0)

[R] suppress output from step function

2009-04-23 Thread Pierre Moffard
st step produces all the output that would have been produced even if I hadn't named the result SS. SS<-step(linear_model,print=FALSE) - this doesnt work either. Does anyone have any ideas? Thank you all in advance. Best, Pierre [[alternat

[R] Re : Running random forest using different training and testing schemes

2009-04-12 Thread Pierre Moffard
you need to include in your code something like: tree<-rpart(result~., data, control=rpart.control(xval=10)). this xval=10 is 10-fold CV. Best, Pierre De : Chrysanthi A. À : r-help@r-project.org Envoyé le : Dimanche, 12 Avril 2009, 17h26mn 59s Objet :

[R] Re : Running random forest using different training and testing schemes

2009-04-12 Thread Pierre Moffard
ining and the test set you can do (n the number of data points): n<-length(data[,1]) indices<-sample(rep(c(TRUE,FALSE),each=n/2),round(n/2),replace=TRUE) training_indices<-(1:n)[indices] test_indices<-(1:n)[!indices] Then, data[train,] is the training set and data[test,] is the test s

[R] subscript out of bounds in eqscplot problem

2009-04-08 Thread Pierre Moffard
highly appreciated. Thanks in advance, Pierre [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

[R] gamm (mgvc) and time-varying coefficient model

2009-03-01 Thread Marie-Pierre Sylvestre
Dear R users, I have repeated measurements on individuals. I want to estimate the time-varying effect of a factor variable X (taking three levels), e.g. a model in the spirit of Hastie and Tibshirani (1993). I am considering using the package "mgvc" which implements generalized additive model

[R] gamm (mgvc) and time-varying coefficient model

2009-03-01 Thread Marie-Pierre Sylvestre
Dear R users, I have repeated measurements on individuals. I want to estimate the time-varying effect of a factor variable X (taking three levels), e.g. a model in the spirit of Hastie and Tibshirani (1993). I am considering using the package "mgvc" which implements generalized additive models,

[R] regexp capturing group in R

2009-02-24 Thread pierre
tried: txt <- "blah blah start=20080101 end=20090224" m <- regexpr("start=(\\d{8}).*end=(\\d{8})", filename, perl=T); dates = substring(filename, m, m+attr(m,"match.length")-1); but I get the whole matching substring... Any idea? ~Pierre _

Re: [R] Event sequence analysis

2009-02-02 Thread Jean-Pierre Müller
tat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Jean-Pierre Müller SSP / Anthropole 4123 / UNIL / CH - 1015 Lausanne Voice:+41 21

Re: [R] Calculate Specificity and Sensitivity for a given threshold value

2008-11-13 Thread Pierre-Jean-EXT.Breton
dea for my problem ? Pierre-Jean -Original Message- From: Frank E Harrell Jr [mailto:[EMAIL PROTECTED] Sent: Thursday, November 13, 2008 5:00 PM To: Breton, Pierre-Jean-EXT R&D/FR Cc: r-help@r-project.org Subject: Re: [R] Calculate Specificity and Sensitivity for a given threshold

[R] creating a package with internal (hiden) S4 classes and methods

2008-09-11 Thread Marie-Pierre Sylvestre
Dear R users, I am writing my first R package and it finally past the R CMD check. Before sending it to CRAN, I would like to check a few issues. 1. I created s4 classes and methods in the package but I want only 1 function (foo) to be available to users and documented in the manual that goe

[R] building a package that contains S4 classes and methods

2008-09-09 Thread Marie Pierre Sylvestre
Hello R users, I am trying to make a my first package and I get an error that I can understand. The package is build out of three files (one for functions, 1 for s4 classes and 1 for s4 methods). Once I source them I run package.skeleton( name="TDC" ) within a R session and I get Creating di

[R] Tidying up code - Warning message: deparse may be incomplete

2008-08-28 Thread Marie Pierre Sylvestre
Dear R users, I am currently writing a R package and to do so I am following the guidelines in manual 'Writing R extensions'. In Section 3.1, it is suggested to tidy up the code using a file containing the following: options(keep.source = FALSE) source("myfuns..R") dump(ls(all = TRUE), file = new

[R] random error with lme for repeated measures anova

2008-08-27 Thread Jean-Pierre Bresciani
Hi, what is the appropriate syntax to get the random error correct when performing repeated measures anova with 'lme'. let's say i have 3 independent variables, with 'aov', i would write something like: aov(dep_var~(indep_var1*indep_var2*indep_var3) + Error(subject/(indep_var1*indep_var2*indep_v

Re: [R] post hoc tests two way repeated measures anova

2008-08-18 Thread Jean-Pierre Bresciani
27;, but those are obviously wrong with my design. i'm quite clueless (and i haven't found any convincing piece of information about how to correctly use 'lme' or 'lmer'). So, any advice along that line is more than welcome. JP Frank E Harrell Jr wrote: > >

[R] post hoc tests two way repeated measures anova

2008-08-15 Thread Jean-Pierre Bresciani
Hi, is there a specific/appropriate function/package to perform post hoc tests when running a two way repeated measures anova? I'm looking for something that will be equivalent to the 'TukeyHSD()' for between subjects anova (with 'aov()'). For one way repeated measures anova, the 'pairwise.t.test

Re: [R] Compiling only some C files

2008-08-11 Thread Lafaye de Micheaux Pierre
other OSes that lack standard compiler tools. If you have some ideas to improve the way i planned to do it, please let me know. I could send you what i have already done. Best, Pierre Duncan Murdoch a écrit : Barry Rowlingson wrote: 2008/8/11 Lafaye de Micheaux Pierre <[EMAIL PROTEC

Re: [R] Compiling only some C files

2008-08-11 Thread Lafaye de Micheaux Pierre
Thank you very much! This is an easy and very effective solution. Thanks again, Pierre Barry Rowlingson a écrit : > 2008/8/11 Lafaye de Micheaux Pierre > <[EMAIL PROTECTED]>: > >> Thank you Barry, >> >> In fact, these are only C files (not C++) and i use the

Re: [R] Compiling only some C files

2008-08-11 Thread Lafaye de Micheaux Pierre
other OSes that lack standard compiler tools. If you have some ideas to improve the way i planned to do it, please let me know. I could send you what i have already done. Best, Pierre Duncan Murdoch a écrit : Barry Rowlingson wrote: 2008/8/11 Lafaye de Micheaux Pierre <[EMAIL PROTEC

Re: [R] Compiling only some C files

2008-08-11 Thread Lafaye de Micheaux Pierre
les in one unique file because many users have to provide their own C file. So i would really like to use some kind of configure or Make script, that should be put in the source directory tree of the R package. Pierre Barry Rowlingson a e'crit: > 2008/8/11 Lafaye de Micheaux Pierre &

[R] Compiling only some C files

2008-08-11 Thread Lafaye de Micheaux Pierre
success... I hope to be clear enough. Best regards, Pierre Lafaye de Micheaux -- Pierre Lafaye de Micheaux Bureau 210, bâtiment BSHM 1251 avenue centrale BP 47 38040 GRENOBLE Cedex 09 FRANCE Tél.: 04.76.82.58.73 / Fax: 04.76.82.56.65 [EMAIL PROTECTED]

[R] Fixing the coefficient of a regressor in formula

2008-05-28 Thread Marie-Pierre Sylvestre
Dear R users, I want to estimate a Cox PH model with time-dependent covariates so I am using a counting process format with the following formula: Surv(data$start, data$stop, data$event.time) ~ cluster(data$id) + G1 + G2 + G3 + G4 + G5 +G6 Gs represent a B-spline basis functions so they s

[R] Un-reproductibility of SVM classification with 'e1071' libSVM package

2008-05-13 Thread Pierre Dangauthier
Hello, When calling several times the svm() function, I get different results. Do I miss something, or is there some random generation in the C library? In this second hypothesis, is it possible to fix an eventual seed? Thank you Pierre ### Example library('e1071') x = rnorm(100) #

[R] assigning and saving datasets in a loop, with names changing with "i"

2007-12-18 Thread Marie Pierre Sylvestre
Dear R users, I am analysing a very large data set and I need to perform several data manipulations. The dataset is so big that the only way I can play with it without having memory problems (E.g. "cannot allocate vectors of size...") is to write a batch script to: 1. cut the data into pieces 2.

[R] Help about nnet library

2007-11-14 Thread Pierre-Henri BUSSIERE
*h2+wo7*h3+wo8*h4+wo9*h5. So I would like to know, where is our error. Is it the sigmoid function that we used ? The transfert function ? Please answer me, I really need to know it. Thanks in advance. Pierre-Henri BUSSIERE -- Pierre-Henri BUSSIERE Ingénieur Informatique -

[R] Identifying and characterizing strings of NA in a vector

2007-10-11 Thread Marie-Pierre Sylvestre
f NA' of length 3. I can easily count the number of NA's per vector, but I am having a hard time counting the number and length of strings of NA's per vector without relying heavily on loops. I will have to perform this task for many vectors. Can somebody help? man

[R] bug: wireframe and tcltk

2007-10-09 Thread Lafaye de Micheaux Pierre
ion about this point? Best regards, Pierre Lafaye de Micheaux -- Pierre Lafaye de Micheaux Bureau 210, bâtiment BSHM 1251 avenue centrale BP 47 38040 GRENOBLE Cedex 09 FRANCE Tél.: 04.76.82.58.73 / Fax: 04.76.82.56.65 [EMAIL PROTECTED] http://www.biostatisticien.eu is.it.a.bug <- fun