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
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
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 :
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
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
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
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
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
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
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
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
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
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
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
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 +
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
. 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
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
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]]
__
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
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
--
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
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
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/!
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
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
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
> 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.
>
>
>
>
>
>
>
> 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
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
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
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
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
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
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
}
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
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
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
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
.
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.
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
.
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
>
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
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
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
;=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
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
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]))
>>
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
[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
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
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.
>
>
>
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
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 é
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
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 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
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
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
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]]
_
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)
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
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 :
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
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
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
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,
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
_
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
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
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
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
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
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
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:
>
>
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
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
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
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
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
&
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]
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
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) #
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.
*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
-
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
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
87 matches
Mail list logo