Re: [R] dependent p.values in R

2016-07-16 Thread Fernando Marmolejo Ramos
Fernando Marmolejo-Ramos Postdoctoral Fellow Gösta Ekman Laboratory Department of Psychology Stockholm University Frescati Hagväg 9A, Stockholm 114 19 Sweden ph = +46 08-16 46 07 website = http://sites.google.com/site/fernandomarmolejoramos

Re: [R] dependent p.values in R

2016-07-10 Thread Fernando Marmolejo Ramos
could simply average those p,values but to my knowledge that approach is biased so i thought, in the same way there is a method to combine independent pvalues (e.g. stouffer method); is there a way to combine dependent pvalues? best f Fernando Marmolejo-Ramos

[R] dependent p.values

2016-07-09 Thread Fernando Marmolejo Ramos
hi all does any one know a method to combine dependent p.values? best Fernando Marmolejo-Ramos Postdoctoral Fellow G�sta Ekman Laboratory Department of Psychology Stockholm University Frescati Hagv�g 9A, Stockholm 114 19 Sweden ph = +46 08-16 46 07

[R] g2 test...

2013-07-04 Thread Fernando Marmolejo Ramos
dear R users is it appropriate to use a Log likelihood ratio (G-test) test of independence when dealing with repeated categorical responses (e.g. 2 by 2 table) instead of the McNemar test? thanks fer [[alternative HTML version deleted]] __

[R] calculating the mode in R...

2011-03-29 Thread Fernando Marmolejo Ramos
Dear R users I’m aware that the package “modest” is useful to find the mode in an array. However, I’d like to know if someone has translated the “mode” function built-in in MATLAB into R language. This function finds the most frequent value in an array

[R] confidence intervals around p-values

2010-09-09 Thread Fernando Marmolejo Ramos
Dear all I wonder if anyone has heard of confidence intervals around p-values... Any pointer would be highly appreciated. Best Fer __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] R-help Digest, Vol 91, Issue 9

2010-09-09 Thread Fernando Marmolejo Ramos
dear all I wonder if anyone has heard of confidence intervals around p-values... Any pointer would be highly appreciated. Best Fer __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] simulating Likert-type data

2009-06-15 Thread Fernando Marmolejo Ramos
Dear R members Could someone indicate me how to simulate Likert-type data using the GAMLSS library and how to fit that data? Let’s say, 50 random numbers in a variable that goes from -4 to +4 in steps of 1 (i.e., -4, -3, -2, -1, 0, 1, 2, 3, 4) with a mean of 2 and SD of 0.15 and belonging to a

[R] transformation and outliers

2009-05-21 Thread Fernando Marmolejo Ramos
Dear R people I ask again… 1. Is there a published reference presenting the normal score transformation? Is there a published paper (in any field) using that transformation in the analysis of data? And this is a new question… 2. The “outliers” library has a function called

[R] normal score transformation...

2009-05-18 Thread Fernando Marmolejo Ramos
Dear R members I’m interested in reviewing some data transformations. I’ve found literature regarding Box-Cox, logarithmic and square-root transformations. I know there is a library called “coin” which contains a transformation called normal score transformation (see “normal_trafo” function).

Re: [R] Fitting weibull and exponential distributions to left censoring data

2008-10-28 Thread Fernando Marmolejo Ramos
Dear Barja Have you looked at gamlss y gamlss.dist libaries? There some functions called WEI, WEI2, and WEI3 Cheers, Fer __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] multiple visualisation of a bootstrap

2008-10-24 Thread Fernando Marmolejo Ramos
Dear R team Jim Lemon recently said regarding new improvements in PrettyR and plotrix: “Remember, it is your whingeing and moaning that have helped to make these packages what they are today.” So, this is my whinge: out there is a freeware called Vista, which runs, among other things, bootstrap

Re: [R] Interpretation of t.test results

2008-10-23 Thread Fernando Marmolejo Ramos
Dear Omar Your t.test is telling that the 2 groups are not different. Now it is not clear if you want to compute the CI for the difference or for each group. If the former, then the outputs provides you with that information. It’d be clearer if you provide the full syntax of your t.test. For

[R] A line for SD, colouring simplevioplots, limiting Y axis in barplots, axis in vioplot

2008-10-23 Thread Fernando Marmolejo Ramos
Dear R people I have 4 questions about graphics. QUESTION 1 2 --- I have the following data, which is also set as a dataframe. Then I plot it using a violin plot… 1. How can a draw a line representing 2 SD for each group? 2. How can I colour the simple violin plots in different

[R] about granova library

2008-10-16 Thread Fernando Marmolejo Ramos
Dear all Recently the granova package was launched. I installed but after when I invoked it in R it requested for other libraries. They were downloaded and install automatically. I tried to run the example syntax of “granova.1w” and “granova.2w” but two things happened: i) either a file called

Re: [R] dixon test

2008-08-12 Thread Fernando Marmolejo-Ramos
hi giov about the dixon test... i just run a simple test with a sample of 40 and I got: Error in dixon.test(x) : Sample size must be in range 3-30 So it seems that most of the test in the outliers package are designed for small samples. See also the Rnews article published in May 2006 (vol

[R] about the 95%CI around the median...

2008-08-04 Thread Fernando Marmolejo Ramos
Dear people I've learnt that by using the boxplot.stats command in the grDevices library I can get the 5-number summaries of a boxplot, plus other important information, like the confidence interval around the median. I'm interested in knowing the actual formula to used in that package to

Re: [R] adding the mean and standard deviation to boxplots

2008-07-23 Thread Fernando Marmolejo-Ramos
Dear users This is a message I was directing to Harold Baize but because I pressed the wrong button the message got lost g!!! So I’m doing it all over again: Lets suppose I have three batches of data: a - rnorm(50,2500,300) b - rnorm(50,3500,250) c - rnorm(50,4000,200) # Now I want to