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
Hello,
First of all, when I tried to use function mic I got an error.
mic(cbind(C, D))
Error in mic(cbind(C, D)) : could not find function "mic"
So I've changed your function myCor and all went well, with a warning
relative to BCa intervals.
myCor <- function(data, index){
mine(data[ind
Hi Rui,
Many thanks. The R code works BUT the results I get are quite weird I guess !
MIC = 0.2650
Normal 95% CI = (0.9614, 1.0398)
The MIC is not inside the confidence intervals !
Is there something wrong in the R code ?
Here is the reproducible example :
##
C=c(2,4,5,6,3,4,5,7,
You need:
myCor <- function(data, index){
mine(data[index, ])$MIC[1, 2]
}
results=boot(data = cbind(C,D), statistic = myCor, R = 2000)
boot.ci(results,type="all")
Look at the differences between:
mine(C, D)
and
mine(cbind(C, D))
The first returns a value, the second returns a symmetric mat
> On Dec 9, 2017, at 8:51 PM, ah a via R-help wrote:
>
> Hello,
>
> As I'm interested to search about the monetary transmission channel in our
> country by MSVAR model,Could you do me favor and tell me How I can run
> different types of MSVAR model (such as MSIAH(2)-VAR(2)) and finding
> impuls
Hi David, Rui,
Thanks for your precious responses. It works !
Best,
De : David L Carlson
.pt>
Cc : "r-help@r-project.org"
Envoyé le : Dimanche 10 décembre 2017 19h05
Objet : RE: [R] Confidence intervals around the MIC (Maximal information
coefficient)
Yo
I am using R(3.4.3), Win 7(extreme edition) 32 bit,
I have 10 excel files data1,xls, data2.xls .. till data10.xls.
I want to create 10 dataframes . How to do ?
regards
Parth
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing li
You could loop over the file names, read each excel file and store the
individual data frames in a list using lapply.
I prefer to read excel files with the package readxl.
The code could be along the lines of
library(readxl)
my_files <- c("file1", "file2")
lapply(my_files, read_excel)
HTH
Ulri
8 matches
Mail list logo