Re: [R] R-dvel [robustness Simulation study of 2 sample test on several combination of factors ]

Hi sst, You could set up your sample sizes as a matrix if you don't want all of the combinations: sample_sizes<-matrix(c(10,10,10,25,25,25,...),nrow=2) and then use one loop for the sample sizes: for(ss in 1:dim(sample_sizes)[2]) { ss1<-sample_sizes[1,ss] ss2<-sample_sizes[2,ss] then step

[R] R-dvel [robustness Simulation study of 2 sample test on several combination of factors ]

hi, i am new in this field. do favorite If I wish to conduct a simulation on the robustness of two sample test by using R language, is that any ways in writing the code? There

Re: [R] Fligner-Killeen test on binary data

Hi emeline, I think there may be a minor language problem. If you mean the "variation" rather than the "variance" in survival, you may simply want a test of proportions. Jim On Mon, Apr 4, 2016 at 7:48 PM, emeline mourocq wrote: > Hello, > > > > I investigate survival

Re: [R] multiple bar plot annotation text labelling

Hi message, What you can do is this: barpos<-barplot(axes=FALSE, ann=FALSE, horiz=TRUE, testbarplot[,2], ylab='group', xlab= '(x values)', xlim=c(0,10),space=c(1,0,0,0, 1,0,0,0)) text(testbarplot[,2],barpos, c('a', 'b', 'c', 'd', 'c', 'e','f', 'g'), pos=4) as I think you want the values

Re: [R] Test for Homoscedesticity in R Without BP Test

I have tried and got the result. Thank you every one. On Tue, Apr 5, 2016 at 12:58 AM, Achim Zeileis wrote: > On Mon, 4 Apr 2016, varin sacha via R-help wrote: > > Hi Deepak, >> >> In econometrics there is another test very often used : the white test. >> The white

Re: [R] Test for Homoscedesticity in R Without BP Test

On Mon, 4 Apr 2016, varin sacha via R-help wrote: 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

[R] Find the dataset(s) that contain(s) non-ASCII characters

Dear list, I’m maintainsing a package containing only datasets (152): http://dutangc.free.fr/pub/RRepos/web/CASdatasets-index.html When R CMD checking the package, I get the following NOTE * checking data for non-ASCII characters

Re: [R] Test for Homoscedesticity in R Without BP Test

On Mon, 4 Apr 2016, Deepak Singh wrote: Respected Sir, I am doing a project on multiple linear model fitting and in that project I have to test Homoscedesticity of errors I have google for the same and found bptest for the same but in R version 3.2.4 bp test is not available. The function is

Re: [R] Test for Homoscedesticity in R Without BP Test

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

Re: [R] Test for Homoscedesticity in R Without BP Test

You might "google Breusch Pagan test r" and find that the test is implemented in lmtest package. On 4 Apr 2016 17:28, "Deepak Singh" wrote: > Respected Sir, > I am doing a project on multiple linear model fitting and in that project I > have to test Homoscedesticity of

Re: [R] Extending the beta period for R 3.3.0 till April 25, Final release on May 3

On 02/04/2016 11:58 AM, Peter Dalgaard wrote: Due to delays in implementing an updated Windows toolchain for CRAN, R Core has found it unsafe to go with our usual 1-week beta testing period. Combined with other scheduling issues, we have decided to postpone the transition to 3.3.0 RC until

Re: [R-es] Problemas con representación de diagrama de torta

Hola, Sí, creo que estás mal interpretando cómo funciona el parámetro "x" de la función "flotaing.pie()". Los números separados por comas, *no* indican el numerador y el denominador. Si no los valores sobre el total que representa la suma de los elementos de "x". Así cuando tienes x=c(13,30) lo

Re: [R] multiple bar plot annotation text labelling

Use only plain text emails. Don't attach file types that will be stripped. See the footer at the bottom of your email for more information. Do give us the data using dput(): > dput(testbarplot) structure(list(V1 = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), V2 = c(3L, 8L, 1L, 9L, 5L, 4L, 2L, 0L), V3 =

Re: [R] Fligner-Killeen test on binary data

Thanks for all those information. I then have to find another way to test difference in variance between my two groups. -Original Message- From: peter dalgaard [mailto:pda...@gmail.com] Sent: Monday, 4 April, 2016 7:11 PM To: emeline mourocq Cc:

Re: [R] Fligner-Killeen test on binary data

Peter: Bravo! On Mon, Apr 4, 2016 at 10:10 AM, peter dalgaard wrote: > That's not an R question but a stats question, but I wouldn't do it. For one > thing: The variance of binary data is a function of the mean, so the research > question is dubious in the first place.

Re: [R] Fligner-Killeen test on binary data

That's not an R question but a stats question, but I wouldn't do it. For one thing: The variance of binary data is a function of the mean, so the research question is dubious in the first place. Secondly, the test is based on ranking and comparing absolute differences from the group median,

Re: [R] Using final sample weight in survey package

On Mon, Apr 4, 2016 at 9:35 AM, Anthony Damico wrote: > hi, probably not.. if your survey dataset has a complex design (like > clusters/strata), you need to include them in the `svydesign` call. > Coercing an incorrect survey design into a replicate-weighted design will >

Re: [R] Evaluating an expression

I think you need to spend some time with an R tutorial or two to learn how to handle such basics yourself. However, if I understand correctly, probably the easiest way to do it is by converting a matrix of 0's to a data frame -- which you shouldn't do at all if you can do your analysis directly

Re: [R] Evaluating an expression

> On Apr 4, 2016, at 9:22 AM, Janszen, Derek > wrote: > > Hi, > > I want to create a data frame similar to the following, but greatly scaled up: > df <- data.frame(aaa= c("a","b","c"), integer(3), integer(3)) > names(df)[2:3] <- paste("var",1:2,sep="")

Re: [R] Using final sample weight in survey package

hi, probably not.. if your survey dataset has a complex design (like clusters/strata), you need to include them in the `svydesign` call. coercing an incorrect survey design into a replicate-weighted design will not fix the problem of failing to account for the sampling strategy On Mon, Apr 4,

[R] Evaluating an expression

Hi, I want to create a data frame similar to the following, but greatly scaled up: df <- data.frame(aaa= c("a","b","c"), integer(3), integer(3)) names(df)[2:3] <- paste("var",1:2,sep="") which yields aaa var1 var2 1 a 0 0 2 b 0 0 3 c 0 0 I would not relish

Re: [R] question about probplot in e1071 package

Dear Thomas, Reading the probplot’s help page, it looks like it is using qqplot underneath. Thus, I think this is what you need. probplot(x, line=FALSE) #probplot(y, line=FALSE) qq.y <- qqnorm(y, plot=F) points(qq.y$y, qq.y$x) I hope this is useful for you. Best

[R] Rcmdr loading issue

Hi, I have a issue loading Rcmdr today. It is very strange, as I could get it to work last week on the same computer. I have tried restarting the computer, and re-installing R and XQuartz, but nothing changes. This is what I get: sessionInfo() R version 3.2.4 (2016-03-10)

[R] Fligner-Killeen test on binary data

Hello, I investigate survival until the following year (0,1) and I wish to test if the variance in survival for two or more groups are significantly different from each other. I read that the Fligner-Killeen test is a non-parametric test which is very robust against departures from

[R] multiple bar plot annotation text labelling

Readers, The attempt is to create a bar plot with text labels adjacent to each datum value. Data file: 1,3,A 1,8,B 1,1,C 1,9,D 2,5,C 2,4,E 2,2,F 2,0,G testbarplot<-read.csv('data1.csv', header=FALSE) barplot(axes=FALSE, ann=FALSE, horiz=TRUE, testbarplot[,2], ylab= 'group', xlab= '(x

[R] R cases on predictive maintenance

Hi Team, Can you please suggest me some good cases where we can use R programming to tackle predictive maintenance problems Many thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,

[R] Using final sample weight in survey package

I have the final sample weight (expansion factor) from a socieconomic survey. I don't know the exact design used in the study ( (probably is a stratified two-stage design). To illustrate my problem I will use the next dataset which have a sample weight (but the design is not specified) and

[R] Question about function DrawDensity3D {VecStatGraphs3D}

Hi, Here is the function DrawDensity3D in package VecStatGraphs3D. My question is: if we use more layers than one, could we change the function in a way that in the final plot only the outmost layer is drawn (the inner layers omitted)? Best regards, Atte Tenkanen function (vectors, Div =

[R] Test for Homoscedesticity in R Without BP Test

Respected Sir, I am doing a project on multiple linear model fitting and in that project I have to test Homoscedesticity of errors I have google for the same and found bptest for the same but in R version 3.2.4 bp test is not available. So please suggest me a test on homoscedesticity ASAP as we

Re: [R] question about probplot in e1071 package

Luisfo, Thank you so much! That does what I need. Best regards, Tom On Mon, Apr 4, 2016 at 10:51 AM, Luisfo Chiroque wrote: > Dear Thomas, > > Reading the probplot’s help page, it looks like it is using qqplot > underneath. > Thus, I think this is what you need. >

[R-es] tentativa de aplicar una no-función

Hola a todos, Tengo este problema con la nueva versión de ggplot. http://stackoverflow.com/questions/36160728/error-calling-a-function-in-r-attempt-to-apply-non-function Alguien sabe cómo reformular la función para que no dé error? Muchas gracias! Aquest correu electrònic i els annexos

[R] question about probplot in e1071 package

Hello! I am using probplot in the e1071 package and want to do something like the following, only with the the 2nd plot overlaying the first. I can't seem to make it work. Any suggestions? *library(e1071) **x <- rnorm(100, mean=5)* *y <- rnorm(100, mean=3)* *probplot(x, line=FALSE) *

Re: [R] minimal attributes to get se.fit

> Does anybody know what are the attributes of a glm fit object that will allow > the "predict call" to produce an se.fit? > > > > I believe the covariance matrix of the coefficients is all that should be > needed > and that is quite small. However, the covariance matrix is not an attribute

[R] minimal attributes to get se.fit

Does anybody know what are the attributes of a glm fit object that will allow the "predict call" to produce an se.fit? I am deleting most of the attributes as the size of the final object is 5Gb and I want to reduce it to under 20Mb, but that causes as error when I ask for an se.fit .

Re: [R] Using R for cURL commands,Message-ID:

On 04/02/2016 05:00 AM, r-help-requ...@r-project.org wrote: Hello, I'm looking for a way in which R can make my live easier. Currently i'm using R convert data from a dataframe to json's and then sending these json's to a rest api using a curl command in the terminal (i'm on a mac). I've

Re: [R] before-after control-impact analysis with R

I searched for analyze BACI data in R and found a couple things that you might find helpful. http://people.stat.sfu.ca/~cschwarz/Stat-650/Notes/MyPrograms/BACI/BACITalk-2012-10-15-UBC/baci.pdf http://people.stat.sfu.ca/~cschwarz/Stat-650/Notes/PDFbigbook-R/R-part013.pdf Jean On Sun, Apr 3,

Re: [R] apply mean function to a subset of data

Thanks David It works perfectly Pedro Sent from my iPhone > On Apr 3, 2016, at 17:44, David L Carlson wrote: > > Here are several ways to get there, but your original loop is fine once it is > corrected: > >> for (i in 1:2) smean[i] <-

[R] Does this code execute the bagging correctly ?!

Hello the code : set.seed(10) y<-c(1:1000) x1<-c(1:1000)*runif(1000,min=0,max=2) x2<-c(1:1000)*runif(1000,min=0,max=2) x3<-c(1:1000)*runif(1000,min=0,max=2) lm_fit<-lm(y~x1+x2+x3) summary(lm_fit) set.seed(10) all_data<-data.frame(y,x1,x2,x3) positions <-

Re: [R] use one way ANOVA to select genes

Hi Your message is scrambled by using HTML post and if I deciphered it correctly you just want to use anova for some data. What is preventing you from doing it? Cheers Petr -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of hehsham alpukhity via R-help