Dear R-experts,
Here below is my code,
I would like to add a vertical line on my plot, showing the median and I would
like to place some values on this graph as well, i.e. 4.3 and -8.4. How can I
do ?
Many thanks for your reply.
A=c(1,2.3,4,3.5,4.3,2.5,6.3,-0.1,-1.5,3.7,-2.3,-3.5,5.4,3.2, -1
lot points, see, well, ?points().
Hope this helps,
Rui Barradas
Em 22-10-2017 16:33, varin sacha via R-help escreveu:
> Dear R-experts,
>
> Here below is my code,
> I would like to add a vertical line on my plot, showing the median and I
> would like to place some values on this graph
Dear R-Experts,
Here below is my R code (reproducible example) to calculate the confidence
intervals around the spearman coefficient.
##
C=c(2,4,5,6,3,4,5,7,8,7,6,5,6,7,7,8,5,4,3,2)
D=c(3,5,4,6,7,2,3,1,2,4,5,4,6,4,5,4,3,2,8,9)
cor(C,D,method= "spearman")
library(boot)
myCor=function(data
ults,type="all")
Hope this helps,
Rui Barradas
On 12/10/2017 3:19 PM, varin sacha via R-help wrote:
> Dear R-Experts,
>
> Here below is my R code (reproducible example) to calculate the confidence
> intervals around the spearman coefficient.
>
> #
second returns a symmetric matrix. Just like
cor()
David L. Carlson
Department of Anthropology
Texas A&M University
-Original Message-
From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of varin sacha via
R-help
Sent: Sunday, December 10, 2017 11:07 AM
To: Rui Barradas ; R-h
Dear R-experts,
I can not install the package "matie". If somebody can tell me where the
problem is.
> install.packages("matie")
Installing package into ‘/Users/Caro/Library/R/3.3/library’
(as ‘lib’ is unspecified)
essai de l'URL
'https://cran.rstudio.com/bin/macosx/mavericks/contrib/3.3/matie_
e my brevity.
On March 28, 2018 3:38:27 PM PDT, Duncan Murdoch
wrote:
>On 28/03/2018 6:27 PM, varin sacha via R-help wrote:
>> Dear R-experts,
>>
>> I can not install the package "matie". If somebody can tell me where
>the problem is.
>>
>>>
Dear R-experts,
Here below my reproducible R code. I want to add many straight lines to a plot
using "abline"
The last fit (fast Tau-estimator, color yellow) will not appear on the plot.
What is going wrong ?
Many thanks for your reply.
##
Y=c(2,4,5,4,3,4,2,3,56,5,4,3,4,5,6,5,4,5,34,21
Many thanks Duncun,
Best,
Le samedi 31 mars 2018 à 18:05:53 UTC+2, Duncan Murdoch
a écrit :
On 31/03/2018 11:57 AM, varin sacha via R-help wrote:
> Dear R-experts,
>
> Here below my reproducible R code. I want to add many straight lines to a
> plot using "ablin
t;)
legend("topright",c("OLS", "L1", "Huber M-estimator", "Tukey", "MM", "fast
tau"),inset=0.02,lwd=2,col=c("black","blue","red","purple","green",
"yellow"),cex=.9)}
##
a as slope
abline( fast$beta, col="yellow" )
legend( "topright"
, c( "OLS"
, "L1"
, "Huber M-estimator"
, "Tukey"
, "MM"
, "fast tau"
)
Dear R-experts,
I am trying to do cross-validation for different models using the cvTools
package.
I can't get the CV for the "FastTau" and "hbrfit". I guess I have to write my
own functions at least for hbrfit. What is going wrong with FastTau ?
Here below the reproducible example. It is a si
Dear R-experts,
Doing cross-validation for 2 robust regressions (HBR and fast Tau). I can't get
the 2 errors rates (RMSE and MAPE). The problem is to predict the response on
the testing data. I get 2 error messages.
Here below the reproducible (fictional example) R code.
#install.packages("MLm
Dear R-experts,
I guess I have a problem with my fast function (fast tau estimator) here below.
Indeed, zero errors look highly suspicious. I guess there is a bug in my R
code. How could I correct my R code ?
# install.packages( "robustbase" )
# install.packages( "MASS" )
# install.packages( "q
Dear R-experts,
Here below the reproducible example. I am trying to get the average of the 100
results coming from the "lst" function. I have tried lst$mean and mean(lst). It
does not work.
Any help would be highly appreciated.
## R script for getting MedAe and MedAeSQ
Many thanks for all of you for your responses.
Best Regards,
SV
Le mardi 8 mai 2018 à 21:58:37 UTC+2, Daniel Nordlund a
écrit :
On 5/8/2018 12:26 PM, varin sacha via R-help wrote:
>
> Dear R-experts,
>
> Here below the reproducible example. I am trying to get the av
Dear R-experts,
I am trying to bootstrap (and average) the median squared error evaluation
metric for a robust regression. I can't get it. What is going wrong ?
Here is the reproducible example.
#
install.packages( "quantreg" )
library(quantreg)
crp <-c(12,14,13,24
Hi Rui and Daniel,
Many thanks for your responses. It perfectly works.
Best,
Le mardi 22 mai 2018 à 12:47:14 UTC+2, Rui Barradas
a écrit :
Hello,
Right!
I copied from the OP's question without thinking about it.
Corrected would be
bootMedianSE <- function(data, indices){
d <- d
Dear R-helpers,
I have fitted a robust regression using lmrob function from robustbase package.
I try to get the different plots for diagnostics of residuals and others. I
can't get them, a window opens but nothing appears on it (the window remains
white, no graph appears) and I get this error
Perfect, many thanks
Best
S
De : David Winsemius
À : varin sacha
Cc : R-help Mailing List
Envoyé le : Vendredi 20 mai 2016 21h13
Objet : Re: [R] Plots for lmrob function
> On May 20, 2016, at 8:28 AM, varin sacha via R-help
> wrote:
>
> Dear R-helpers,
>
&g
Dear R-experts,
I am trying to calculate the bootstrapped (BCa) regression coefficients for a
robust regression using MM-type estimator (lmrob function from robustbase
package).
My R code here below is showing a warning message ([1] "All values of t are
equal to
22.2073014256803\n Can not cal
t working
> On 03 Jul 2016, at 13:47 , varin sacha via R-help
> wrote:
>
> Dear R-experts,
>
> I am trying to calculate the bootstrapped (BCa) regression coefficients for a
> robust regression using MM-type estimator (lmrob function from robustbase
> package).
&g
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Tue, Jul 5, 2016 at 2:51 PM, varin sacha via R-help
wrote:
> Dear Professor Dalgaard,
>
> I really thank you lots for your response. I have solved my problem. Now, I
> have tried to do the s
t;
>
> Bert Gunter
>
> "The trouble with having an open mind is that people keep coming along
> and sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
>
> On Tue, Jul 5, 2016 at 2:51 PM, varin sacha via
Bootstrapped regression coefficients from lmrob function
not working
On 06/07/2016 7:56 AM, varin sacha via R-help wrote:
> Dear Professor Dalgaard,
>
> Okay, this is what I was afraid of.
> Many thanks for your response. I know that my next question is off-topic here
> (on this web
Dear R-experts,
I have fitted a MARS regression and am trying now to plot/draw the BCa
confidence bands around the 3 fitted curves (QUALITESANSREDONDANC,
competitivite and innovation).
Here is the reproducible example.
Dataset =
data.frame(PIBparHab=c(43931,675
e? If not, why
not? If so, show us the code that failed.
Or am I missing the point?
Cheers,
Bert
Bert Gunter
"The trouble with having an open mind is that people keep coming along
and sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip
Dear R-experts,
I am having trouble while doing crossvalidation with a MARS regression
including an interaction term between a factor variable (education) and 1
continuous variable (age). How could I solve my problem ?
Here below my reproducible example.
###
install.packages("ISLR")
lib
County" comic strip )
On Mon, Oct 29, 2018 at 1:46 PM varin sacha via R-help
wrote:
>
> Dear R-experts,
> I am having trouble while doing crossvalidation with a MARS regression
> including an interaction term between a factor variable (education) and 1
> continuous vari
l de facto contains both education as a numeric predictor and
as.factor(education) as well as the interaction term age:as.factor(education).
Does that make sense modelling-wise??
-pd
> On 29 Oct 2018, at 23:50 , varin sacha via R-help
> wrote:
>
> Hi Bert,
>
> Many thanks,
Dear R-experts,
I am trying to get the bootstrapped confidence intervals of Mean squared error
(MSE) for a (G)AM model. I get an error message.
Here below the reproducible R code. Many thanks for your response.
install.packages("ISLR")
library(ISLR)
install.packages("mgc
pe <- c("norm","basic", "stud", "perc")
> boot.ci(results, type = ci.type)
>
>
>
> Hope this helps,
>
> Rui Barradas
>
> Às 20:36 de 22/11/2018, varin sacha via R-help escreveu:
>> Dea
pe = ci.type)
# view results
results
plot(results)
# get 95% confidence intervals
boot.ci(results, type="all")
# # # # # # # # # # # # # # # # # # # # # # # # # #
Le vendredi 23 novembre 2018 à 00:20:41 UTC+1, varin sacha via R-help
a écrit :
Great, many thanks Rui, it perfe
Dear R-experts,
We can fit MARS regression using the packages "earth" and/or "mda" or others
packages.
However, I am wondering if it is possible to fit a CMARS (Conic multivariate
adaptive regression splines) using R ?
I have googled "conic MARS with R software", I did not get anything, so Goog
any R
> software. Maybe you can do better.
>
> -- Bert
>
> Bert Gunter
>
> "The trouble with having an open mind is that people keep coming along and
> sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip
Dear R-experts,
I have fitted a model with 2-way and 3-way interactions.
I would like, for the 3-way interaction (year,age,by=education), to obtain
3D-plots. How could I do that ?
Many thanks for your response.
Here is the reproducible example:
#
install.packages("ISLR")
library
C, data=mydata,
groups=A,
panel=panel.bwplot.superpose, ## take control of colors of the boxes
horizontal=FALSE,
between=list(x=1, y=1))
)
On Mon, Jan 7, 2019 at 7:59 PM David Winsemius wrote:
>
> On 1/7/19 3:35 PM, varin sacha via R-help wrote:
>> Dear R-experts
Dear R-experts,
The reproducible example is below.
I am trying to use the visreg package for 2 ways-interactions and for 3
ways-interactions.
For 2 ways-interactions, everything goes fine except that the legend on the top
(2.HS Grad ; 3. Some College ; 4. College Grad) is not entire/complete, b
Dear R-Experts,
Thanks to Prof. Bonnett, I have got an R script working to calculate confidence
intervals around the semipartial correlation coefficients.
Now, I would like to calculate BCa bootstrap CIs using the boot library and the
boot.ci(results, type="all") function. How could I modify my
Hi Deepak,
In econometrics there is another test very often used : the white test.
The white test is based on the comparison of the estimated variances of
residuals when the model is estimated by OLS under the assumption of
homoscedasticity and when the model is estimated by OLS under the assump
Dear R-helpers,
I would like to know if it is possible using R (earth packages or another one)
to realize robust MARS regression RMARS or RCMARS (C mean "conic").
Some authors have done it using MOSEK combined with Salford-MARS and special
MATLAB programs.
Best Regards,
Sacha
___
Dear R-experts,
Here below you will find my R sessionInfo. I have a problem with the
library(Rcmdr). The error message is here below. How can I solve that problem ?
Many thanks.
sessionInfo()
R version 3.3.2 (2016-10-31)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: macOS Sierr
Cc : R-help Mailing List
Envoyé le : Samedi 5 novembre 2016 20h42
Objet : Re: [R] Problem with downloading library(Rcmdr)
> On Nov 5, 2016, at 11:27 AM, varin sacha via R-help
> wrote:
>
> Dear R-experts,
>
> Here below you will find my R sessionInfo. I have a problem with t
u can read here
below.
What am I doing wrong ?
> R CMD INSTALL -1 rlme_0.4.tar.gz
Erreur : unexpected symbol in "R CMD"
- Mail original -
De : Jeff Newmiller
À : varin sacha ; varin sacha via R-help
; David Winsemius
Cc : R-help Mailing List
Envoyé le : Dimanche 6 novemb
Mark,
Many thanks for your response.
Best,
S
De : Mark Sharp
Cc : R-help Mailing List
Envoyé le : Lundi 30 janvier 2017 23h04
Objet : Re: [R] Save a generated .txt (or .csv) file on the desktop
You can define the "file" argument in your call to write.csv(
Dear R-experts,
Here below my R code. The warning message is not a problem to me but there is
an error message more problematic. I understand the error message but I don't
know if it is possible to correct the error and if yes, how to correct it.
Many thanks.
n <- 60
b <- runif(n, 0, 5)
a <-
bs)
Thanks.
Md
On Tue, Oct 27, 2020 at 7:21 PM Sarah Goslee wrote:
> Hi,
>
> a is of length 60.
> b is of length 60.
> z is of length 57.
>
> What do you expect to have happen when you create y_model ? What
> happens to those other 3 observations?
>
> Sarah
>
Many thanks Duncan,
It works !
Best.
Le mardi 27 octobre 2020 à 20:49:25 UTC+1, Duncan Murdoch
a écrit :
On 27/10/2020 3:06 p.m., varin sacha via R-help wrote:
> Dear R-experts,
>
> Here below my R code. The warning message is not a problem to me but there is
&
Dear R-helpers,
Here below my R code showing warnings and error messages I don't understand.
What is going wrong ?
install.packages("devtools")
library(devtools)
install_github("ProcessMiner/nlcor")
library(nlcor)
A=c(505, 530, 419, 486, 608, 468, 519, 486, 532, 289, 5
Dear R-experts,
Here below my reproducible example. No error message but I can not get a
result. I get "NaN" as a result. I don't understand what is going on. Many
thanks for your precious help, as usual.
# # # # # # # # # # # # # # # # # # # # # # # # #
x<-c(499,491,500,517,438,495,501,525,5
rs. However this would probably break some existing
code, since there are various functions that return NULL instead of numeric(0).
-Bill
On Fri, Nov 13, 2020 at 2:05 PM varin sacha via R-help
wrote:
> Dear R-experts,
>
> Here below my reproducible example. No error message but I can not g
Dear R-Experts,
Here below my reproducible R code.
How can I get the AIC of my model (robust GAM) ?
Best Regards,
y<-c(499,491,500,517,438,495,501,525,516,494,500,453,479,481,505,465,477,520,520,480,477,416,502,503,497,513,492,469,504,482,502,498,463,504,495)
x<-c(499,496,424,537,480,484,503,57
cking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Wed, Jan 6, 2021 at 12:17 PM varin sacha via R-help
wrote:
> Dear R-Experts,
>
> Here below my reproducible R code.
> How can I get the
Dear R-experts,
Here below my R code. What is happening with my green curve ? How to solve the
problem ?
Many thanks,
##
#DATA
y<-c(499,491,500,517,438,495,501,525,516,494,500,453,479,481,505,465,477,520,520,480,477,416,502,503,497,513,492,469,504,482,502,498,463,504
ry
o <- order(x)
lines(x[o], yfit[o], col="green", lwd=2)
Duncan Murdoch
On 09/01/2021 2:58 p.m., varin sacha via R-help wrote:
>
Dear R-experts,
I dont get any error message but I can't get the green curve/line on the plot.
What is going on ? How to solve my problem ?
Here is the R code :
PIB.hab<-c(12000,34000,25000,43000,12500,32400,76320,45890,76345,90565,76580,45670,23450,34560,65430,65435,56755,8
Dear R-experts,
I have tried to reach the maintainer of the rgam package. Until now, no
response.
Since I'm in a bit of a hurry, I try to reach you because as I try to install
the rgam package using the command :
install.packages("rgam")
I get this warning message :
Warning message:
package
y to respond to these kinds of questions.
Also, I note the most obvious choice for GAM-based models in R, is
mgcv, which ships with R.
But there are other packages, if mgcv doesn't meet your needs.
On Wed, Jan 13, 2021 at 8:44 AM varin sacha via R-help
wrote:
>
> Dear R-experts,
>
>
Murdoch wrote:
>
> On 12/01/2021 2:43 p.m., varin sacha via R-help wrote:
>> Dear R-experts,
>> I have tried to reach the maintainer of the rgam package. Until now, no
>> response.
>> Since I'm in a bit of a hurry, I try to reach you because as I try to
&g
Dear R-experts,
Here below my reproducible R code. I get an error message (end of code) I can't
solve.
Many thanks for your help.
##
#Data
y=c(34000,45000,19000,48900,65000,67000,78000,9,51000,32000,54000,85000,38000,76345,87654,90990,78654,67894,56789,65432,18998,789
at the top. They should
remain untouched during your program.
However in the loop you redefine y and then use the redefined y as an
argument to robustgam() the next time through
the loop. This looks like a serious error.
HTH,
Eric
On Sun, Jan 17, 2021 at 12:20 PM varin sacha via R-help
w
thing like
> ypred = pred.robustgam( fit18, data.frame(X=Testing)
>
> 2. your logic is wrong. You define the vectors x and y at the top. They
> should remain untouched during your program.
> However in the
Dear R-experts,
Here below my R code working but I would like to get 3 values not only 1. The
value I get is, according to my R code, the variance value. My goal is to get 3
values : the bias value, the variance value and the MSE value. How to solve my
problem ?
Many thanks.
#
Dear R-experts,
I am trying to get the 8 graphs like the ones in this paper :
https://statweb.stanford.edu/~tibs/reshef/comment.pdf
My R code does not show any error message neither warnings but I d'on't get
what I would like to get (I mean the 8 graphs), so I am missing something.
What's it ?
(1,2,3), col = c("black","green","blue","red"))
plot((1:30)/10, power.cor[7,], ylim = c(0,1), main = "Circle", xlab = "Noise
Level", ylab = "Power", pch = 1, col = "black", type = 'b')
points((1:30)/10,
Rui,
Angelo,
I found it :=)
Many thanks
S.
Le dimanche 9 mai 2021 à 23:49:41 UTC+2, Angelo Canty
a écrit :
Have you looked in the pdf file (power.pdf) to which you instructed R to
send the plots?
On 2021-05-09 5:27 p.m., varin sacha via R-help wrote:
> Dear Rui,
>
> I
The comments only line is your last code line.
The result is attached.
Hope this helps,
Rui Barradas
Às 19:39 de 09/05/21, varin sacha via R-help escreveu:
> Dear R-experts,
>
> I am trying to get the 8 graphs like the ones in this paper :
> https://statweb.stanford.edu/~tibs/
ng an open mind is that people keep coming along and
> sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
>
> On Sun, May 9, 2021 at 2:59 PM varin sacha via R-help
> wrote:
>
>> Rui,
>>
>> The cr
", type = 'b')
points((1:30)/10, power.cors[6,], pch = 2, col = "green", type = 'b')
points((1:30)/10, power.cork[6,], pch = 3, col = "blue", type = 'b')
points((1:30)/10, power.dcor[6,], pch = 4, col = "red", type = 'b')
po
gt; Were you expecting image files? I don't see any plot device e.g. pdf()
> in your code.
>
> Jim
>
>> On Wed, May 12, 2021 at 6:34 PM varin sacha via R-help
>> wrote:
>>
>> Dear Experts,
>>
>> My R code was perfectly working since I decide to
Dear R-experts,
Here below a toy R code example. I would like some countries (not all of them)
"Italy", "Canada", "Greece" and "Norway" to appear in red color. The others
remaining black. How can I do that without big changes in my R code ? Indeed, I
would like my R code to remain like this as
Dear R-experts,
Here below my R code. I would need your help to improve my graph/plot.
- The x-axis to be longer not to stop at 500 value
- All the name on the y-axis to appear not only a few of them and the name
(Fribourg(f), Appenzell Rhodes Intérieures,...) to appear entire, not to be cut
Ma
ot;Schaffhouse","Berne(f)",
"Thurgovie","Valais(f)","Argovie","Appenzell Rhodes Extérieures",
"Genève","Zoug","Tessin","Neuchâtel","Vaud","Uri","Nidwald",
"Appenzell Rhodes Extérieures","Genève","Zoug","Tessin",
"Neuchâtel","Vaud","Uri","Nidwald","Berne(d)","Zurich","Obwald",
"Saint-Gall","Soleure","Lucerne","Glari
Dear R-experts,
Here is a toy example. How can I get the bootstrap confidence intervals working
?
Many thanks for your help
library(DescTools)
library(boot)
A=c(488,437,500,449,364)
dat<-data.frame(A)
med<-function(d,i) {
temp<-d[i,]
HodgesLehmann(A)
}
boot.out<-
Hi,
I really thank you a lot for your response.
Le samedi 6 novembre 2021, 02:37:46 UTC+1, David Winsemius
a écrit :
On 11/5/21 1:16 PM, varin sacha via R-help wrote:
> Dear R-experts,
>
> Here is a toy example. How can I get the bootstrap confidence intervals
> workin
Hi,
Many thanks for your responses.
Le samedi 6 novembre 2021, 08:39:22 UTC+1, Rui Barradas
a écrit :
Hello,
Às 01:36 de 06/11/21, David Winsemius escreveu:
>
> On 11/5/21 1:16 PM, varin sacha via R-help wrote:
>> Dear R-experts,
>>
>> Here is a toy ex
Dear R-experts,
Here below my R code. I am trying to do 2 things :
1) I would like to repeat this R code 1 times
2) Out of the 1 repetitions I would have liked R to tell me how many times
the "true" mean value of the population called "s" in my R example here below
is included in the 1
;- 1e3, not 1e4.
set.seed(2021)
N <- 1e3
out <- replicate(N, {
boot.out <- boot(data = dat, statistic = med, R = 1)
boot.ci(boot.out, type = "bca")$bca[, 4:5]
})
mean(out[1,] < mean(s) & mean(s) < out[2,])
Hope this helps,
Rui Barradas
Às 10:47 de 11/12/21, vari
t; N <- 1e3
> out <- replicate(N, {
> boot.out <- boot(data = dat, statistic = med, R = 1)
> boot.ci(boot.out, type = "bca")$bca[, 4:5]
> })
> mean(out[1,] < mean(s) & mean(s) < out[2,])
>
>
> Hope this helps,
>
> Rui Barradas
&g
Dear R-experts,
Here below my R code. Using "bca" in the boot.ci function at the end of my R
code works perfectly. Using "perc" and "basic" perfectly works as well. But
using "norm", my R code does not work anymore.
Is it possible to solve that problem and to make my R code work ?
#
ue in ?boot.ci.
So just change
$norm[, 4:5]
to
$norm[, 2:3]
Hope this helps,
Rui Barradas
Às 22:24 de 16/12/21, varin sacha via R-help escreveu:
> Dear R-experts,
>
> Here below my R code. Using "bca" in the boot.ci function at the end of my R
> code works perfect
Dear R-experts,
Here below my R code working but really really slowly ! I need 2 hours with my
computer to finally get an answer ! Is there a way to improve my R code to
speed it up ? At least to win 1 hour ;=)
Many thanks
library(boot)
;- sample(s, size = 5)
boot.out <- boot(data = a, statistic = med, R = 1)
boot.ci(boot.out, type = "stud")$stud[, 4:5]
})
mean(out[1, ] < mean(s) & mean(s) < out[2, ])
#[1] 0.952
Hope this helps,
Rui Barradas
Às 11:45 de 19/12/21, varin sacha via R-help escreveu:
a")$bca[, 4:5]
+ })
> #coverage probability
> mean(out[1, ] < m & m < out[2, ])
[1] 0.952
I hope this helps,
John
--
John Fox, Professor Emeritus
McMaster University
Hamilton, Ontario, Canada
Web: http://socserv.mcmaster.ca/jfox/
On 2021-12-29, 2:09 PM, "R-help o
Dear R-experts,
Here below my R code for the percentile bootstrap confidence intervals with an
error message.
Is there a way to make my R code work ?
Many thanks for your help and time.
library(boot)
s=rnorm(10,0,1)
(m<-median(s))
N <- 100
n<-
John Fox, Professor Emeritus
McMaster University
Hamilton, Ontario, Canada
Web: http::/socserv.mcmaster.ca/jfox
> On Jan 8, 2022, at 12:04 PM, varin sacha via R-help
> wrote:
>
> Dear R-experts,
>
> Here below my R code for the percentile bootstr
kages that appear to provide that
information. Have you tried any of them?
On Sun, Feb 13, 2022 at 12:44 PM varin sacha via R-help
wrote:
>
> Dear R-experts,
>
> Here below my R code. I am trying to get the exact 95% confidence intervals
> around the meanlog for lognormal, around the
Dear R-experts,
I have used the EnvStats package and the elnorm function (p. 248).
I would like to calculate the exact 95% confidence intervals around the mean,
not around the meanlog.
Here below my R code, how can I get the exact 95% CIs around the mean ?
Many thanks.
library(EnvStats)
x=rlnor
Dear R-experts,
Here below my sessionInfo( ). I cannot download the pcalg package/library. Is
my R version too old ?
sessionInfo()
R version 4.1.2 (2021-11-01)
Platform: x86_64-apple-darwin17.0 (64-bit)
Running under: macOS High Sierra 10.13.6
Matrix products: default
BLAS:
/Library/Framewo
Many thanks for your response.
Le mercredi 16 mars 2022, 21:11:06 UTC+1, varin sacha via R-help
a écrit :
Dear R-experts,
I have used the EnvStats package and the elnorm function (p. 248).
I would like to calculate the exact 95% confidence intervals around the mean,
not around
Ivan,
Here is what I get :
download.file('https://stat.ethz.ch/CRAN/src/contrib/PACKAGES','PACKAGES')
essai de l'URL 'https://stat.ethz.ch/CRAN/src/contrib/PACKAGES'
Erreur dans download.file("https://stat.ethz.ch/CRAN/src/contrib/PACKAGES", :
impossible d'ouvrir l'URL 'https://stat.ethz.ch/
djust(r$p.value, method = "holm")) :
# no non-missing arguments to min; returning Inf
p.dagitty.param.correct
#[1] Inf
Hope this helps,
Rui Barradas
Às 10:58 de 31/03/2022, varin sacha via R-help escreveu:
> Dear R-experts,
>
> Here below my R code (toy example). No
Dear R-experts,
Here below my R code working but I don't know how I can get the graph
(estimated DAG).
If you could help me to get the graph, many thanks.
###
library(pcalg)
x1<-c(508,413,426,500,568,372,484,512,529,322,544,586,480,561,567,488,450,548,526,561,435,567,537,521,5
Dear R-experts,
Here below my R code working but I would like to get the names of the variables
to appear on the graph instead of 1 ; 2 ; 3 and 4.
Is it possible ? I have tried something with colnames(x) but I get an error
message.
Many thanks for your help.
#
librar
Dear R-experts,
While comparing groups, it is better to assess confidence intervals of those
differences rather than comparing confidence intervals for each group.
I am trying to calculate the CIs of the difference between the two Cramer's V
and not the CI to the estimate of each group’s Cramer'
Dear R experts,
Google is very often my friend but this time it does not !
Are you aware of an R package in which the directed causal discovery algorithm
called the Information Geometric Causal Inference (IGCI) of (Daniusis et al.,
2010) is implemented ?
Best,
Sacha
Envoyé de mon iPhone
_
https://fentechsolutions.github.io/CausalDiscoveryToolbox/html/index.html
>
> HTH,
> Eric
>
>
>> On Mon, Nov 21, 2022 at 9:00 PM varin sacha via R-help
>> wrote:
>> Dear R experts,
>>
>> Google is very often my friend but this time it does not !
>> Are you
Dear R-experts,
Here below the R code working (page 8
http://www2.uaem.mx/r-mirror/web/packages/FactoMineR/FactoMineR.pdf).
But I am trying to get all the labels (the writes) : comfort, university,
economic, world, ... smaller.
How could I do that ?
Many thanks.
library(FactoMineR)
data(chi
ersions. Probably simple ... in function definition and print(plot(res, axes
> = axes, ...)) addition could do the trick, but I am not sure.
>
> Cheers
> Petr
>
>> -Original Message-
>> From: R-help On Behalf Of varin sacha via R-
>> help
>>
1 - 100 of 187 matches
Mail list logo
