similar to: Logistic regression with more than two choices

I want to fit a multinomial model with logit link. For example let this matrix to be analyzed: male female aborted factor 10 12 1 1.2 14 14 4 1.3 15 12 3 1.4 (this is an example, not the true data which are far more complex...) I suppose the correct function to analyze these data is polr from MASS library. The data have been

Dear members, I have a series of values in a vector and some value are missing and replaced with NA. For example: a <- c(27, 25, NA, NA, 24, 26, 27, NA, 26) I would like to replace the NAs with the value taken from the previous value that is non-NA. The output would be in this case: 27 25 25 25 24 26 27 27 26 Now I do that with a for loop, but I try to eliminate all the loops to gain in

Dear all For a model I need to combine several netCDF files into one (which works fine). For better overview I'd like to delete/remove some of the attributes. Is there a simple way doing this? I'm using the package netcdf4, which creates an object of class(nc) = "ncdf4". It seems that for earlier versions of netcdf objects, there was the function att.delete.nc{RNetCDF}. But

I define priors in jags within r using a gamma distribution. I would like to control the shape but I have problems. Any help will be usefull. From help of dgamma ___________________ The Gamma distribution with parameters shape = a and scale = s has density f(x)= 1/(s^a Gamma(a)) x^(a-1) e^-(x/s) and rate=1/scale From jags user manual ____________________ dgamma(r, mu) has a density of

Dear All: I am interested in multinomial logit models (function multinon, library nnet) but I'm having troubles in choose whether to define the predictors as factors or not. I had posted earlier this example (thanks for the reply ronggui): worms<- data.frame(year= rep(2000:2004, c(3,3,3,3,3)),age=rep(1:3,5),

Hi all, I know there are some howto's about Powerware 5115 and Ubuntu, but I can't still get that thing working. Can someone give some light... 1. Powerware 5115, serial port (ttyS0) 2. Ubuntu 6.06.1 as LTSP-server As reference: http://gentoo-wiki.com/HOWTO_NUT http://keystoneit.wordpress.com/2006/09/25/network-ups-tools-nut-on-ubuntu/ --- a. I try to start driver as root. root at

Hi all, I know there are some howto's about Powerware 5115 and Ubuntu, but I can't still get that thing working. Can someone give some light... 1. Powerware 5115, serial port (ttyS0) 2. Ubuntu 6.06.1 as LTSP-server As reference: http://gentoo-wiki.com/HOWTO_NUT http://keystoneit.wordpress.com/2006/09/25/network-ups-tools-nut-on-ubuntu/ --- a. I try to start driver as root. root at

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20120318/d68edf35/attachment.pl>

Hello, I don't know if the problem originates from R-devel 3.2 or Rcurl itself. I post this message to the R-devel list and to the author of RCurl (duncan at r-project.org). > library("RCurl") Le chargement a n?cessit? le package : bitops > print(sessionInfo()) R Under development (unstable) (2015-02-03 r67717) Platform: x86_64-apple-darwin13.4.0 (64-bit) Running under:

I use R 2.15.2 on MacosX 10.8.2 I have a problem checking for a package that I want to submit in CRAN. When I check like that, no problem or warning are detected. system("R CMD check '[path]/phenology'") I create the .tar.gz package with that. The file phenology_3.43.tar.gz is generated. I can install the package and it works. system("R CMD build

Hi, I found that using different contrasts (e.g. contr.helmert vs. contr.treatment) will generate different fitted probabilities from multinomial logistic regression using multinom(); while the fitted probabilities from binary logistic regression seem to be the same. Why is that? and for multinomial logisitc regression, what contrast should be used? I guess it's helmert? here is an example

Dear R users, Does anyone know if it is possible to use multinom to do a polychotomous fit using one categorical and one numeric variable as response. The doc. for multinom states that for formula , response can be K>2 classes. Is this 2 and more, or as I have understood it only greater than 2. I have tried fitting my data, but have only encountered error messages. On another note, Is it

Hi, I'm new to R. I have some daily soil moisture data for the year 1979 in netcdf format such as these sm19790101.1.nc sm19790102.1.nc . . . sm19791231.1.nc I need to average a variable called "sm" to monthly resolution. I've done these days = formatC(1:31, width=2, flag="0") ncfiles = lapply(days, function(d){ filename = paste("sm197901", d,

When ploting a timeseries, the months are shown with abbreviation in the current language of the system. For example, x <- seq(from=as.Date("2000-04-01"), length.out=100, by=1) y <- rnorm(length(x), 5,2) plot(x, y) Show for me "avi mai jui jul" as I use R with French language localization. I see in the Windows FAQ how to completely change the language of R, but I

If one estimates a model using multinom, is it possible to perform the omnibus LR test ( the analogue to omnibus F in linear models ) using the output from multinom ? The residual deviance is there but I was hoping I could somehow pull out the deviance based on just using an intercept ? Sample code is below from the CAR book but I wasn't sure how to do it based on that example. Thanks

Hi, all I have a medium sized data, 6 years. Each observation is a case with a date variable, such as '2004-08-02'. Some of the months didn't occur a case. I want to plot the 6 years data by month, and the Y_axis is the freqency of cases for each month, meaning 12*6=72 bars or points in the figure. I though of a method, 1st, using the months function, then ploting. But I need to

Hi, I try to make my first package? The HelloWorld.R file is: #### HelloWorld.R #### #' showHello est une fonction R permettant d'afficher le message #' "Hello World!" sur la console. #' @title la fonction showHello() showHello <-function(){ cat("Hello World!\n") } I use the following procedure to get the tar: # set the working directory where the file is

I would like to create a function to convert data based on a class using as.numeric Similarly, using a plot, I do: plot.essai <- function(x, ...) {return(x*2)} d <- 10 class(d) <- "essai" plot(d) It works: [1] 20 attr(,"class") [1] "essai" Now same with as.numeric: as.numeric.essai <- function(x, ...) {return(x*2)} as.numeric(d) It does not work:

Dear all, I was wondering how to read DNA sequences in R, is there a specific function and/or a specific package for that? Thank you very much in advance, Gian ** [[alternative HTML version deleted]]

Hello all. I am using method "bag()" to develop a bagging with my dataset. When I do bag( dataP, dataP$Score, B=10) I obtain this error: error en bag.default(Score ~ ., data = dataP, B = 10) : entrada en evaluacion: recursivo por defecto o problemas anteriores? ( in english, evaluating input: default recursive or previous problems?) Any idea? Thanks! [[alternative HTML