similar to: 2-parameter MLE problems

Hello! I have 2 variables - predictor "pred" and response variable "DV": pred<-c(439635.053, 222925.718, 668434.755, 194242.330, 5786.321, 115537.344, 100835.368, 7133.206, 159058.286, 4079991.629, 3380078.060, 2661279.136, 2698324.478, 1245213.965, 1901815.503, 1517019.451, 1396857.736, 1034030.988, 432249.574, 342329.325, 1831335.792, 2209578.859, 1641709.754,

Hello there, I'm trying to do lots of plots in one for-loop. But somehow ggplot does not evaluate arguments as expected. Here is an example: library(lattice) library(ggplot2) pl <- list() pl2 <- list() cDat <- as.data.frame(cbind(x1=0:100,x2=0:10,x3=1:20)) for(obs in c("x1", "x2")) { pl[[obs]] <- xyplot(cDat[,obs] ~ cDat[,"x3"], main=obs)

I'm puzzled by the following problem, which appears when attempting to run an analysis on part of a dataset: If I try: csubset <- dat$Diagnosis==0 cont <- lme(fixed=cform, random = ~1|StudyName, data=dat,subset=csubset,na.action=na.omit) Then I get: Error in eval(expr, envir, enclos) : Object "csubset" not found But if I do

Dear R experts, I have a labeled data set. Each data is assigned a binary label 0 or 1. Assume that I use some clustering algorithm to group the data by clusters (using some features of the data). Now I want to know how many data are labeled as 0/1 in each cluster. For example, assume that I have 9 labeled data grouped into three clusters. The ids of the clusters are 1, 2, and 3. The dataset is

Hello all, I what to print the banner plot that is output from the mona method in the cluster package but the problem is I dont want to print all that red ink. Here is an example: data(animals) ma <- mona(animals) ma ## Plot similar to Figure 10 in Struyf et al (1996) plot(ma) I can change the bar color by using the argument col=c(0,0) - plot(ma,col=c(0,0)) - but then the variable labels also

Hi, I have a table for an 1 week exam result for many classes in school, like this: Day Class_ID Test Result 1 Monday 1 Paper Passed 2 Tuesday 1 Oral Passed 3 Friday 1 Paper Passed 4 Monday 3 Paper Passed 5 Sunday 3 Oral Passed 6 Monday 3

Hi, I want to insert a key into each panel of a trellis plot, which I can do with a custom panel function that calles draw.key. The problem arises because I want the top right hand corner of the key to start in the top right hand corner of the panel. If you run my code below, you can see that the key appears in the center of each panel. This is because the default viewport in draw.key is the

I am having a bit of a problem getting AMD to work on a new server. On my regular office server it works like a charm. I am running Asterisk 1.4.13, Zaptel 1.4.5.1 on both machines. Both servers run CentOS 5 and I am using a SIP trunk to send out calls (the same one on both servers). Here is the output of a call on my office server: -- Attempting call on Local/0445540881644 at CC2 for

We have an Asterisk server with a small outgoing call center. We use AMD and it usually works very well on Zap channels (E1 PRI). We added a couple of SIP trunks to reduce long distance costs but now AMD gets stuck when the call goes out through the SIP channels. Here is an example call using a SIP line: -- Executing [016566275538 at CC2:1] Set("Local/016566275538 at CC2-dad7,2",

Hi all, If I run LDA on the same data (2-class classification) with default(no priors specified in the lda function) vs. "prior=c(0.5, 0.5)", the results are different. The (0.5, 0.5) priors give better 1-classify-to-1 rate, and the proportional priors(default, no priors specified in the lda function) give better 0-classify-to-0 rate, for both training and testing data sets. However,

In the following code, the only difference between the two plots is the order the variables are plotted. In this case, the plot of "cdata.den" in plot #1 is different from its plot in #2. Specifically, "cdata.den" spans the x-axis from -5 to 30 in plot #1 and from 0 to 20 in plot #2. Does anyone understand why these two plots do not yield the same result? #Make density

Hello, I'm trying to develop a box plot of time series data to look at the range in the data values over the entire period of record. My data initially starts out as a list of hourly data, and then I've been using this code to make this data into the final ts array. # Read in the station list stn.list <- read.csv("/home/kbennett/fews/stnlist3", as.is=T, header=F) # Read in

Hallo all, I have the following glm model: f1 <- as.formula(paste("factor(y.fondi)~", "flgsess + segmeta2 + udm + zona.geo + ultimo.prod.", "+flg.a2 + flg.d.na2 + flg.v2 + flg.cc2", " +(flg.a1 + flg.d.na1 + flg.v1 + flg.cc1)^2", " + flg.a2:flg.d.na2 + flg.a2:flg.v2 +

Dear Bert, Thank you so much for your kind and valuable feedback. I tried finding the starting values using the approach you mentioned, then did the following to fit the nonlinear regression model: nlregmod2 <- nls(y ~ theta1 - theta2*exp(-theta3*x), start = list(theta1 = 0.37, theta2 = exp(-1.8), theta3 =

Hi all, please see https://reviews.llvm.org/D85603 and its related changes for our most recent and hopefully final attempt at putting the `convergent` attribute on a solid theoretical foundation in a way that is useful for modern GPU compiler use cases. We have clear line of sight to enabling a new control flow implementation in the AMDGPU backend which is built on this foundation. I have

Hi Hal, On Mon, Aug 17, 2020 at 2:13 AM Hal Finkel <hfinkel at anl.gov> wrote: > Thanks for sending this. What do you think that we should do with the > existing convergent attribute? My preference, which is implicitly expressed in the review, is to use `convergent` both for the new and the old thing. They are implicitly distinguished via the "convergencectrl" operand

On Mon, Aug 17, 2020 at 7:14 PM Hal Finkel <hfinkel at anl.gov> wrote: > On 8/17/20 11:51 AM, Nicolai Hähnle wrote: > > Hi Hal, > > > > On Mon, Aug 17, 2020 at 2:13 AM Hal Finkel <hfinkel at anl.gov> wrote: > >> Thanks for sending this. What do you think that we should do with the > >> existing convergent attribute? > > My preference, which

Hi all, I have a has_many relationship, Cookbook :has_many Recipes If I want to create a new recipe in my cookbook I do: cookbook.recipes.create(data) It would be nice if I could then easily get the id of the new recipe, but the create() method doesn''t seem to return anything. Any ideas how? Thanks in advance! Cameron -- Posted via http://www.ruby-forum.com/.

Hi The fix is very simple: the defaulting of this macro to 0 in gcc/config/arm/arm.h, present in 2.2 has been omitted in 2.5. The attached patch puts the defaulting clause back, the same as it was in 2.2 M -------------- next part -------------- A non-text attachment was scrubbed... Name: 0004-arm-default-MACHO_DYNAMIC_NO_PIC_P.patch Type: text/x-diff Size: 385 bytes Desc: not available

Hello >   The fix is very simple: the defaulting of this macro to 0 in > gcc/config/arm/arm.h, present in 2.2 has been omitted in 2.5. > > The attached patch puts the defaulting clause back, the same as it was in 2.2 This is known problem which was fixed after 2.5 release. -- With best regards, Anton Korobeynikov Faculty of Mathematics and Mechanics, Saint Petersburg State