similar to: Problem with multinom ?

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

Hi Marc That's a workaround I can use. Thanks. I'm a newbie regarding netCDF data. Is there any information I'm losing when switching between the packages? Raphael Von: Marc Girondot [mailto:marc.girondot at u-psud.fr] Gesendet: Mittwoch, 2. August 2017 15:13 An: Felber Raphael Agroscope <raphael.felber at agroscope.admin.ch> Betreff: Re: AW: [R] Remove attribute from netcdf4

Dear Marc Thanks for your remark. I don't want to use both packages. I mentioned the package RNetCDF to show that there is a similar function I' d like to use. Raphael Von: Marc Girondot [mailto:marc.girondot at u-psud.fr] Gesendet: Mittwoch, 2. August 2017 14:51 An: Felber Raphael Agroscope <raphael.felber at agroscope.admin.ch>; r-help at r-project.org Betreff: Re: [R] Remove

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

Dear Pavel Michna [maintainer of the package] (and copy to R-help Mailing List) When I try to load the package RNetCDF (after install from CRAN) I get this message: > library("RNetCDF", lib.loc="/Library/Frameworks/R.framework/Versions/3.0/Resources/library") Error : .onLoad a ?chou? dans loadNamespace() pour 'RNetCDF', d?tails : appel : NULL erreur : I/O error

I just try the new R version 3.0.0. All is good until now. However, I had the surprise to see that all packages are not available. Until now I found one missing for OSX among these I used: maptools Is there a list of packages that were available in CRAN for 2.15 version and no more available ? And thanks a lot for all the programmers and maintainers for this great software. Sincerely, Marc

I have a problem with the RNetCDF package in MacOSX 10.8.3, R3.0.0. If you have a solution, it would be great ! Thanks a lot. Marc Girondot > install.packages("RNetCDF") essai de l'URL 'http://cran.at.r-project.org/bin/macosx/contrib/3.0/RNetCDF_1.6.1-2.tgz' Content type 'application/x-gzip' length 2071758 bytes (2.0 Mb) URL ouverte

Update taking into account the answer of Prof Ripley. ___________________________________________________ Below is a reproducing example that produces the crash. If I do exactly the same in R 3.1.2, it works perfectly. I have the lastest libcurl library (curl @7.40.0) installed So, perhaps indeed that the problem originates from RCurl, but the problem appears only in R 3.2. But before to

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

Le 02/08/2017 ? 12:03, raphael.felber at agroscope.admin.ch a ?crit : > 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

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

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

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 R-users, I am a new user of R and I am trying to build a discrete choice model (with more than two alternatives A, B, C and D) using logistic regression. I have data that describes the observed choice probabilities and some background information. An example below describes the data: Sex Age pr(A) pr(B) pr(C) pr(D) ... 1 11 0.5 0.5 0 0 1 40 1 0 0 0 0 34 0 0 0 1 0 64 0.1 0.5 0.2 0.2 ...

Dear list users, I have trouble with covariance analysis. I measured nitrate concentrations in the soil (NO3) and the percentage of legumes (LEG, continuous), affected by 2 different CO2 concentrations (CO2, discrete). I suspect that CO2 has an effect on LEG and NO3, but also that LEG has an effect on NO3, so this is the formula I wrote to test this: NO3 ~ CO2 + LEG + CO2:LEG Will LEG be

Hi, do you know, a method to convert an decimal value (integer) to the corresponding hexadecimal value ? thinks for help. Romain -- Lorrilli?re Romain UMR 8079 Laboratoire Ecologie, Syst?matique et Evolution B?t. 362 Universit? Paris-Sud 91405 Orsay cedex France tel : 01 69 15 56 85 fax : 01 69 15 56 96 mobile : 06 81 70 90 70 email : romain.lorrilliere at ese.u-psud.fr

2 years ago at one Meetup in Paris, Marc Girondot, professor at University Paris-Saclay, presented his cluster built with 12 Odroid (equivalent to Rapsberry Pi). The stack was almost the same size than your picture as there was no fan and a narrower distance between each PCB card. There were many more cables. For a total cost of less than 1000 ? and an electrical consumption reduced by a

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),

Dear R users, The nnet package includes the multinom method for the confint function. The R Help file (?confint) for the generic function in the stats package and the help files for the glm and nls methods in the MASS package indicate that one can use the "parm" argument as "a specification of which parameters are to be given confidence intervals, either a vector of numbers or

I am trying to create a LaTex table based on a multinom (nnet) object using the stargazer command. I have created a small data frame to demonstration the problem: data <- data.frame(age=1:21, hight=20:40, ed=c(1,2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,2,3)) data$ed <- as.factor(data$ed) I then make a multinomial model using the command multinom from the nnet package: model <- multinom(ed ~