search for: eqns

Full_Name: Ross Boylan Version: 2.2.0 OS: Linux Submission from: (NULL) (65.175.48.58) \eqn{{\bf\beta}_j}{b(j)} in my .Rd file produces this error -------------------------------------------- ! Missing $ inserted. <inserted text> $ l.7 \eqn{{\bf\beta}_j}{\bf\beta}_ jnormal-bracket5bracket-normal{b(j)} -- ! Missing $ inserted. <inserted

I'm seeing various things fail when I try to next braces inside \eqn. This source \eqn{{\bf\beta}_j}{b(j)} is the vector produces this error ---------------------------------------------- [4] ! Missing $ inserted. <inserted text> $ l.258 \eqn{{\bf\beta}_j}{\bf\beta}_ j{{b(j)} is the vector of coefficients fo... I've inserted a

I'm not sure who is responsible for quantile(), but I assume they read this list. Ivan Frohne and I have produced a revision of the quantile.default() function which enables the computation of alternative sample quantile definitions. The code and .Rd file are attached. This enables the user to produce quantiles that are equivalent to those in various statistics package. There is a type

Hi, I am doing a small experiment to test if I can use MathJax ( official site: http://www.mathjax.org ) in R's html help pages (i.e. options(help_type='html')) and it seems working with some minor modifications. The screenshot (rendered by Firefox 4.0 beta and with STIX fonts) of the help page in html format with MathJax enabled is at the following url:

Surprisingly, I played around with some test code and below actually creates equations that look correct. tempmat<-matrix(10,nrow=6,ncol=6) restrictmat<-diag(6) colnames(tempmat)<-c("AUD.l1","CHF.l1","CAD.l1","GBP.l1","EUR.l1","JPY.l 1")

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 have been working with the pls procedure and have problems getting the procedure to work with matrix or frame data. I suspect the problem lies in my understanding of frames, but can't find anything in the documentation that will help. Here is what I have done: I read in an 10000 x 8 table of data, and assign the first four columns to matrix A and the second four to matrix B pls <-

I should'nt have sent the last mail so fast. Same problem with \eqn{u_j = a_j + b x + c x^2, \quad j=1, \ldots, r-1} {u[j] = a[j] + b*x + c*x^2 j = 1,\dots,r-1} I thought the problem in the first case could have to do with the use use of \mbox{} (with the braces) within the arguments of \eqn, but here there are none braces in the arguments of \eqn{}{}. Another case giving the

...umn, and perform some operation to test to see if the column is all zeros and tally up a character vector with the names of the columns to obtain the terms in equation i eqn.terms <- vector() for( v in 1:length( est$estimate ) ) { j <- attr( eval( deriv( as.formula( eqns[[i]] ), names( startvals ) ) ), "gradient" ) if( qr( j[,v] )$rank > 0 ) { eqn.terms <- rbind( eqn.terms, name <- names( est$estimate )[v] ) } } derivs[[i]] <- deriv( as.formula( eqns[[i]] ), eqn.te...

...code of linear.hypothesis. When I type "linear.hypothesis" I get: function (model, ...) UseMethod("linear.hypothesis") <environment: namespace:car> but when I type "ftest.systemfit", I do see the actual code. Why? Anyway, here are the results in more detail: eqns<-list(eq = y ~ trend+ x1 + x[,1] + x[,2] + x[,3] + x[,4] + x[,5] + x[,6] + x[,7] + x[,8] + x[,9] + x[,10] + x[,11] + x[,12] + x[,13] Rrestr10<-matrix(0,10,16);Rrestr10[1,16]=Rrestr10[2,15]=Rrestr10[3,14]=Rrestr10[4,13]=Rrestr10[5,12]=Rrestr10[6,11]=Rrestr10[7,10]=Rrestr10[8,9]=Rrestr10[9,8]=...

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 ) ) { val <- theta[i] storage.mode( val ) <- "double" assig...

...ing 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( names( theta ) ) # returns NULL ## get the values of the parameters for( i in 1:length( theta ) ) { val <- theta[i] storage.mode( val ) <- "double" assig...

Hi all, Don't know if this is a known issue, but I couldn't find anything so I report anyway. When checking eg ?qr in both RStudio and the naked R IDE, the help page is rendered incorrectly. More specifically, any use of \bold{...} is printed as is, rather than interpreted as bold. Same happens on ?svd. According to the manual Writing R Extensions, this should still be recognized. When I

Here is a modified version of qraux.Rd, an edited version of R-2.14.0/src/library/base/man/qraux.Rd This gives some details and an example for the case of pivoting. In this case, it is not true that X = QR; rather X[, pivot] = QR. It may save some other people bugs and time to have this information. Tim Hesterberg -------------------------------------------------- % File

...with OLS using 'lm', with 2SLS using 'systemfit' and it didn't accept those transformations like 'lm' does. Here's an example: this is, of course, OK: > data(kmenta) > demand <- q ~ p + d > instr <- ~ d + f > fit1 <- systemfit("2SLS", eqns=list(demand), inst=instr, data=kmenta) But, now if I'd like to estimate a model with logarithm of p as a regressor, an error occurs: > demand2 <- q ~ log(p) + d > fit2 <- systemfit("2SLS", eqns=list(demand2), inst=instr, data=kmenta) Error in log(p) : Object "p&quo...

I want to create a clean .rpm package for R built with MKL and ICC. I follow Fedora instrcutions[0] to create the package. As a base, I use the R-3.2.2.src.rpm. I am left with this error: ------------------------------------------ installing R info pages ... updating '/usr/share/info/dir' ... make[1]: Leaving directory '/home/poisonivy/rpmbuild/BUILD/R-3.2.2/doc/manual' + mv

--=-hiYzUeWcRJ/+kx41aPIZ Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: 8bit Hi, on March 10 I filed a wishlist bug report asking for the inclusion of some changes to factanal() and the associated print method. The changes were originally proposed by John Fox in 2005; they make print.factanal() display factor correlations if factanal() is called with rotation =

It seems that MASS suggest to judge on the basis of sum(residuals(mode,type="pearson"))/df.residual(mode). My question: Is there any rule of thumb of the cutpoiont value? The paper "On the Use of Corrections for Overdispersion" suggests overdispersion exists if the deviance is at least twice the number of degrees of freedom. Are there any further hints? Thanks. -- Ronggui

I'm trying to solve a system of 3 equations as part of a sub-routine in R, ie first eqn a/x-b*sqrtx+c=log(1/dx+1/e(sqrtx); snd eqn (f*y)/z-g/y-h=-log(2/x+(z/y)/(i*x) and third eqn is of the form zz=x/(j-k(z/y) where a..k inclusive are constants, x,y,z and zz are inputs. How can this be done in R? [[alternate HTML version deleted]]

'Writing R Extensions' tells you that $ needs to be escaped in Rd files outside \code and similar. So I was surprised to find that ca 80 CRAN packages have constructions like (from ISwR) \item{\code{folate}}{ a numeric vector, folate concentration ($\mu$g/l). } This does not render sensibly in non-latex conversions, and it is what we have \eqn{} for. That $ needs to be