similar to: negative log likelihood

Dear All I have three data sets, X1, X2 and Y. X1 is data, X2 and Y were generated in (different) R programs. All three vectors have one column of 60 data points. I am using the code lm(Y~X1)$coef and lm(Y~X2)$coef. The first returns two values, an intercept and a slope, but the second returns 60 values. I suspect there is something in the "type" of X2 such that it forces the regression

Hello, I use R a lot, one thing bugs me is that when I try the following > x<- -8 > x^(1/3) [1] NaN However, it is fine with -8^(1/3). Priority goes to the power. Can you help me out for this? Thanks. Best, Zhiyuan J. ZHENG Ph.D. Candidate Economic Department Virginia Polytechnic Institute and State University Phone: 540-231-5120 , Blacksburg, VA, 24060

Dear R experts, I had a question which may not be directly relevant to R but I will be grateful if you can give me some advices. I ran a two-level multilevel model for data with repeated measurements over time, i.e. level-1 the repeated measures and level-2 subjects. I could not get convergence using lme(), so I tried MLwiN, which eventually showed the level-2 variances (random effects for

Dear forum members, How can I force a negative slope in a linear regression even though the slope might be positive? I will need it for the purpose of determining the trend due reasons other than biological because the biological (genetic) trend is not positive for these data. Thanks. Julia Example of the data: [1] 1.254 1.235 1.261 0.952 1.202 1.152 0.801 0.424 0.330 0.251 0.229

I am running my code in a loop and it does not work but when I run it outside the loop I get the values I want. n <- 1000; # Sample size fitglm <- function(sigma,tau){ x <- rnorm(n,0,sigma) intercept <- 0 beta <- 0 ystar <- intercept+beta*x z <- rbinom(n,1,plogis(ystar)) xerr <- x + rnorm(n,0,tau) model<-glm(z ~ xerr, family=binomial(logit))

I am new to R and I am trying to do a monte carlo simulation where I generate data and interject error then test various cut points; however, my output was garbage (at x equal zero, I did not get .50) I am basically testing the performance of classifiers. Here is the code: n <- 1000; # Sample size fitglm <- function(sigma,tau){ x <- rnorm(n,0,sigma) intercept <- 0 beta

Dear all, After a thorough research, I still find myself unable to find a function that does linear regression of 2 vectors of data using the "total least squares", also called "orthogonal regression" (see : http://en.wikipedia.org/wiki/Total_least_squares) instead of the "ordinary least squares" method. Indeed, the "lm" function has a

Hi, I am trying to construct a logistic regression model from my data (104 patients and 25 events). I build a full model consisting of five predictors with the use of penalization by rms package (lrm, pentrace etc) because of events per variable issue. Then, I tried to approximate the full model by step-down technique predicting L from all of the componet variables using ordinary least squares

Hello, I'm a newcomer to R so please forgive me if this is a silly question. It's that I have a linear regression: fm <- lm (x ~ y) and I want to test whether the slope of the regression is significantly less than 1. How can I do this in R? I'm also interested in comparing the slopes of two regressions: fm1 <- lm (x ~ y) fm2 <- lm (a ~ b) and asking if the slope of fm1 is

Dear List: A few weeks ago I posted some questions regarding data simulation and received some very helpful comments, thank you. I have modified my code accordingly and have made some progress. However, I now am facing a new challenge along similar lines. I am attempting to simulate 250 datasets and then run the data through a linear model. I use rm() and gc() as I move along to clean up the

I have a very simple problem. When using glm to fit binary logistic regression model, sometimes I receive the following warning: Warning messages: 1: fitted probabilities numerically 0 or 1 occurred in: glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, 2: fitted probabilities numerically 0 or 1 occurred in: glm.fit(x = X, y = Y, weights = weights, start = start,

The Star Internet anti-virus service, powered by MessageLabs, discovered a possible virus or unauthorised code (such as a joke program or trojan) in an email sent by you. This email has now been quarantined and was not delivered. Please read this whole email carefully. It explains what has happened to your email, which suspected virus has been caught, and what to do if you need help. To

Dear R help My problem is very similar to the analysis detailed here. If we use the mayo dataset provided with the survivalROC package the estimate for Somer's Dxy is very negative -0.56. The Nagelkerke R2 is positive though 0.32. I know there is a difference between explained variation and predictive ability but I am surprised there is usch a difference given that even a non predictive model

I have recently upgraded my firmware from v6.3 to v7.4. Now when the phone is booted or rebooted, the initial screen "Initializing Vlan" takes forever to initialize before it initializes IP. Any ideas/Thoughts? (Trying not to Revert back to v6.3). -C -------------- next part -------------- An HTML attachment was scrubbed... URL:

Man, I really feel stupid, but after banging my head on a brick wall for several hours ... I need help! I've gotten the BLF to work on my shiny new GXP2010. The GXP is exten 5707, and I've got an xlite on 5608. When I make a call from an outside line, I dial SIP/5608. The little blinky light on the GXP that's monitoring 5608 goes, well, "blink blink". :) I then press

Hello. Consider the response-variable of data.frame df is constant, so analytically perfect fit of a linear model is expected. Fitting a regression line using lm result in residuals, slope and std.errors not exactly zero, which is acceptable in some way, but errorneous. But if you use summary.lm it shows inacceptable error propagation in the calculation of the t value and the corresponding

Has anyone in the UK got ANI to work on an inbound call ? Using asterisk 1.4 trunk and zaptel 1.4 trunk, with a Euro-ISDN 30 Julian ______________________________________________________________________ This email has been scanned by the MessageLabs Email Security System. For more information please visit http://www.messagelabs.com/email

Hi All, Having searched the archives, I can see there has been much discussion at various points regarding capture of caller id information from good old BT. If I understand correctly, it seems that not only do the drivers not currently support it, but my X101P possibly/probably can't do it anyway due to hardware? So, that leaves me with the modem route, which seems more and more unlikely,

Dear list, I'm trying to fit a glm model using family=poisson(link = "identity"). The problem is that the glm function fits a model with a negative intercept, which sounds like a nonsense to me, being the response a Poisson variable. >From a previous discussion on this list I've understood that the glm function uses IRLS for the fitting without any constraint so it is

Have asterisk@home 1.2.1 The server is on an internal network eg 10.10.10.10 It is NAT'd 1:1 via Checkpoint firewall to external public IP eg 50.50.50.50 The remote IAX2 phone (ATCOM320) is configured to call 50.50.50.50 on extension 1055. Outbound calls to 1055 work perfectly. Inbound calls from 1055 get picked up as if it were an external call (see below) and goes straight to the ring