search for: diml

Hi My code is as following: #include <R.h> #include <Rinternals.h> //* the Projector part *// void Projector(double *L, int *dimL, double *G, int *dimG, double *W, int *dimW, int *xymod, int *dimxy, double *modif, int *dimif, double *Lsum) { ...} //* the interface part *// #define getDim(A) INTEGER(coerceVector(getAttrib(A,R_DimSymbol), INTSXP)) SEXP Projector5(SEXP L, SEXP G, SEXP W, SEXP xymod, SEXP modif) { //* diges...

...IC(G));nprot++; PROTECT(W=AS_NUMERIC(W));nprot++; PROTECT(xymod=AS_INTEGER(xymod));nprot++; PROTECT(modif=AS_NUMERIC(modif));nprot++; double *lptr; lptr=REAL(L); double *gptr; gptr=REAL(G); double *wptr; wptr=REAL(W); int *xyptr; xyptr=INTEGER(xymod); double *ifptr; ifptr=REAL(modif); int *dimL; dimL=INTEGER(GET_DIM(L)); int *dimG; dimG=INTEGER(GET_DIM(G)); int *dimW; dimW=INTEGER(GET_DIM(W)); int *dimxy; dimxy=INTEGER(GET_DIM(xymod)); int *dimif; dimif=INTEGER(GET_DIM(modif)); SEXP ans; PROTECT(ans=allocMatrix(REALSXP, dimG[0], dimL[1])); nprot++; double *ansptr; ansptr=REAL(ans);...

>From one of my students' simulations: m <- matrix(1, 0, 0) # 1 to force numeric not logical eigen(m) and segfault in TRED2 in src/appl/eigen.f Easy to fix, but I wonder what else might have been overlooked? (svd is protected). --please do not edit the information below-- Version: platform = sparc-sun-solaris2.7 arch = sparc os = solaris2.7 system = sparc, solaris2.7 status =

Suppose I have an n-diml array A and I want to extract the first "row" -- ie all elements A[1, ...] Interactively if I know 'n' I can write A[1,,,,,] with (n-1) commas. How do I do the same more generally, eg in a script? (I can think of doing this by converting A to a vector then extracting the appro...

Why the warning messages (2:4)? > x <- rexp(1000,0.2) > fitdistr(x,"exponential",list(rate=1)) rate 0.219824219 (0.006951308) Warning messages: 1: one-diml optimization by Nelder-Mead is unreliable: use optimize in: optim(start, mylogfn, x = x, hessian = TRUE, ...) 2: NaNs produced in: dexp(x, 1/rate, log) 3: NaNs produced in: dexp(x, 1/rate, log) 4: NaNs produced in: dexp(x, 1/rate, log) and with respect to this function, it is curious that only...

Hi dear, I want to estimate d.f for Chi-squared distribution: est.chi[i,]<-c( fitdistr(as.numeric(data2[,i]),"chi-squared",start=list(df=1))$estimate)Warning message:In optim(x = c(7.86755, 7.50852, 7.86342, 7.70589, 7.70153, 7.58272, : one-diml optimization by Nelder-Mead is unreliable: use "Brent" or optimize() directly Who can help me to solve this problem? Best whishes, Soheila [[alternative HTML version deleted]]

i am running one script for many files, Script file atttached:- http://r.789695.n4.nabble.com/file/n4556770/AR_anarkali_aic.r AR_anarkali_aic.r & warnings are given as below what is meaning of these warnings? Warning messages: 1: In optim(coef, err, gr = NULL, hessian = TRUE, ...) : one-diml optimization by Nelder-Mead is unreliable: use "Brent" or optimize() directly is there anything wrong with my code or my approach to calculate AIC values ? how i know due to which function warnings is generated? if i want to print warnings with file name then ...........? -- View this m...

...essage: bounds can only be used with method L-BFGS-B in: optim(7, fn = fun, u = 9) while R> fun<-function(x,y) (x-y)^2 R> optim(7,fn=fun,y=9) $par [1] 8.999854 $value [1] 2.145767e-08 $counts function gradient 28 NA $convergence [1] 0 $message NULL Warning message: one-diml optimization by Nelder-Mead is unreliable: use optimize in: optim(7, fn = fun, y = 9) -Christina

Hello, I want to estimate the parameters of a binomial distributed rv using MLE. Other distributions will follow. The equation system to solve is not very complex, but I've never done such work in R and don't have any idea how to start... The system is: (1) n*P = X (2) [sum {from j=0 to J-1} Y{j} /(n-j)] = -n * ln (1-X / n) where * only X is given (empirical mean)

...The function was originally written for an older version of SPlus. Can anyone give me some insights into where the problem is? Thanks R 2.4.1 on MAC OSX 2mb ram Mark Grant markg at uic.edu > attach(Aspirin.frame) > hblm(Diff ~ 1, s = SE) Error in solve.default(R, rinv) : 'a' is 0-diml > traceback() 6: .Call("La_dgesv", a, b, tol, PACKAGE = "base") 5: solve.default(R, rinv) 4: solve(R, rinv) 3: summary.blm(fit) 2: eb.calc(rho[i], X, Y, s.e., df.se, corrs, prior, ...) 1: hblm(Diff ~ 1, s = SE) > > hblm function(formula, s.e., df.se = Inf, corrs = F,...

Dear List, I am not sure how should i optimize a log-likelihood numerically: Here is a Text book example from Statistical Inference by George Casella, 2nd Edition Casella and Berger, Roger L. Berger (2002, pp. 355, ex. 7.4 # 7.2.b): data = x = c(20.0, 23.9, 20.9, 23.8, 25.0, 24.0, 21.7, 23.8, 22.8, 23.1, 23.1, 23.5, 23.0, 23.0) n <- length(x) # likelihood from a 2 parameter Gamma(alpha,

...d to down columns. Now par(mfrow/mfcol = c(nr, nc)) followed by par(mfg=c(i,j)) is the preferred style. o plot( <aov.object> ) failed after 2nd plot, PR#279 o Several patches for systems with long != int. o approx failed if inputs contained NAs. o weighted lm fits with 1 diml x failed if any weight was zero. o lm.influence failed when there were undetermined coefficients, PR#280. o problem with function()... inside saved functions because arg count increased to accommodate stored source. (Incorrect number of arguments to "lambda") o improvemen...

...d to down columns. Now par(mfrow/mfcol = c(nr, nc)) followed by par(mfg=c(i,j)) is the preferred style. o plot( <aov.object> ) failed after 2nd plot, PR#279 o Several patches for systems with long != int. o approx failed if inputs contained NAs. o weighted lm fits with 1 diml x failed if any weight was zero. o lm.influence failed when there were undetermined coefficients, PR#280. o problem with function()... inside saved functions because arg count increased to accommodate stored source. (Incorrect number of arguments to "lambda") o improvemen...

R-1.7.0 built on NetBSD 1.6, but the validation test suite failed: Machinetype: Intel Pentium III (600 MHz); NetBSD 1.6 (GENERIC) Remote gcc version: gcc (GCC) 3.2.2 Remote g++ version: g++ (GCC) 3.2.2 Configure environment: CC=gcc CXX=g++ LDFLAGS=-Wl,-rpath,/usr/local/lib make[5]: Entering directory `/local/build/R-1.7.0/src/library' >>> Building/Updating