similar to: Question in R

Hi, I tried to calculate the MLE estimation of a parameter of a gamma distribution of a data set, which is called "amount".. I did it with the command: fitdistr(amount,"gamma") and got the error message: Error in optim(x = c(363.36, 551.15, 603.91, 170.99, 225.21, 218.01, 1816.82, : non-finite finite-difference value [2] Warning message: NaNs wurden erzeugt

https://bugs.freedesktop.org/show_bug.cgi?id=92991 Bug ID: 92991 Summary: Nouveau module loading regression in 4.4-rc1 on optimus system Product: xorg Version: unspecified Hardware: Other OS: All Status: NEW Severity: normal Priority: medium Component: Driver/nouveau

Hi all, I'm trying to apply loess regression to my data and then use the fitted model to get the *standardized/studentized residuals. I understood that for linear regression (lm) there are functions to do that:* * * fit1 = lm(y~x) stdres.fit1 = rstandard(fit1) studres.fit1 = rstudent(fit1) I was wondering if there is an equally simple way to get the standardized/studentized residuals for a

Hello, I am using the Ljung Box test in R to compute if the resiudals of my fitted model is random or not. I am not sure though what the results mean, I have looked at various sources on the internet and have come up with contrasting explanations (mainly because these info deal with different program languages, like SAS, SPSS, etc). I know that my residuals should appropriate white noise( is

Dear list, I would like to extract from a matrix all those rows, that are unique. By unique, I don't mean the unique that is accomplished by the function unique(), though... Consider the following example: > h [,1] [,2] [1,] 4 4 [2,] 1 4 [3,] 4 1 Now unique(h) returns exactly the same - because 1 4 and 4 1 is not the same for that function. What I would like to

"The R Inferno" is now on the Burns Statistics website at http://www.burns-stat.com/pages/Tutor/R_inferno.pdf Abstract: If you are using R and you think you're in hell, this is a map for you. Also, I've expanded the outline concerning R on the Burns Statistics 'Links' page. Suggestions (off-list) for additional items are encouraged. Patrick Burns patrick at

There's a new page on the Burns Statistics website http://www.burns-stat.com/pages/Tutor/spreadsheet_addiction.html that looks at spreadsheets from a quality assurance perspective. It presents R as a suitable alternative to spreadsheets. Also there are several specific problems with Excel that are highlighted, including the status of statistical functionality in Excel. Patrick Burns Burns

'assign' does not give a warning if 'x' has length greater than 1 -- it just uses the first element: assign(c('a1', 'a2'), 1:2) One way of thinking about this is that people using 'assign' get what they deserve. The other is that it is used seldom enough that adding a warning isn't going to slow things down appreciably. Patrick Burns patrick at

Hello, I don't understand the description / help-text for the numericDeriv() function. Why is there a new environment used? And what is meant with an environment here? Is it similar or the same as a local workspace, like an environment in functional languages? And why is it needed here? numericDeriv could just calculate the difference bewtween two values and divide this difference by the

Is there a wildcard value for vector values in r? For instance: > M <- *wildcard > (M==1) TRUE >(M=="peanut butter") TRUE >is.na(M) FALSE thanks, Francis -- Francis Smart (406) 223-8108 cell

Hi R-users Can the users let me know how to do MGARCH estimate (Bivariate GARCH) and volatility forecast for 2 variables in R. thanks and regards snvk

Hello, I have a data frame that has three fields. Resp# ActCode ProdUsed 100 3 2 100 3 2 100 4 3 100 4 3 101 3 6 102 2 1 102 3 1 103 5 1 103 5 1 103

You may recall that there was a discussion of a technical report from the statistical consulting group at UCLA. I have a draft of a comment on that report, which you can get from http://www.burns-stat.com/pages/Flotsam/uclaRcomment_draft1.pdf I'm interested in comments: corrections, additions, deletions. Patrick Burns patrick at burns-stat.com +44 (0)20 8525 0696 http://www.burns-stat.com

I have a function declared thus. FirstEigenvectorBoundary.Training <- function(InputFileName='C:/Samples2/PT_Amp.txt', Header=TRUE, Colour="red") Inside the function, I have the following call out<-list(x=Eigenvectors[2:(NumMetricsSelected+1),1], y=-0.8, z=NumMetricsSelected); NumMetricsSelected has the value 2 and Eigenvectors has the following form [,1]

Hi Is it possible to subset an n-dimensional array by a vector of n dimensions? E.g. assume I have > x <- array(1:24, dim=2:4) > x[1,1,2] [1] 7 > dims <- c(1,1,2) I would like a function that I can supply x and dims as parameters to, and have it return 7. Also, I would like to do something like: > x[1,1,] [1] 1 7 13 19 > dims2<- c(1,1,NA) And have a function of x and

Hi, I'm working on a new ogg file chopper (porting 'hogg chop' to C). These work pagewise, ie. they just copy the pages that are needed from the input to the output bitstream. I'm modifying the last page in the output bitstream so that it has the eos flag set. The last page is always one that has a granulepos set (that of the last packet finishing on that page). But let's say

Hi, I have a binary (o/1 - coded) data set and want to test it's autocorrelation structure. Is that function implemented in R? Can I use the ACF - funtion with binary data? Thanks for your help, Daniel

Unix versions of R CMD BATCH have reported proc.time() unless the script ends in q(). E.g. if the input is 'search()' the output is > invisible(options(echo = TRUE)) > search() [1] ".GlobalEnv" "package:stats" "package:graphics" [4] "package:grDevices" "package:utils" "package:datasets" [7]

I have read some very good reviews comparing R (or Splus) to SAS. Does anyone know if there are any reviews comparing R (or Splus) to Stata? I am trying to get others to try R in my department, and I have never used Stata. Regards, -Cody Cody Hamilton, Ph.D Institute for Health Care Research and Improvement Baylor Health Care System (214) 265-3618 This e-mail, facsimile, or letter

i need to know how to estimate a linear regression whose coefficients sum to zero