similar to: rcorr p-values for pearson's correlation coefficients

Hello All, I was wondering if it's possible to reshape data from long to wide in R without using a "timevar". I've pasted some sample data below along with some code. The data are sorted by Subject and Drug. I want to transpose the Drug variable into multiple columns in alphabetical order. My data have a variable called "RowNo" that functions almost like a

Hi, I'm a new R user and I'm having trouble with the read.csv command. It somehow treats the first column as a row name field even though it's not a row name. there are no missing columns/entries and i'm not sure how to resolve this. the format of my data is A, B, C, D,......(3984 columns) 12, 13, 41,......(all numeric) it either treats column A as rownames or if I explicitly

Dear R-helpers, I recently switched from SAS to R, in order to model the occurrence of rare events through logistic regression. Is there a package available in R to calculate the Goodman-Kruskal Gamma? After searching a bit I found a function "rcorr.cens" which should do the job, but it is not clear to me how to define the input vectors? Is "x" a vector with the fitted

Hi, I try to export the outputs of rcorr into excel. but I got error message,"cannot coerce class "rcorr" into a data.frame". Actually i just need export part of results of this analysis,e.g. p-values or stat-values. Does anyone have sort of exprience before or you can help on how to export subset of results of analysis? Many Thanks! Xin

Hello R-User! I have a data.frame with 82 variables (columns) and 290 rows. The variables are set to classes factor, ordered factor and numeric. I used the following code Matrix.My.data<-as.matrix(Df.My.Data[2:82]) Matrix.My.data.rcorr<-rcorr(Matrix.My.data, type="spearman") and got the following message: Fehler in storage.mode(x) <- if (.R.) "double" else

I am using the rcorr.cens function from the Hmisc package and the time-dependent ROC curve obtained using tdrocc in the survcomp package. I understand that the C statistic from rcorr.cens has to be subtracted from 1 if high values of the risk variable lower survival. Given that I wonder what the connection is between that C statistic and the AUC from the tdrocc object. If they are substantially

Hi, I've just run an rcorr on some data in Spearman's mode and it's just produced the following values; [,1] [,2] [1,] 1.00 -0.55 [2,] -0.55 1.00 n= 46 P [,1] [,2] [1,] 0 [2,] 0 I presume this means the p-value is lower than 0.00005, but is there any way of increasing the number of significant figures used? How should I interpret this value? Cheers Jim

Dear List, I would like to calculate the Goodman-Kruskal gamma for the predicted classes obtained from an ordinal regression model using lrm in the Design package. I couldn't find a way to get gamma for predicted values in Design so have found previous positings suggesting to use : Rcorr.cens(x, S outx = TRUE) in the Hmisc package My question is, will this work for predicted vs

Dear r-helpers! I want to classify my vegetation data with hierachical cluster analysis. My Dataset consist of Abundance-Values (Braun-Blanquet ordinal scale; ranked) for each plant species and relev?. I found a lot of r-packages dealing with cluster analysis, but none of them is able to calculate a distance measure for ranked data. Podani recommends the use of Goodman and Kruskals' Gamma for

Dear list, I get the following message when I use rcorr in library "Hmisc" ------------------------------------------------------ > rcorr(lskPox0t30, type=c("spearman")) Error in "[<-.data.frame"(*tmp*, is.na(x), value = 1e+30) : matrix subscripts not allowed in replacement ------------------------------------------------------ I do not understand

Dear all ?Hmisc::rcorr states that it takes as main argument "a numeric matrix". But is it normal that it fails in such an ugly way on a data frame? (See below.) If the function didn't attempt any conversion to a matrix, I would have expected it to state that in the error message that it didn't accept 'data.frame' objects in its input. Also, I vaguely remember having used

Hi, I'm trying to calculate 95% confidence interval of C statistic of logistic regression model using rcorr.cens in rms package. I wrote a brief function for this purpose as the followings; CstatisticCI <- function(x) # x is object of rcorr.cens. { se <- x["S.D."]/sqrt(x["n"]) Low95 <- x["C Index"] - 1.96*se Upper95 <- x["C

Dear R-help, This is an example in the {Hmisc} manual under rcorr.cens function: > set.seed(1) > x <- round(rnorm(200)) > y <- rnorm(200) > round(rcorr.cens(x, y, outx=F),4) C Index Dxy S.D. n missing uncensored Relevant Pairs Concordant Uncertain 0.4831 -0.0338 0.0462 200.0000

Hi all, Sorry about posting a really novice question. I was able to run rcorr after converting the list to a matrix by your help. I'm though wondering if there is any way to find out an exact p value as the output only gave me 0 for P value as shown below. I've added options(digits=10), which doesn't seem to help at all. Any help would be appreciated. P D Prime T

Hello, Also, the other file, NPA.csv, is not in tabular form. Can you please reformat it? Rui Barradas Citando ruipbarradas at sapo.pt: > Hello, > > Please keep this on the list, always cc r-help. > One of the files in your attachment is empty: > > y <- read.csv(file.choose("GT.csv")) > Error in read.table(file = file, header = header, sep = sep, quote =

Hi all, I try to receive an ipv6 address for my PPP link via autoconf (against Cisco machine), but there is a problem I'm unable to solve so far. I can see cisco's Router-Advertisement containing prefix etc., but no IP adress is assigned to tun0 interface. Anyone succesful in IPv6 over dial-up PPP connection? When I manually assign IPv6 address, connection works. My system is

Hello, Please keep this on the list, always cc r-help. One of the files in your attachment is empty: y <- read.csv(file.choose("GT.csv")) Error in read.table(file = file, header = header, sep = sep, quote = quote,? : ? no lines available in input Rui Barradas ? Citando Chaitanya Ganne <Chaitanya.Ganne at jefferson.edu>: > Thank you so much for your input. > > I am

Hi! All, I am trying to fetch rows from a data frame which matches to first 2 columns of another data frame. Here is the example what I am trying to do: > ptable=read.table(file="All.txt",header=T,sep="\t") > ptable=as.matrix(ptable) > dim(ptable) [1] 9275 6 > head(ptable) Gene1 Gene2 PCC PCC3 PCC23 PCC123 [1,]

Dear colleagues, This is not strictly a R question, but I hope it is ok to ask on the list. I fed a vector of p-values from about 20 million anova tests to the package q-value and obtained this output: > qsummary(asso_p.qvalue) Call: qvalue(p = asso_p.vec) pi0: 1 Cumulative number of significant calls: <1e-04 <0.001 <0.01 <0.025 <0.05 <0.1 <1

We need to monitorate the E1 with nagios, somebody did this? any ideia? Thanks in advance! -- Atenciosamente, ------------------------------------------------------- Dario Quiroz (71) 9275-9080 gtalk: darioquiroz at gmail.com ------------------------------------------------------- -------------- next part -------------- An HTML attachment was