similar to: : Lp norm estimation

Hi all, I want to use "which.max" to identify the maximum in a 2D array/matrix and I want "argmin" and return the row and column indices. But "which.max" only works for vector... Is there any convinient way to solve this problem? Thanks a lot! [[alternative HTML version deleted]]

Dear all, I am facing the following problem in optimization: w = (d, o1, ..., op, m1, ..., mq) is a 1 + p + q vector I want to determine: w = argmin (a - d(w))' A (a - d(w)) where a is a 1xK marix, A is the covariance matrix of vector a, d(w) is a 1xK vector which parameters are functions of parameters d, o1 .. op, m1 .. mq. Is there some function to solve this problem easily? I know

Dear list, hi ! I am a R beginner and I have a function to optimize . alpha = argmin{ f(x,alpha) } I want alpha to be in [0,1]. Is there any function that can work? I use nlm() but i can't fix the domain of alpha. thanks in advance _______________________ Jiang Peng, Ph.D. Candidate Department of Mathematics & Antai college of Economics and Management Shanghai Jiao

hi all i would like to use r and winbugs in order to undertake seemingly unrelated regression. i am using the R2WINBUGS library. i just simulated a simple example (sample size is 25) in order to get the correct code. i suspect that the problem is in my definition of the prior. it wants a multivariate node. the example could be extended by including more than 1 X variable in the system of

I've found a problem when using categorical variables in lmp() from package lmPerm According to help(lmp): "This function will behave identically to lm() if the following parameters are set: perm="", seq=TRUE, center=FALSE.") But not in the case of including categorical variables: require(lmPerm) set.seed(42) testx1 <- rnorm(100,10,5) testx2 <-

Summary: The directory /usr/local/lib/R/site-library can be used for site-wide multi-user installations of R. This is even the default on Debian and Ubuntu (following a suggestion by Kurt and Fritz a decade ago over beers -- what could be better). However, R enforces wrong permissions on directories and files there (removing group-write) which effectively defeats this purpose,

Hi, I'm having trouble using the constrOptim function to generate the 9-component vector argmin of the function ELfsds: ELfsds <- function(pvechat){ LG=0 for(i in 1:9){ LG=LG+log(pvechat[i]) } return(-LG) } with accompanying gradient function: gradfunc <- function(thetavec){ g=1/(9*thetavec) return(g) } The constraints on the optimization problem are: 1 - components of

Since this has come up already: It is a good idea to re-install all packages for a minor-version increment of R, e.g. 2.5.1 -> 2.6.0 (it is major.minor.patchlevel). This is most easily done by > update.packages(checkBuilt=TRUE, ask=FALSE) However, if you don't want to do that yet, be aware that - Certain S4-using packages must be reinstalled, and using old versions can make R

Hi guys i found this encryption decryption in php and try to convert it in rails but i am unable to successfully convert it. So plz help me. I you write the whole conversion code then it will be great. PHP code is like this ================ # #/********************************************** #** #** MD5 block cypher #** #** Author..: leapinglangoor [

Dear List-members, I found a strange behaviour in the lqs function. Suppose I have the following data: y <- c(7.6, 7.7, 4.3, 5.9, 5.0, 6.5, 8.3, 8.2, 13.2, 12.6, 10.4, 10.8, 13.1, 12.3, 10.4, 10.5, 7.7, 9.5, 12.0, 12.6, 13.6, 14.1, 13.5, 11.5, 12.0, 13.0, 14.1, 15.1) x1 <- c(8.2, 7.6,, 4.6, 4.3, 5.9, 5.0, 6.5, 8.3, 10.1, 13.2, 12.6, 10.4, 10.8, 13.1, 13.3, 10.4, 10.5, 7.7, 10.0, 12.0,

Sorry, maybe it's easy but I haven't found anything useful: how can I obtain a list C that contains all the members in the list B that are not in list A? This are lists of nanes, not numbers! Thank you Gabriele Zoppoli, MD Ph.D. Fellow, Experimental and Clinical Oncology and Hematology, University of Genova, Genova, Italy Guest Researcher, LMP, NCI, NIH, Bethesda MD Work: 301-451-8575

Dear r-help, Is there any available description of the components of lqs objects found in the package "lqs"? > names(slts) [1] "crit" "sing" "coefficients" "bestone" [5] "fitted.values" "residuals" "scale" "terms" [9] "call"

Hello, the R people. I look for robust regression in R. This method is available in S, its name is rreg. Colud anyone teach me ? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the

Hi how can I find, in a vector of characters, which is the most frequent one? Thanks Gabriele Zoppoli, MD Ph.D. Fellow, Experimental and Clinical Oncology and Hematology, University of Genova, Genova, Italy Guest Researcher, LMP, NCI, NIH, Bethesda MD Work: 301-451-8575 Mobile: 301-204-5642 Email: zoppolig at mail.nih.gov

I tried building R-1.6.1 under Mac OS X 10.2.2 and experienced the following build failure: <...stuff omitted...> cc1: warning: changing search order for system directory "/usr/local/include" cc1: warning: as it has already been specified as a non-system directory gcc -bundle -flat_namespace -undefined suppress -L/usr/local/lib -o ctest.so ansari.o chisqsim.o d2x2xk.o

Is there a way of getting a list of indices or something indicating which set of points lqs when doing a regression? The $bestone, seems to return only 2 pts. no matter how good or bad the fit. I checked the archive for the past few months, but didn't find anything. thanks in advance, Mark Hall -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

Dear r-devel members, I tried to build R packages on a PC running Windows XP but experience problems. However, it is ok when there is no inst directory in a package. Any help would be appreciated. The following is an example, C:\work>R CMD check VR_7.2-19.tar.gz * checking for working latex ... OK * using log directory 'C:/work/VR.Rcheck' * using R version 2.1.1, 2005-06-20 *

hi all is there a package that undertakes subset selection but BASED ON AIC or any other information criteria. i've seen the "subselect" and the "leaps" package but i have not played around with them yet. thanx

Dear list, (this is a follow up from a previous query) Why does R CMD INSTALL write most of its messages to stderr? If it wrote to stdout, then we could capture its output within an R session when calling sink("stdout.txt", type="output") install.packages("MASS", type="source") sink() As it stands, the stderr messages can't be captured via

I would like to code the following in R: a1(b1+b2+b3) + a2(b1+b3+b4) + a3(b1+b2+b4) + a4(b1+b2+b3) or in summation notation: sum_{i=1, j\neq i}^{4} a_i * b_i I realise this is the same as: sum_{i=1, j=1}^{4} a_i * b_i - sum_{i=j} a_i * b_i would appreciate some help. Thank you. -- View this message in context: http://r.789695.n4.nabble.com/summation-coding-tp4646678.html Sent from the R