search for: deter088

Hello, Simple question that I am stuck on and can't seem to find an answer in the help files currently. I have a list which contains repeated ID's. I would like to have R count the number of ID's. For example: ID=c(1,1,1,1,2,2,2,2,3,3,3,3) as.data.frame(ID) Clearly, there are 3 groups. How would I have R give me the summary: ID 3 Many thanks, Charles [[alternative HTML

Greetings, I have been exploring the use of the caret package to conduct some plsda modeling. Previously, I have come across methods that result in a R2 and Q2 for the model. Using the 'iris' data set, I wanted to see if I could accomplish this with the caret package. I use the following code: library(caret) data(iris) #needed to convert to numeric in order to do regression #I

A general question that I have been pursuing for some time but have set aside. When finishing some analysis, I can have multiple matrices that have specific column names. Ideally, I would like to combine these separate matrices for a final output as a csv file. A generic example: Matrix 1 var1A var1B var1C x x x x x

Hello, I am trying to reformat some data so that it is organized by group in the columns. The data currently looks like this: group X3.Hydroxybutyrate X3.Hydroxyisovalerate ADP 347 4 4e-04 3e-04 5e-04 353 3 5e-04 3e-04 6e-04 359 4 4e-04 3e-04

Hello R users, This is more of a convenience question that I hope others might find useful if there is a better answer. I work with large datasets that requires multiple parsing stages for different analysis. For example, compare group 3 vs. group 4. A more complicated comparison would be time B in group 3 of group L with B in group 4 of group L. I normally subset each group with the

Greetings R users, I have a curious problem. I read in a csv file (subset shown below) as normal data=read.table("C:/Users/Chaz/Desktop/test.csv",sep=",",header=TRUE, na.strings=".") However, the numbers from the dataset are not registered as numeric: is.numeric(data$Mesh) [1] FALSE When I try as.numeric, it converts all the values to different integers. This

This is more of a general question without data. After doing 'survdiff', from the 'survival' package, on strata including four groups (so 4 curves on a Kaplan Meier curve) you get a chi squared p-value whether to reject the null hypothesis or not. Is there a method to followup with pairwise testing on the respective groups? I have searched the library but have come up with

Generic question... I am familiar with generic power calculations in R, however a lot of the data I primarily work with is multivariate. Is there any package/function that you would recommend to conduct such power analysis? Any recommendations would be appreciated. Thank you for your time, Charles [[alternative HTML version deleted]]

Hello, I am trying to set up a loop that can run the survdiff function with the ultimate goal to generate a csv file with the p-values reported. However, whenever I try a loop I get an error such as "invalid type (list) for variable 'survival_data_variables[i]". This is a subset of my data: structure(list(time = c(1.51666666666667, 72, 72, 25.7833333333333, 72, 72, 72, 72, 72,

To R users, I am trying to use cut2 function from the 'Hmisc' library. However, when I try and run the function on the following variable, I get an error message (displayed below). I suspect it is because of the NA but I have no idea how to address the error. Many thanks to any insights. structure(list(var1 = c(97, 97, 98, 98, 97, 99, 97, 98, 99, 98, 99, 98, 98, 97, 97, 98, 99, 98,

Hello everyone, I'm having some trouble interpreting the results of a Wilcoxon (Mann-Whitney U) test. Hope you can help. This is the R script that I am running: a <- c(1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 1, 5, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 3, 1, 1) b <- c(1, 2, 1, 1, 2, 3, 2, 2, 1, 2, 1, 1, 1, 2) wilcox.test(a, b, alternative="t", mu=0, exact=FALSE, paired=FALSE) #1st

Greetings R users, I am trying to renumber my groups within the file shown below. The groups are currently set as 8,9,10,etc. I would like to renumber this as 1,2,3,etc. I have searched the help files and only come across using the rownames to renumber the values but I need to match values. Any assistance is always appreciated, Regards, Charles structure(list(Group = structure(c(1L, 1L, 1L,

Greetings R users, I have been working on running plsda and I would like to have the R2 and Q2 values. I know the function R2 from the 'pls' package will generate both R2 and Q2 but they are for each separate class. Is there a way to get the cumulative R2 and Q2 for the whole model? R2(pls.new, estimate="all") Response: B (Intercept) 1 comps 2 comps train

Greetings R users, I have been hoping someone would be familiar with this topic. I understand fully everything on this list is from the good graces of those who wish to help. Thanks to those who have helped in multiple circumstances. However, I wanted to post this question once more. I hope someone will read it and perhaps have some additional insight into what I may be able to do. Thank

Greetings R users, My interest in the Q2cum score comes my endeavor to replicate SIMCAP PLS-DA analysis in R. I use the exact same dataset. After doing the analysis in R, I can get the exact same R2Ycum. However, the Q2cum is significantly off. Adding the Q2cum of the 1st and 2nd component comes close but that seems unconvincing and I don't understand why the 3rd component Q2cum would be

Greetings R users, My interest in the Q2cum score comes my endeavor to replicate SIMCAP PLS-DA analysis in R.? I use the exact same dataset.? After doing the analysis in R, I can get the exact same R2Ycum.? However, the Q2cum is significantly off.? Adding the Q2cum of the 1st and 2nd component comes close but that seems unconvincing and I don't understand why the 3rd component Q2cum would be

Greetings, I am working on modifying some code from the nlme package. I have had many discussions on the mixed models mailing list and have been directed to simply 'hack' the source code to have the degrees of freedom generated by one function to use in the output of another function that doesn't generate them. My current holdup is an error regarding a .c file called

Hello, Does anyone on this list know what inner_perc_table is or where it is typically found? I am trying to modify some source code and it is used with the .C() function. When I try and run it, it states that 'inner_perc_table is not found'. It is only called in such a way and isn't defined at any point in the previous code in which the function works. Should this be in some other

Greetings, I am exploring some random forest analysis methods and have come upon one aspect I don't fully understand from any manual. The code of interest is as follows from the randomForest package: myiris=cbind(iris[1:4], matrix(runif(508*nrow(iris)),nrow(iris),508)) This would be following by the rfcv() function for cross-validation but I am confused about the former syntax. My

Greetings, I am familiar with the function cite('packageName') which provides the output generated from the DESCRIPTION file. In most cases this is sufficient but I was wondering if there are contributing authors (in addition to the primary) also listed on the CRAN page. Is there a proper way to account for them or are they generally not listed? Regards, -- Charles Determan