On 18/10/2013 15:01, Vincent Guyader wrote:
Hi everyone,
If I compute a Ordered Logistic or Probit Regression with the polr
function from MASS package. the summary give me : coefficients, Standard
error and Tvalue.. but not directly the p.value.
I can compute manualy the Pvalue, but Is there
require(MASS)
n- 10
n
[1] 1e+05
w- rgamma(n,rate=0.5,shape=5);head(w)
[1] 14.135679 18.839336 6.623502 17.165505 2.847572 7.511726
fit.w- fitdistr(w,weibull)
Warning messages:
1: In densfun(x, parm[1], parm[2], ...) : NaNs produced
2: In densfun(x, parm[1], parm[2], ...) : NaNs
I want to estimate the following fixed effect model:
y_i,t = alpha_i + beta_1 x1_t + beta_2 x2_i,tx2_i,t = gamma_i + gamma_1
x1_t + gamma_2 Z1_i + gamma_3 Z2_i
I can use ivreg from AER to do the iv regression.
fm - ivreg(y_i,t ~ x1_t + x2_i,t |x1_t + Z1_i + Z2_i,
data =
On Sat, 19 Oct 2013, Ignacio Martinez wrote:
I want to estimate the following fixed effect model:
y_i,t = alpha_i + beta_1 x1_t + beta_2 x2_i,tx2_i,t = gamma_i + gamma_1
x1_t + gamma_2 Z1_i + gamma_3 Z2_i
I can use ivreg from AER to do the iv regression.
fm - ivreg(y_i,t ~ x1_t + x2_i,t |
Le vendredi 18 octobre 2013 à 13:27 -0400, Earl Brown a écrit :
Thanks Duncan. However, now I can't get the Spanish and Portuguese accented
vowels to come out correctly and still keep the indents in the saved
document, even when I set encoding = UTF-8:
library(XML)
concepts - c(español,
Thanks.
One more question. When I do summary to see the coefficients how can i tell
R to print just some of the coefficients?
Thanks again
On Sat, Oct 19, 2013 at 7:29 AM, Achim Zeileis achim.zeil...@uibk.ac.atwrote:
On Sat, 19 Oct 2013, Ignacio Martinez wrote:
I want to estimate the
The summary currently doesn't provide this feature. But you can use something
like coeftest(obj)[1:3,] or coeftest(obj)[myreg,] or something like that...
(sent from mobile phone)
Original message
From: Ignacio Martinez ignaci...@gmail.com
Date:
To: Achim Zeileis
At 19:05 18/10/2013, William Dunlap wrote:
However, I do not understand the
substitute(...())
idiom. Would you care to explain it? (No is an acceptable answer!).
I don't completely understand it either, I treat it as an idiom. I
saw it on this list once.
Bill Dunlap
Spotfire, TIBCO
On Oct 19, 2013, at 2:26 AM, Charles Thuo wrote:
require(MASS)
n- 10
n
[1] 1e+05
w- rgamma(n,rate=0.5,shape=5);head(w)
[1] 14.135679 18.839336 6.623502 17.165505 2.847572 7.511726
fit.w- fitdistr(w,weibull)
Warning messages:
1: In densfun(x, parm[1], parm[2], ...) : NaNs produced
2013/10/19 Prof Brian Ripley rip...@stats.ox.ac.uk
On 18/10/2013 15:01, Vincent Guyader wrote:
Hi everyone,
If I compute a Ordered Logistic or Probit Regression with the polr
function from MASS package. the summary give me : coefficients, Standard
error and Tvalue.. but not directly the
On Sat, 19 Oct 2013, Michael Dewey wrote:
At 19:05 18/10/2013, William Dunlap wrote:
However, I do not understand the
substitute(...())
idiom. Would you care to explain it? (No is an acceptable answer!).
I don't completely understand it either, I treat it as an idiom. I saw it
on
Hi everyone!
I am making some time series models, and as i want to compare a lot of
models I thought it would be smart to compare the AIC, AICc and BIC values
from the models. My question is, how can I extract the BIC and AICc from the
model?
As an example:
kings -
On 13-10-18 1:54 PM, Bert Gunter wrote:
Yes, similar, but better, as match.call() will get unwanted named
arguments, too.
However, I do not understand the
substitute(...())
idiom. Would you care to explain it? (No is an acceptable answer!).
I suspect it's a bug, though I can't see it's one
On Oct 19, 2013, at 3:44 AM, Chris89 wrote:
Hi everyone!
I am making some time series models, and as i want to compare a lot of
models I thought it would be smart to compare the AIC, AICc and BIC values
from the models. My question is, how can I extract the BIC and AICc from the
model?
Dear R users,
Suppose I want to randomly generate some data, in matrix form, randomly
swap some of the elements and calculate trace of the matrix for each of
these stages. If the value of trace obtained in the later is bigger than
the former, drop the latter matrix and go back to the former
On Oct 19, 2013, at 10:32 AM, David Winsemius wrote:
On Oct 19, 2013, at 3:44 AM, Chris89 wrote:
Hi everyone!
I am making some time series models, and as i want to compare a lot of
models I thought it would be smart to compare the AIC, AICc and BIC values
from the models. My question
Hello,
Em 19-10-2013 19:07, David Winsemius escreveu:
On Oct 19, 2013, at 10:32 AM, David Winsemius wrote:
On Oct 19, 2013, at 3:44 AM, Chris89 wrote:
Hi everyone!
I am making some time series models, and as i want to compare a lot of
models I thought it would be smart to compare the
Hello,
Seems simple.
fun - function(n = 10){
matd - matrix(sample(1:30,30, replace=FALSE), ncol=5, nrow=6)
res - list(mat = NULL, Design_best = matd, Original_design = matd)
trace - sum(diag(matd))
res$mat - rbind(res$mat, c(trace = trace, iterations = 0))
Vincent,
I believe Prof. Ripley is referring to this:
http://www.stats.ox.ac.uk/pub/MASS4/
Contact
Details:---
Contact me: tal.gal...@gmail.com |
Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) |
Thank you so very much!
It works like a charm !!!
Regards,
Laz
On 10/19/2013 5:03 PM, Rui Barradas wrote:
fun - function(n = 10){
matd - matrix(sample(1:30,30, replace=FALSE), ncol=5, nrow=6)
res - list(mat = NULL, Design_best = matd, Original_design = matd)
trace - sum(diag(matd))
Dear useRs,I have the following data z of two variables x(z[,1]) and
y(z[,2]). dput(z)structure(c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14,
15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34,
35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51,
Dear useRs,
I have the following data z of two variables x(z[,1]) and y(z[,2]).
dput(z)
structure(c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,
20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39,
40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50,
Follow the posting guide! Code code code.
(What you call a dotplot I think is what most of us call a
scatterplot. And there are many ways of doing this in R and packages,
and we have no idea what you used.)
Cheers,
Bert
On Sat, Oct 19, 2013 at 2:29 PM, eliza botto eliza_bo...@hotmail.com
Dear R users,
There is a small problem which I don't know how to sort it out, based on
the former example I had explained earlier own.
I am calling my own three functions which are based on simulations as below:
library(gmp)
library(knitr) # load this packages for publishing results
Dear R users,
Dear R users,
(I had not included two more functions in the previous mail. This
version is complete)
There is a small problem which I don't know how to sort it out, based on
the former example I had explained earlier own.
I am calling my own functions which are based on
Duncan, like what happens with Milan, when I print to the console both the
accents and the indents come through correctly. However, when I save to a file,
the indents disappear:
# when saved to a file the indents disappear
saveXML(doc, file = test.xml, encoding = UTF-8)
# contents of text.xml
Hi,
No problem.
dat1 -
subset(head(read.table(temp1.txt,skip=6,sep=,,header=FALSE,fill=TRUE,stringsAsFactors=FALSE),-1),select=1:4)
dat1$V1 - as.numeric(dat1$V1)
str(dat1)
#'data.frame': 96409 obs. of 4 variables:
# $ V1: num 891231 900101 900101 900101 900101 ...
# $ V2: int 2400 100
I know I cannot be the only one who is not able to install the XML package
from CRAN (zip or tar file) Many packages depend on this XML package. Can
someone help me either access the source for a good binary? I have received no
assistance from the author/developer of the package.
Steven
Even when I attempt to unzip the file on windows, I get an error that says its
not a compressed file.
From: Rolf Turner [mailto:rolf.tur...@xtra.co.nz]
Sent: Tuesday, October 01, 2013 4:23 PM
To: Steven Dwayne Randolph
Cc: r-help@r-project.org
Subject: Re: [R] 'XML' package cannot be un-zipped
On 10/20/2013 08:29 AM, eliza botto wrote:
Dear useRs,
...
I made a dot plot between them. what i want to do is to bold the size of
dots,compared to that of others, for the following ranges
---80x100 and 2y4
---40x60 and 8y10
---100x120 and 6y8
Then, i also want to draw square
Dear Denis and jim,Thankyou for your help.
eliza
Date: Sun, 20 Oct 2013 10:55:05 +1100
From: j...@bitwrit.com.au
To: eliza_bo...@hotmail.com
CC: r-help@r-project.org
Subject: Re: [R] bold dot size and name in plot
On 10/20/2013 08:29 AM, eliza botto wrote:
Dear useRs,
...
I made
On 13-10-19 3:47 PM, Steven Dwayne Randolph wrote:
I know I cannot be the only one who is not able to install the XML package
from CRAN (zip or tar file) Many packages depend on this XML package. Can
someone help me either access the source for a good binary? I have received no
assistance
On Sat, Oct 19, 2013 at 8:23 PM, Duncan Murdoch
murdoch.dun...@gmail.com wrote:
On 13-10-19 3:47 PM, Steven Dwayne Randolph wrote:
I know I cannot be the only one who is not able to install the XML
package from CRAN (zip or tar file) Many packages depend on this XML
package. Can someone
Dear All,
For a report that I am writing, I need to create a histogram plot with
x-axis fixed as very bad, bad, fair, good and very good - i.e. the order
not changed. Can someone give me an example? For sample purposes, I am
giving the following data: 159, 374, 3765, 11388, 6708. Thanks,
On 10/20/2013 11:54 AM, kan...@iitb.ac.in wrote:
Dear All,
For a report that I am writing, I need to create a histogram plot with
x-axis fixed as very bad, bad, fair, good and very good - i.e. the order
not changed. Can someone give me an example? For sample purposes, I am
giving the
Hi,
May be this also helps. Using ?ggplot()
indx - 1 + 2*((z[,1]80 z[,1] 100) (z[,2]2 z[,2] 4)) + 4*((z[,1]40
z[,1] 60) (z[,2]8 z[,2] 10)) + 8*((z[,1]100 z[,1] 120)
(z[,2] 6 z[,2] 8))
indx[indx1] - 2
z1-cbind(z,indx)
z2 - as.data.frame(z1)
colnames(z2)[1:2] - c(x,y)
rect1 -
Dear All,
opt = list()
opt$aa - TRUE
Error in opt$aa - TRUE : object 'opt' not found
Why?
--
Regards,
Taiyun
--
Taiyun Wei weitai...@gmail.com
Homepage: http://blog.cos.name/taiyun/
Phone: +86-15201142716
Renmin University of China
__
37 matches
Mail list logo