similar to: Efficient multinom probs

Displaying 20 results from an estimated 793 matches similar to: "Efficient multinom probs"

2008 Apr 07
tcltk issue remains
Dear R-help, I''m trying to load the fGarch package and keep running into problems with tcltk: After succesfully instaling fGarch (and dependencies) I get: >library(fGarch) Loading required package: fBasics Loading required package: fImport Loading required package: fSeries Loading required package: robustbase Loading required package: fCalendar Loading required package: MASS
2008 Jun 03
splitting data frame based on a criteria
Hi, I have a data frame that I want to split into two based on the values of a variable in it. The variable Y has numeric values ranging between 0 through 70. I want to plot the frequencies of another variable X in two different cases: - When Y = 0 and - When Y > 0 How does one go about doing this? In general, I want to do several analyses with this data frame that are a variation of the
2006 Feb 28
any more direct-search optimization method in R
Hello list, I am dealing with a noisy function (gradient,hessian not available) with simple boundary constraints (x_i>0). I''ve tried constrOptim() using nelder mead to minimize it but it is way too slow and the returned results are not satisfying. simulated annealing is so hard to tune and it always crashes R program in my case. I wonder if there are any packages or functions can
2011 Mar 07
Teaching R: To quote, or not to quote?
Hi All, When I teach an intro workshop on R, I''ve been minimizing "quote confusion" by always using quotes around package names in function calls. For example: install.packages("Hmisc") update.packages("Hmisc") library("Hmisc") citation("Hmisc") search() # displays package names in quotes detach("packages:Hmisc") # just as
2012 Aug 27
Can anybody, please, explain me how many parameter are estimated using randomLCA? For examples, model "dentistry.lca2random" estimate 1 scale (or variance, b_j) parameter and 2 position parameters (a_cj)? Doesn''t it? Do I need at least 4 diagnostic tests for such a model? What happens if I specify options blocksize and byclass? How many diagnostic tests (or rater) I
2004 Oct 18
meta.summaries and se''s of effect sizes
Hi All, I would like to use meta.summaries from package rmeta to do a meta-analysis. I have available effect sizes as r''s (which could be easily transformed to effect sizes in terms of d''s). My problem is that I''m not sure what the se''s of these r''s should be ... The r-values are themselves computed from F-tests and t-tests for various studies. Are
2006 Mar 21
Hessian from optim()
Hello! Looking on how people use optim to get MLE I also noticed that one can use returned Hessian to get corresponding standard errors i.e. something like result <- optim(<< snip >>, hessian=T) result$par # point estimates vc <- solve(result$hessian) # var-cov matrix se <- sqrt(diag(vc)) # standard errors What is actually Hessian representing here?
2010 Nov 18
Hello All, I''m trying to run a maximum likelihood analysis using dmultinomial (i''m avoiding dmultinom as I''d like to run it with vectors for the ML stuff). However, I''m having a hard time getting even the simplest example running. Any help would be greatly appreciated. > library(mc2d) > dmultinomial(x=c(0,0,1),prob=c(1,1,1),size=1,log=TRUE) Error in if
2005 Sep 26
hidden markov models
Dear R community, I am looking for an R package or other software to study hidden Markov models. I need to be able to incorporate multivariate emissions and covariates for the transition probabilities. The msm package seems almost perfect for my purpose, but I do not think it allows multivariate emissions. I will be grateful for your suggestions. All the best, -- Emilio A. Laca One
2005 May 06
bivariate normal cdf
-- R Help List -- I am looking for a bivariate normal cdf routine in R. I have some fortran routines for this, which appear to be based on 15-point quadrature. Any guidance/suggestions on making these in loadable R-functions would be appreciated. Thanks, Dan =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Daniel A. Powers, Ph.D. Department of Sociology University of Texas at
2005 Oct 04
Animation of Mandelbrot Set
Hi, I was playing with Mandelbrot sets and come up with the following code, I thought I would share: library(fields) # for tim.colors library(caTools) # for write.gif m = 400 # grid size C = complex( real=rep(seq(-1.8,0.6, length.out=m), each=m ), imag=rep(seq(-1.2,1.2, length.out=m), m ) ) C = matrix(C,m,m) Z = 0 X = array(0, c(m,m,20)) for (k in 1:20) { Z =
2006 Mar 01
library file for R''s nmath routines
Hi, I am wondering where the library file for R''s nmath routines are? Doing a search on libR gave me the following: /usr/lib/ /usr/lib/ /usr/lib/R/lib/ /usr/lib/R/lib/ /usr/lib/ /usr/lib/ None of these have the functions in
2005 Jul 03
code for model-averaging by Akaike weights
Dear all, does anyone have r code to perform model-averaging of regression parameters by Akaike weights, and/or to do all-possible-subsets lm modelling that reports parameter estimates, AICc and number of parameters for each model? I have been looking for these in the archive but found none. (I am aware that many of you would warn me against these methods advocated by Burnham and
2008 May 09
K-Means Clustering
Hello, I am hoping you can help me with a question concerning kmeans clustering in R. I am working with the following data-set (abbreviated): BMW Ford Infiniti Jeep Lexus Chrysler Mercedes Saab Porsche Volvo [1,] 6 8 2 8 4 5 4 4 7 7 [2,] 8 7 4 6 4 1 6 7 8 5 [3,] 8 2 4
2017 Jun 29
package to fit mixtures of student-t distributions
Hello! I am new to R (before used python exclusively and would actually call the R solution for this issue inside a python notebook, hope that doesn?t disqualify me right of the batch). Right now I am looking for a piece of software to fit a 1D data sample to a mixture of t-distributions. I searched quite a while already and it seems to be that this is a somehwat obscure endeavor as most
2004 Oct 21
Cluster Analysis: Density-Based Method
Hi people, Does anybody know some Density-Based Method for clustering implemented in R? Thanks, Fernando Prass _______________________________________________________
2012 Apr 20
Dear helpers, is there any possible that transition (in depmixS4) is in scale of two variable, e.g transition=~scale(x1,x2)? If it can be, how transition of two variable (covariate time) can be worked in depmixS4-hidden markov model for time series. Many thanks, nglthu -- View this message in context: Sent from
2010 Aug 24
Rotate x-axis label on log scale
Hi I''d appreciate some help with plotting odds ratios. I want to rotate the labels on the x-axis by 45 degrees. The usual way of doing this, using text - e.g. text(1, par(''usr'')[3]-2.25..) - gives no result when the y-axis is a log scale. I guess this is because, as the par help says, for a logarithmic y-axis: y-limits will be 10 ^ par("usr")[3:4]
2005 Mar 23
non-derivative based optimization and standard errors.
Hi AlL, I ahve this problem that my objective function is discontinous in the paramaters and I need to use methods such as nelder-mead to get around this. My question is: How do i compute standard errors to a problem that does not have a gradient? Any literature on this is greatly appreciated. Jean,
2011 Mar 22
Popularity of R, SAS, SPSS, Stata, Statistica, S-PLUS updated
Greetings, I''ve just put out the latest version of "The Popularity of Data Analysis Software" at This update includes complete data for 2010, the addition of number of blogs for each software, more coverage of Statistica, and, where possible, measures regarding the implementations of the SAS Language: Carolina and the World Programming