similar to: error in fitdistr

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

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,

Hi I've written the following program: #include <R.h> #include <Rdefines.h> #include "Projector.h" SEXP Projector2(SEXP L, SEXP G, SEXP W, SEXP xymod, SEXP modif){ int nprot=0; PROTECT(L=AS_NUMERIC(L));nprot++; PROTECT(G=AS_NUMERIC(G));nprot++; PROTECT(W=AS_NUMERIC(W));nprot++; PROTECT(xymod=AS_INTEGER(xymod));nprot++; PROTECT(modif=AS_NUMERIC(modif));nprot++;

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 approp elements then reshaping it to an array, but I wonder if there isn't a

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

Full_Name: Christina Merz Version: R version 2.5.0 (2007-04-23) OS: mingw32 Submission from: (NULL) (213.70.209.132) R> version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status

Dear All I can not download R for Linux. I do not know which file I should install? Best Regards, Soheila

Dear All Is there any function for "median test" in R? Best Regards, Soheila [[alternative HTML version deleted]]

Dear All How can I transform 1 column to 2 columns in R? 211217_s_at GO:0005249 211217_s_at GO:0005251 211217_s_at GO:0005515 211217_s_at GO:0015271 211217_s_at GO:0030955 211217_s_at GO:0005249 211217_s_at GO:0005251 211217_s_at GO:0005515 211217_s_at GO:0015271 211217_s_at GO:0030955 Best Wishes, Soheila [[alternative HTML version deleted]]

Dear All I want to arrange our output in 3 columns in R. Example: N 04320 03010 00514 00604 00563 03430 03022 V 22.000000 0.000000e+00 15.0000000 15.0000 13.0000 1.300000e+01 51.0000000 p 0.640625 2.910383e-11 0.4257812 0.4375 0.6875 3.356934e-04 0.2434998 N 03060 04140 03050 00562 00565 V 20.0000000 62.00000000 273.0000000

var1 var2 var3 var4 var5 var6 var7 var8 var9 var10 gold 2 3 1 2 4 0 1 4 4 3 2 2 4 2 4 3 4 2 4 4 4 2 3 3 0 0 4 1 0 2 4 4 2 1 4 0 3 2 0 0 2 4 4 2 3 4 0 2 2 0 0 0 3 4 2 2 2 3 2 2 0 0 0 2 4 2 2 4 1 1 2 0 0 3 3 3 2 3 4 1 4 0 0 0 0 3 4 2 3 1 0 2 2 1 0 2 3 3 2 0 3 1 1 1 1 2 1 2 3 2 1

Dear All I have 54000 plots in R, How can I observe them? If I‌ have to save them one-by-one? Soheila [[alternative HTML version deleted]]

Dear All I want to chose just spacial columns in R. (read table) data1<- read.table("/home/Documents/data.txt",header=F,sep = "\t", as.is =F) data.2<-data1[sub.data[,1],] The dimension of data.2 is correct but my data are transformed NA head(data.2) V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18 V19 V20 V21 V22 V23 V24 NA <NA> NA NA NA

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)

Dear sirs, I have tried to find a command to fit(plot) a regression line but I haven't found any command on the html.help. My question is if I am using a regression for instance reg<-lm(y~x) and I would like to make a regression line, to see how well my data is fit. Can you please help me with this?? Yours Sincerely Soheila Zad

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,

Dear ALL I have this function in R: func_LN <- function(data){ med_ge <- matrix(c(rep(NA,nrow(data)*ncol(data))), nrow = nrow(data), ncol=ncol(data), byrow=TRUE) T <- matrix(c(rep(NA,length(n)*ncol(data))), nrow = length(n), ncol=ncol(data), byrow=TRUE) Tdiff<- matrix(c(rep(NA,length(n)*ncol(data))), nrow = length(n), ncol=ncol(data), byrow=TRUE) T1<- c(rep(NA,ncol(data)))

>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 =

We were looking at conducting a pilot program with Dell to virtualize our lab environments. We are specifically looking at Dell/Citrix provisioning servers. The OS in this setup is streamed from the server on every boot. I was curious if the R Project had any licensing stipulations in this type of setup. This is in the beginning phases and we wanted to answer some questions about licensing.

Prof Brian Ripley wrote: > > Even better, use the function fitdistr in package MASS... > I had a look in my MASS library (from the package VR_6.2-6) and couldn't find this function. Is there a newer version available? Thanks for any help. Jim This email message and any accompanying attachments may contain confidential information. If you are not the intended recipient, do not