similar to: R-source code of a function

I need read a table in a string with special format. I used read.csv and textConnection function. But i am confuse about textConnection by follow code. case A: It is OK£¡ str0 <- '{"abc",{"def","X,1&Y,2&Z,3"}}' str1 <- strsplit(str0,'"')[[1]][6] str2 <- gsub("&","\n", str1) con <-

I am needing some expertise with regard to the S-Plus command lmRobMM and its R counterpart rlm(formula,data,method="MM") I have used lmRobMM(formula,data) in S-Plus on the Stackloss data and obtained for my residuals 6.217777 1.150717 6.427946 8.174019 -0.6713005 -1.248641 -0.4236203 0.5763797 -1.057899 0.3593823 11 12 13 14 15 16

I estimate two competing simple regression models, A and B where the LHS is the same in both cases but the predictor is different ( I handle the intercept issue based on other postings I have seen ). I estimate the two models on a weekly basis over 24 weeks. So, I end up with 24 RSquaredAs and 24 RsquaredBs, so essentally 2 time series of Rsquareds. This doesn't have to be necessarily thought

Hello, I'm trying to annotate my plots nicely and am running into trouble. This example contains two problems: a) the 'text()' arguments do not show the conditional behavior I'm trying to give them. I try to test for the intercept of my regression and reformat the output accordingly ('+ intercept' in the >= 0 case and '- sqrt(intercept^2)' in the other case),

Dear all, I have a question about the knowing for which row I have the max value of one of my variables. I calculated the Rsquared for different columns and made a list to gather them. I unlisted this list to create a vector with this values. I want to know for which column I have the max value of Rsquared. The columns were always named in the same way. They always start with results4$depth_

hello mailing list! i still consider myself an R beginner, so please bear with me if my questions seems strange. i'm in the field of biology, and have done consecutive hydraulic conductivity measurements in three parallels ("Sample"), resulting in three sets of conductivity values ("PLC" for percent loss of conductivity, relative to 100%) at multiple pressures

Hi, I had a problem with LLVM not emitting local variable info, even though I had calls to llvm.dbg.define. After some tracking I found that the llvm.dbg.declare (and probably value too) have to have a !dbg !nr after them to get emitted at all. maybe someone can adjust http://llvm.org/docs/SourceLevelDebugging.html#format_common_declare that the two llvm.dbg. functions need a !dbg line info

I calculate an anova test in the following way: expdata<-read.table("/home/dorien/UA/meta-music/optimuse/optimuse1-build-desktop/results/results_processedCP", header=TRUE)

Hello, I'm with a conceptual doubt regarding Rsquared of both lm() and postResample(library caret). I've got a multiple regression linear model (lets say mlr) with anR² value of 67.52%. Then I use this model pro make predictions with predict() function using the same data as input , that is, use the generated model to predict the value associated with data that I used as input. Next, if

Hi Mark, I don't know wether you recived a sufficient reply or not, so here are my comments to your problem. Supressing the constant term in a regression model will probably lead to a violation of the classical assumptions for this model. From the OLS normal equations (in matrix notation) (1) (X'X)b=X'y and the definition of the OLS residuals (2) e = y-Xb you get - by

Hello, I'm doing some linear regression: >lm<-lm(osas~alp,data) >summary(lm) However, the Rsquared in the output of summary() is not the same as the "standard" Rsquared calculated by spreadsheets, and outlined in statistical guidebooks, being SSR/SSTO. The output says "multiple Rsquared", but it is no multiple regression... What's the difference? Thanks,

Dear all, I have been using rlm() for robust regression. Could someone please suggest an appropriate measure of goodness-of-fit [1]? All I've found after trawling the web, literature databases, and previous r-help posts, is the "robust R^2" on pp. 362-363 of the S-plus manual, which is available at http://web.mit.edu/afs/athena/software/splus_v7.0/www/statman1.pdf (7.57 MB)

Hi, I'm using glm function to do logistic regression and now I want to know if it exists a kind of R-squared with this function in order to check the model. Thank you.

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

All, Happy World Cup and Wimbledon. This morning finds me with the first of my many daily questions. I am running a ks.test on residuals obtained from a regression model. I use this code: > ks.test(Year5.lm$residuals,pnorm) and obtain this output One-sample Kolmogorov-Smirnov test data: Year5.lm$residuals D = 0.7196, p-value < 2.2e-16 alternative hypothesis: two.sided Warning

Dear R-helpers, I am developing a Mixed-Effects model for a study of immunoassays using 'lme4' library. The study design is as follows: 10 samples were run using 7 different immunoassays, 3 times each, in duplicates. Data attached. I have developed the following model: c.lme <- lmer(Result~SPL + (SPL|Assay/Run) -1, data=data) This model has excellent predictions - the Rsquared of

Hello, In my lab we use a four parameter logistic fit model for our ELISA data (absorbance values). We are currently testing the use of different solvents and need to find a way to add a correlation value (such as an R squared or something similar) so we can test different solvents in making this standard curve. We currently use the drc package and this is our script for the 4 parameter:

I'm attempting to modify a parser generator to emit LLVM code instead of C. So far the experience has been trivial, but I am now running into an error regarding labels that I can't seem to solve. Situation 1: A label is used immediately after a void function call (l6 in this case): <snip> %tmp26 = load i32* @yybegin, align 4 %tmp27 = load i32* @yyend, align 4 call void

Hello, I have a question concerning ?for loops? on multiple columns. I made 91 columns with results (all made together with a for loop) and I want to us lm to fit the model. I want to compare the results of all these calculated columns (91) with one column with observed values. I use the function lm to fit the model and calculate r.squared. I manage to do this for each column separately: For

I have two variables, minutereturnsa which can be thought of as my independent variable and minutereturnsb which can be thought of as my dependent variable. When I run correlations on the two variables, depending on which of the three methods I use, I get values of between -.15 through -.19. Then, when I do a regression, I get an rsquared of .004 which is more in line with my intuition. In other