similar to: optima in unimode

I'm attempting to do model selection with AIC, using a glm and a lognormal distribution, but: fit1<-glm(BA~Year,data=pdat.sp1.65.04, family=gaussian(link="log")) ## gives the same result as either of the following: fit1<-glm(BA~Year,data=pdat.sp1.65.04, family=gaussian) fit1<-lm(BA~Year,data=pdat.sp1.65.04) fit1 #Coefficients: #(Intercept) Year2004 # -1.6341

Dear R users, the first 1.0 version of the rPref package is now on CRAN. rPref allows to select the Pareto-optimal tuples from a data set, also called Skylines in the database community. For example, optimal tuples from mtcars according to "high(mpg) * high(hp)" (where "*" is the Pareto operator) are those cars, for which no other dominating car exists. There, a car

Hello, I'm trying to make a lattice plot (using xyplot()). I have included a "layout=c(3,4)" statement, giving me 12 plots per page and an "as.table=TRUE" statement, directing the way the plots are laid out. I have 18 plots altogether and so 6 of them end up on the second page. Everything looks fine for the first page, but the x-axis labels (e.g. 1993, 1994...) are all

About two years ago there was a thread about this which suggested that at that time nobody had these coefficients ready to go. (a) has anyone in the meanwhile programmed them? (b) I think I can see how to do the tetrachoric one with mvtnorm on similar lines to an example on the help page so will try that if nobody else already has (c) looking at the polychoric one makes me realise yet again

Per the Posting Guide, why didn't you post the reproducible R code example? On April 24, 2018 8:22:15 PM PDT, Donald Macnaughton <donmac at matstat.com> wrote: >I'm drawing a paneled histogram using the lattice package. I've >succeeded in >adding minor tick marks to the vertical axis, but I can't get the >desired >number of minor tick marks between the major

I'm drawing a paneled histogram using the lattice package. I've succeeded in adding minor tick marks to the vertical axis, but I can't get the desired number of minor tick marks between the major tick marks. I've attached a self-contained program to illustrate the problem. Thanks for your help, Don Macnaughton Here's my sessionInfo: R version 3.4.3 (2017-11-30) Platform:

Hello R Help, I am trying solve an MLE convergence problem: I would like to estimate four parameters, p1, p2, mu1, mu2, which relate to the probabilities, P1, P2, P3, of a multinomial (trinomial) distribution. I am using the mle2() function and feeding it a time series dataset composed of four columns: time point, number of successes in category 1, number of successes in category 2, and

Thanks Jeff, I attached a file with the program to my earlier email because the posting guide seemed to imply that non-binary attachments would work. But I see that the file was stripped off. I installed the program file on a web site, but when I downloaded it, the line breaks were stripped out. So I've included the program below: ------------------------------------------------------- #

Dear List, I have a richness data distributing across 20 N to 20 S latitude. (120 E-140 E longitude). I would like to draw the richness in the north hemisphere and a regression line in the plot (x-axis: latitude, y-axis: richness in the north hemisphere). The above demand is done using plot. Then, south hemisphere richness and regression are required to be generated using the same y-axis above

hi, i am writing my master thesis and i am dealing with 146474 observations (panel data), i have just learned the R so i am a beginner!! i am trying to use the "plm" package and i have a duplication problem; i have written the following commands to read my data and create my model >dsn<-plm.data(ds, c("stno", "year")) ds=name of my data, stno=individual

Hi, I have a share in my smb.conf file that I want to share to members of a supplemental group, but I keep getting a network error. here is the share: [mp3s] comment = MPEG3 Encoded Audio path = /usr/share/mp3 public = no browsable = no valid users = +mp3users writable = yes create mask = 0644 I have a group called mp3users that inside my /etc/group file looks like this:

Hi, everyone, Night. I have three questions about multiple linear regression in R. Q1: y=rnorm(10,mean=5) x1=rnorm(10,mean=2) x2=rnorm(10) lin=lm(y~x1+x2) summary(lin) ## In the summary, 'Residual standard error: 1.017 on 7 degrees of freedom', 1.017 is the estimate of the constance variance? Q2: beta0=lin$coefficients[1] beta1=lin$coefficients[2] beta2=lin$coefficients[3]

Hi all, I try to interpolate a data set in the form: time Erg 0.000000 48.650000 1.500000 56.080000 3.000000 38.330000 4.500000 49.650000 6.000000 61.390000 7.500000 51.250000 9.000000 50.450000 10.500000 55.110000 12.000000 61.120000 18.000000 61.260000 24.000000 62.670000 36.000000 63.670000 48.000000 74.880000 I want to get smoothed splines by using the class gam The first way I tried , was

Dear All, A method for the estimation is univariate unimodal densities (with unknown mode) is described in "Statistical Inference under Order Restrictions" by Barlow et al.. Would anyone know whether there is an R-implementation (preferably with reference) for the estimation of univariate discrete unimodal densities (with unknown mode)? Thanks in advance for your help. Kind

Buenas a todos, Lo primero agradecer todas las respuesta sque tuve en el tema de Bootstrap dataframe, que por estar de baja no he podido agradecer. De aquel tema salió una sugerencia que me parece muy interesante y que a dia de hoy no soy capaz de hacer de una manera optima. Lo que quiero hacer es coger un dia al azar de todo el periodo, y a partuir de ese dia, coger por ejemplo los 20 dias

I had Wine-20050310 installed from FreeBSD ports. This version worked fine for me. But after upgrading to Wine-20050524 some applications don't work properly. E.g. Intranet Chat: http://www.softpedia.com/get/Internet/Chat/LAN-Messaging-Clients/Intranet-Chat.shtml. Running this applications I see a following message: fixme:system:SystemParametersInfoW Unimplemented action: 8193

Dear All, Does anybody know if there is a functionality in R to break histograms that show a clear bi-modal (or multi-modal) distribution into a series of unimodal histograms that added up result in the original histogram? I was thinking of using QQ-plots (for which tools are available in R), and then observing the number of times the observed quantiles cross the 1:1 line, but this only gives an

Dear R helpers, I wish to use the "sum" operator for each row of a data frame. However, it appears that the operator acts on the entire data frame, over all columns. What is the best way to obtain row- wise operation? The following code shows my attempts so far, and their problems:- test1=array(rbinom(120,1,0.5),c(20,3)) test1[,3]=NA sum(test1[,1:2]) test1[,3][sum(test1[,1:2])>=2]=1

Hi Helpers, Does anyone know how to read jpeg, bmp or png files into R? I have some photos of brain scans and I want to quantify some aspects of their size. I might be able to do this with the 'locator' function, if I could figure out how to read the files in and make them into an image that I can display on the R Windows device. I have experimented with readBin, using a simple

Jonathan, Try making a list instead of an array. See ?list. Also, did you look into random forests? I'm not sure what you want to do, but there might be methods there to do some of the work for you. Sean On 3/19/04 1:12 PM, "Jonathan Williams" <jonathan.williams at pharmacology.oxford.ac.uk> wrote: > I would like to collect the trees grown by rpart fits in an array,