similar to: combining P values using Fisher's method

Hi, This is more of a general question than a pure R one, but I hope that is OK. I want to combine one-tailed independent p-values using the weighted version of fisher's inverse chi-square method. The unweighted version is pretty straightforward to implement. If x is a vector with p-values, then I guess that this will do for the unweighted version: statistic <- -2*sum(log(x)) comb.p <-

Hi to all is there a fisher's post hock test or fisher's combination test in R available? Maybe not -- its very easy to do it in excel ... Kind regards Knut

Is there any function to show the points like this example of SPSS? http://biostatistic.de/temp/reg.jpg The point size should represent the number of data at this point. with regards Knut Krueger

Is there any possibility to get high resolution plots in a windows xp system? I tried it with the device function png(filename = "c:/r/highresplot%d.png",pointsize=12, res=900) but when I try to set: width = 480, height = 480 or pointsize = 12, the text is not scaled in the same way as the plots. with regards Knut Krueger http://www.biostatistic.de

Hi to all I need a string like temp <- paste("m1","m2","m3",sep=",") But i must know how many items are in the string,afterwards the other option would be to use a vector temp <- c("m1","m2","m3") No problem to get the count of items but I must get afterwards the string "m1,m2,m3" No problem to build the

Hi to all, test <- data.frame("X"=c(1:4),"Y"=c(5:8),"Z"=c(8:11)) test <- test[,-2] Is there a way to specify the col name "Y" to delete instead the number? Kind regards Knut

Dear R list, For some reason I am unable to access neither search.r-project.org, nor http://finzi.psych.upenn.edu/ so I cannot search the archives for a possible answer (I Googled for this but didn't find anything). Is it possible to draw barplots using a texture instead of colors, for a black and white printer? TIA, Adrian -- Adrian Dusa Arhiva Romana de Date Sociale Bd. Schitu

If there is a helpfile open (f.e ?glm) and it is the top window, then an exception error occurs (closing RGUI) when I hit the open file button. If the helpfile is not the top window (of the RGUI) I am able to open a new script without any error. The RGUI is not closing complete there is a blank screen left which I have to close with the X Button or Taskmanager Windows XP Home - German Version

Hi to all, I am just looking for more efficient ways ;-) is there a better way instead a loop to decrease x if greater y test <- c(1,3,5,7,9) decrease if greater 1 to test2 <- c(1,2,4,6,8) Kind regards Knut

Hi to all, how could I to rotate automatically a data sheet which was imported by read.xls? x1 x2 x3 .... xn y1 1 4 7 ... xn/y1 y2 2 5 8 .... xn/y2 y3 3 6 9 ....xn/y2 yn ... ... ... Xn/Yn to y1 y2 y3 .... yn x1 1 2 3 ..... Yn/x1 x2 4 5 6 .... Yn/x2 x3 7 8 9 .... Yn/x2 xn ... ... ... ..... Yn/xn Kind regards Knut

. I am looking for a histogram or box plot with the adding normal distribution curve I think that must be possible, but I am not able to find out how to do. Regards Knut

Hi to all I found an formular of an ** ***p-Value Calculator for the Chi-Square test* *http://www.danielsoper.com/statcalc/calc11.aspx* *with the formula* *http://www.danielsoper.com/statkb/topic11.aspx* *what's the gamma function of this formula in r?* *df=5* *ch2=25.50878* *the following code does not give the result <0.001 for the values above * *p=

Hi to all, does anybody knows whether there is a package to plot those http://www.equine-science.de/temp/graph.jpg graphs. the thickness of the points and/or the lines should be represent the numbers of behaviours With kind regards Knut

I was wondering why it is possible to read long labels from the CVS files but not from the SPSS files. I did not have much time to search for the code but I found: in foreign_0.8-10 source file var.h.in > /* Definition of the max length of a short string value, generally > eight characters. */ > #define MAX_SHORT_STRING ((SIZEOF_DOUBLE)>=8 ? ((SIZEOF_DOUBLE)+1)/2*2 > : 8)

I am trying to send a manuscript to a journal. One of the figures build by R is in the right orientation and 4 are rotated clockwise 90 deg in the preview. I used the right click save to PS option and I used the command line postscript("c:/temp/fig04.eps",bg="transparent",onefile = TRUE ,pointsize=20,paper = "letter",height=8,width=8,horizontal=FALSE,family =

I searched the help list and the manual, but I do not find my mistake. Switch is working with character , with integer, but not in the third example Regards Knut count1 <- 0 test

Hi to all is there any construct to sum data=data.frame(row1=c(1,1,3,1,2,3,2,2,1,3,4,5,2,3,2,1) , row2=c(2,2,1,1,1,2,1,2,1,1,1,1,2,2,2,1) ) Means I would like to get all y of row1 if in row2 of the data.frame is an x f.e row1=3 and row2=2 so I would like to get 6 And is there another construct to get the count of pairs where row1=3 and row2=2 means the result should be

hi, is it possible to draw a plot inside another plot f.e in the upper right corner. I do not mean the possbility par(mfrow = c(2,2). Kind Regards Knut

does anybody know why there are the two warnings in the example above? Regards Knut > day_4 [1] 540 1 1 1 1 1 1 300 720 480 > day_1 [1] 438 343 1 475 1 562 500 435 1045 890 > is.vector (day_1) [1] TRUE > is.vector (day_4) [1] TRUE > wilcox.test(day_4 ,day_1,paired=TRUE,alternative="two.sided",exact=TRUE,conf.int=TRUE) Wilcoxon

Where can I found hepl about latex errors: the R CMD check tells me less, but: * checking PDF version of manual ... WARNING LaTeX errors when creating PDF version. This typically indicates Rd problems. * checking PDF version of manual without hyperrefs or index ... ERROR Re-running with no redirection of stdout/stderr. Hmm ... looks like a package Error in texi2dvi(file = file, pdf = TRUE, clean