Displaying 20 results from an estimated 2000 matches similar to: "Error message with glm"
2011 Dec 22
0
Finding predicted probabilities
I ran three logit models in R with the Zelig package and I'm trying to
compute the
predicted probabilities for a number of different values on the independent
variable.
My dep variable was accepted or decline and my indep variable is bid
amount, and varies.
So for a bid amount of 3, what's the expected probability of winning.
For a bid amount of 5, what's the expected probability of
2012 Aug 13
3
Using the effects package to plot logit probabilities
I'm trying to run a logit model and plot the probability curve for a number
of the important predictors. I'm trying to do this
with the Effects package.
df=data.frame(income=c(5,5,3,3,6,5),
won=c(0,0,1,1,1,0),
age=c(18,18,23,50,19,39),
home=c(0,0,1,0,0,1))
str(df)
md1 = glm(factor(won) ~ income + age + home,
data=df,
2012 Jul 19
3
Removing values from a string
So I have the following data frame and I want to know how I can remove all
"NA" values from each string, and also
remove all "|" values from the START of the string. So they should
something like "auto|insurance" or "auto|insurance|quote"
one = data.frame(keyword=c("|auto", "NA|auto|insurance|quote",
"NA|auto|insurance",
2012 Jul 05
2
Plotting the probability curve from a logit model with 10 predictors
I have a logit model with about 10 predictors and I am trying to plot the
probability curve for the model.
Y=1 = 1 / 1+e^-z where z=B0 + B1X1 + ... + BnXi
If the model had only one predictor, I know to do something like below.
mod1 = glm(factor(won) ~ as.numeric(bid), data=mydat,
family=binomial(link="logit"))
all.x <- expand.grid(won=unique(won), bid=unique(bid))
y.hat.new
2011 Nov 17
1
Error When Installing the RODBC Package
I'm running R in Ubuntu 10.10 and am trying to install the RODBC package.
However, I get the following error message:
ERROR: configuration failed for package ‘RODBC’
* removing ‘/home/amathew/R/i686-pc-linux-gnu-library/2.13/RODBC’
The downloaded packages are in
‘/tmp/RtmpekzPOQ/downloaded_packages’
Warning message:
In install.packages() :
installation of package 'RODBC' had
2012 Aug 02
1
Naive Bayes in R
I'm developing a naive bayes in R. I have the following data and am trying
to predict on returned (class).
dat = data.frame(home=c(0,1,1,0,0), gender=c("M","M","F","M","F"),
returned=c(0,0,1,1,0))
str(dat)
dat$home <- as.factor(dat$home)
dat$returned <- as.factor(dat$returned)
library(e1071)
m <- naiveBayes(returned ~ ., dat)
m
2011 Dec 21
1
Predicting a linear model for all combinations
Lets say I have a linear model and I want to find the average expented
value of the dependent variable. So let's assume that I'm studying the
price I pay for coffee.
Price = B0 + B1(weather) + B2(gender) + ...
What I'm trying to find is the predicted price for every possible
combination of values in the independent variables.
So Expected price when:
weather=1, gender=male
weather=1,
2012 Feb 09
1
Grouping together a time variable
I have the following variable, time, which is a character variable and it's
structured as follows.
> head(as.character(dat$time), 30) [1] "00:00:01" "00:00:16" "00:00:24" "00:00:25" "00:00:25" "00:00:40" "00:01:50" "00:01:54" "00:02:33" "00:02:43" "00:03:22"
[12]
2012 Aug 08
1
Calculating percentages across multiple columns
I have the following data and am trying to find the percentage of bid
values purchased for that price.
So let's say I have a bid of 5 and it's sold 2 times for $3 and $5. Since
the original bid was $5, the
percentage of times that that bid value results in a sold purchase AT that
specific bid level was
1/3 because of the three time where the bid was three, it ended up being
sold for $5
2011 Dec 15
1
Reordering a numeric variable
I'm running a linear model in R using the car package.
I have a variable education, which i have recoded and regrouped to my
wishes.
However, R seems to place each element of that variable in alphabetical
order.
When I am running the model, don't I need the model order from lowest to
highest to make an inference that
a one unit change in one variable produced a one unit change in
2012 Sep 11
1
Plotting every probability curve
I don't have a logistic regression model and am trying to generate
probability curves for all possible combinations of
the variables. My logit model has 5+ variables, and I want to draw curves
for every scenario.
See code below. When home_owner is 0 and 1, I want curves. The same goes
for all other variables categories, so that
I have permutations for all possible combinations.
I've
2012 Aug 07
2
Re-grouping data in R
I have a data frame with a column of values that I want to bucket (group)
into specific levels.
> str(dat)'data.frame': 3678 obs. of 39 variables:
$ id : int 23 76 129 156 166 180 200 214 296 344 ...
$ final_purchase_amount : Factor w/ 32 levels
"\\N","1082","1109",..: 1 1 1 1 1 1 1 1 1 1 ...
So I ran the following to
2011 Dec 16
1
Zellig Error Message
I'm trying to calculate predicted probabilities in R with Zelig and keep
getting the following error.
Can anyone help?
> x.low <- setx(mod, type=1)Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] :
incorrect number of dimensions
When I ran the model, I ran everything but the explanatory variable as a
numeric variable. Now, I'm trying everything and no
2012 Sep 26
1
Specifying a response variable in a Bayesian network
I'm trying to teach myself about Bayesian Networks and am working with the
following data and the bnlearn package.
I understand the conceptual aspects of BNs, but I'm not sure how to specify
the response variables in R when constructing
a dag plot. I've cecked ?hc and done numerous google searches without luck.
Can anyone help?
library("bnlearn")
2011 Jul 23
2
xml2-config issues
I'm trying to install the XML package on Ubuntu 10.10, and I keep getting
a warning message the XML could not be found and had non-zero exit
status. How can I fix this problem?
> install.packages()
Loading Tcl/Tk interface ... done
--- Please select a CRAN mirror for use in this session ---
Installing package(s) into ‘/home/amathew/R/i686-pc-linux-gnu-library/2.13’
(as ‘lib’ is
2011 Dec 16
0
Error constructing probabilities in Zelig
I've run an ordered logistic regression model in R with Zelig and am
looking to calculate predicted probabilities. Zelig has a series of simple
one line commands to generate the information I want on first differences
and so forth. Unfortunately, I keep getting an error when running the zelig
function and was wondering if there was a quick alternative for generating
predicted probabilities for
2011 Dec 16
0
Incorrect Number of Dimensions in Zelig with setx()
I'm running an ordered logit in R with the Zelig package and am trying to
calculate some predicted probabilities. However, I get the following error
message.
> x.low <- setx(mod, cars=1)Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] :
incorrect number of dimensions
I googled this problem and couldn't find anything, minus a question by
me on this same
2010 Mar 29
1
Question about 'logit' and 'mlogit' in Zelig
I'm running a multinomial logit in R using the Zelig packages. According to str(trade962a), my dependent variable is a factor with three levels. When I run the multinomial logit I get an error message. However, when I run 'model=logit' it works fine. any ideas on whats wrong?
## MULTINOMIAL LOGIT
anes96two <- zelig(trade962a ~ age962 + education962 + personal962 + economy962 +
2009 Mar 15
3
read.xls question
I'm an R newbie and had a question about the read.xls function. I've heard that this is often not a reliable function to use for importing data. However, I have created numerous xls files which contain information about voter turnout and macroeconomic indicators in India. I'm writing a paper on the relationship between economic growth and voter turnout.
This is the command I use:
2013 Feb 21
0
Odd Error message with rare events logit
I'm running a rare events logit model in R using the Zelig package and am
getting some of the oddest error messages that I can't figure out.
y = rnorm(100)
x = c(rep("0",1), rep("1",99))
d = data.frame(won=x, bid=y)
d
mod1 <- zelig(y~x, model="relogit", data=d, tau=1/100,
case.correct="prior", bias.correct=TRUE, robust=TRUE)
> mod1 <-