similar to: trying to build matrix for tests...

I've trying to check for model bias. It was recommended to me that I should break the residuals into n groups and test to see if the mean is different from zero. What is the best way to generate the groups. I've generated a matrix that contains the residuals and groups, but the values all have quotes around them. 1) how do I get rid of the quotes? 2) is this an acceptable method for

R-folks, Is there a function, like aggregate, that allows users to bin values? I've got to break down a data frame into classes of 5cm (or something like it), and I only know how to do it using code like, signif <- symnum( stems$dbh, corr = FALSE, na = FALSE, cutpoints = c(0,10,20,30,40,999), symbols = c(0,10,20,30,40) )

Hello, I am trying to export the results from symnum() while maintain their readability. I tried using sink to text file and also copying and pasting but the results end up looking like this: > symnum(c5.s) bC bED bEN bLP bLS bPA bPD bPR p bbContag 1 bbED + 1 bbENN_MN + B 1 bbLPI , , , 1 bbLSI + B B , 1 bbPAFRAC , * * , * 1 bbPD , B B

Hi there, I am really new to statistics in R and statistics itself as well. My situation: I ran a lot of OLS regressions with different independent variables. (using the lm() function). After having done that, I know there is endogeneity due to omitted variables. (or perhaps due to any other reasons). And here comes the Hausman test. I know this test is used to identify endogeneity. But what I

Hello, Is there a way to put multiple data.frames or objects into a single .csv ? I know I can use sink() to do this for a text file. Is there way to do it for a .csv? e.g. for sink() sink("c5.k.s.p.92.LN.sl.vs.mean.txt") symnum(c5.k.92.LN.sl.vs.mean) symnum(c5.s.92.LN.sl.vs.mean) symnum(c5.p.92.LN.sl.vs.mean) sink() Could I put these all into a single .csv, in a similar format to

Hi, I apologize if the solution is right in front of me, but I can't find anything on how to convert a class of 'noquote' to 'matrix'. I've tried as.matrix and I've tried coercing it by 'class(x)<-matrix', but of course that didn't work. I've been using the function 'symnum' which returns an object of class 'noquote'. Here is an

I've been trying to figure out how to get the names of the parameter vector variables when inside the function that nlm calls to return the objective function value: knls <- function( theta, eqns, data, fitmethod="OLS", instr=NULL, S=NULL ) { ## print( names( theta ) ) # returns NULL ## get the values of the parameters for( i in 1:length( theta ) )

I've figured out how to use optim (barely) to estimate 2 parameter = weibull distributions. I can't get over how easy this is. What I need to = do is use a weight in the observations..... For example,=20 the tree diameters and weights are are=20 4.70 , 100 6.00, 98 7.10, 75.0 8.10, 86.3 8.60, 80.456 8.90, 20.5 9.50, 16.6 11.40, 12.657 11.80, 12.47 14.50,

I thought I had plotted a map of oregon and washington and overlayed some points using R. Is that possibe since I can't seem to find the source file I used or commands for plotting maps and points? If not, sorry for the bother... Thanks, Jeff. Jeff D. Hamann Hamann, Donald & Associates, Inc. PO Box 1421 Corvallis, Oregon USA 97339-1421 Bus. 541-753-7333 Cell. 541-740-5988 jeff_hamann

I'm having some trouble installing the add on packages. I would like to install the RMySQL_0.4-3.tar.gz add on package. I've downloaded it and unzipped the package. The INSTALL file says R CMD INSTALL [-l<alternate-dir>] RMySQL_<version>.tar.gz does this mean in the R system enter the following? > INSTALL -l /root/RMySQL_0.4-3.tar.gz (?) when I do this I get... Error:

Is it possible to write the output/results (redirect) to a latex file? Jeff. Jeff D. Hamann Hamann, Donald & Associates, Inc. PO Box 1421 Corvallis, Oregon USA 97339-1421 Bus. 541-753-7333 Cell. 541-740-5988 jeff_hamann at hamanndonald.com www.hamanndonald.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I'm trying to make sure my R package will install properly (build help files, etc) and when I try to run the rcmd check (or any rcmd option) I get the following error. Does dcf.pm stand for perl module? If not, why do I get this message? C:\systemfit\r_code>rcmd check --help Can't locate R/Dcf.pm in @INC (@INC contains: c \PROGRA~1\r\rw1041\share\perl; /usr/lib/perl5/5.6.1/cygwin-multi

I'm trying to build a vector (fitting systems of equations) and I'm having a little trouble getting the response variables from a list of formulas... bdgmodel <- lbdg ~ d1sqr + ld1 + lhg hgmodel <- lhg ~ ht2 + lht + whc inst <- ~ d1sqr + ld1 + ht2 + lht + whc systemeq <- list( bdgmodel, hgmodel ) # manually generate the y matrix... y <- cbind( matrix( lbdg ),

I've been working on a new package and I have a few questions regarding the behaviour of the nlm function. I've been (for better or worse) using the nlm function to fit a linear model without suppling the hessian or gradient attributes in the objective function. I'm curious as to why the nlm requires 31 iterations (for the linear model), and then it doesn't work when I try to add

I would like to plot lots of plots (2x2 plots) of a set of data for many sets. I can do it just fine for one set of data, but I would like to generate a small mountain of plots for others to examine. Basically, I would like to generate the same par( mcol=c(2,2) ) plot for about 1000 datasets. I've been iterating through them in the plot window, but I would like to send all the plots to a eps

I've developed a dll (win32) of a bunch of functions prototypes as: /* for building on *nix */ #ifndef WIN32 #define __stdcall /*nothing*/ #endif void __stdcall dbh_2_height( const unsigned long *func_idx, const unsigned long *metric, const unsigned long *species, const double *dbh, double *pred_height ); and have created a R interface

I'm frightfully new at developing r packages. I'm developing an econometrics package, and each time I update my beta package, I close R, delete the package dir in the library dir, rebuild the zip file, rerun R, reinstall the package and run the functions. Is there a more efficient way? Jeff. Jeff D. Hamann Hamann, Donald & Associates, Inc. PO Box 1421 Corvallis, Oregon USA 97339-1421

Dear R programming folks: I'm trying to finish off a package for non-linear simultaneous system estimation and I've been trying to figure out how to get the names of the parameter vector variables when inside the function that nlm calls to return the objective function value: knls <- function( theta, eqns, data, fitmethod="OLS", instr=NULL, S=NULL ) { ## print(

Hi, I've had a look the bug list and searched though the R documentation, email lists etc. but didn't see anything on this: when I do: summary(species.glm1, correlation = TRUE) I get a correlation matrix like this: Correlation of Coefficients: ( p I(H C pH * 1 I(pH^2) * B 1 Ca . . 1 I(Ca^2) . . B attr(,"legend") [1] 0 ` ' 0.3 `.'

I would like to produce a set of values, on the same chart, from an sql table that is structured like... species dbh ht expf DF 1.2 8.9 10.0 DF 2.4 17.3 12.4342 DF 3.1 20.9 56.76 PP 2.3 16.9 100.0 PP 12.8 97.3 40.3 PP 8.2 63.0 98.34 . . . SS blah, blah, blah... is it possible to, using a single command in the plot command to plot the different groups on the same plot or will I have to iterate