Hi,
May be this helps.
dat1- structure(list(trialnumber = 1:9, Control = c(0L, 1L, 0L, 0L,
1L, 0L, 1L, 0L, 0L), HC = c(0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L
), LC = c(1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L), pleasantness = c(3L,
3L, 2L, 2L, 3L, 4L, 3L, 3L, 2L)), .Names = c(trialnumber, Control,
HC,
On 10/01/2013 08:17 AM, Xianwen Chen wrote:
Dear fellows,
The two questions are on sm.density.compare(). I compare kernel density
estimates of two arrays of data.
I'd like to scale up y-axis so that I can show better the differences in
y values. English is not my first language so I'll try to
Hi,
Not sure if this helps.
#example dataset
set.seed(24)
mat1-matrix(sample(1:50,20*20,replace=TRUE),ncol=20)
set.seed(49)
colnames(mat1)- sample(rep(LETTERS[1:4],5),20)
rownames(mat1)- colnames(mat1)
mat2-combn(LETTERS[1:4],2)
res- sapply(split(mat2,col(mat2)),function(x){
Hi
I'm ChungHoon Cho and graduate student of Yonsei University in south korea.
I just wanna know one thing.
Can I use R for implementing GBI(Graph-Based Induction)?
and is there any sample code for that?
I will wait for your response.
thank you.
[[alternative HTML version
I am fitting a model with both nested and non-nested factors. I would really
appreciate any comments on what is the correct code to fit the model with lmer.
The data structure is a bit complicated. My unit of analysis is the
village-organization-dyad. The outcome is observed at the level of the
Hi Jean,
nbsp;
I tried to follow the introductions on statmethods... Didn#39;t really work.
nbsp;
Do you think it might have sth to do with hownbsp;my dataset is built up?
nbsp;
Thanks for your help!
nbsp;
Tobias
nbsp;
nbsp;
nbsp;
Gesendet:nbsp;Montag, 30. September 2013 um
On 10/01/2013 05:53 PM, happyR wrote:
...
Do you think it might have sth to do with hownbsp;my dataset is built up?
Yes, happyR, it is definitely your dataset. I think this one is far
superior:
spineless-matrix(c(6,3,9,5,4,6,7,1,3,8,7,6,10,4,3,7,6,11),
nrow=2)
Dear R-Users,
A new book Learning Regression Analysis
by Simulation (Springer) is out.
http://www.springer.com/statistics/statistical+theory+and+methods/book/978-4-431-54320-6
There is a slight problem with that second link. It is the same as the first.
It is easier to sent the data (or a decent sample of it) in the email by using
dput. Type ?dput for more information. Essentially you just do dput(myfile)
copy it and paste it into the email. This allow the reader
I have this weird data set with 2 predictors and one dependent variable -
attached.
predictor1 has all zeros except for one 1.
I am runnning a simple logistic regression:
temp-read.csv(x data for reg224.csv)
myreg - glm(dv~predictor1+predictor2,data=temp,
family=binomial(logit))
google complete separation logistic
-- Bert
On Tue, Oct 1, 2013 at 7:41 AM, Dimitri Liakhovitski
dimitri.liakhovit...@gmail.com wrote:
I have this weird data set with 2 predictors and one dependent variable -
attached.
predictor1 has all zeros except for one 1.
I am runnning a simple
Hi,
I did have warning messages about convergence issues using binomial GLM
with logit link with my data in the past
Do you detect separation using the function separation.detection{brglm}?
Regards,
Xochitl C.
Xochitl CORMON
+33 (0)3 21 99 56 84
Doctorante en sciences
Thanks so much for your help everyone. This really helped me a lot.
On Sun, Sep 29, 2013 at 11:42 PM, Barry Rowlingson
b.rowling...@lancaster.ac.uk wrote:
On Sun, Sep 29, 2013 at 10:48 PM, john doe anon.r.u...@gmail.com wrote:
I am having trouble understanding how classes in R work. Here
Thank you very much, Bert - it's very helpful.
This post says that R issues a warning:
Warning message:*glm.fit: fitted probabilities numerically 0 or 1 occurred *
However, in my case there is no warning. How could I detect complete
separation in my data? I need to be able to flag it in my
Xochitl CORMON
Le 01/10/2013 17:29, Dimitri Liakhovitski a écrit :
Thank you very much, Bert - it's very helpful.
This post says that R issues a warning:
Warning message:
*glm.fit: fitted probabilities numerically 0 or 1 occurred
*
Actually the warning message should be something
Ah, thank you very much - I did not understand first brglm was the name of
a package!
Dimitri
On Tue, Oct 1, 2013 at 11:34 AM, Xochitl CORMON
xochitl.cor...@ifremer.frwrote:
Xochitl CORMON
Le 01/10/2013 17:29, Dimitri Liakhovitski a écrit :
Thank you very much, Bert - it's very
Le 01/10/2013 17:41, Dimitri Liakhovitski a écrit :
Ah, thank you very much - I did not understand first brglm was the name
of a package!
Dimitri
My bad!
If there is separation you should see it in the way the coefficient
diverges from one (it's pretty exponential). You can increase the
hey lobachevsky - thanks!
but: Error in barplot.default (matrix, col = c(quot;yellowquot;,
quot;purplequot;), beside = TRUE) :
#39;height#39; must be a vector or a matrix
what to do? :)
Gesendet:nbsp;Dienstag, 01. Oktober 2013 um 10:42 Uhr
Von:nbsp;quot;Jim Lemon
Hello,
I have a question concerning the output of leveneTest. I don't understand
the 7 in my output
Levene's Test for Homogeneity of Variance (center = median)
Df F value Pr(F)
group 2 0.0173 0.9829
7# Where does this number come from?
Thanks.
Hermann
res
group.ID
Can anyone please guide me on any useful links or resource regarding text
mining in R?
--
Regards
Umesh Khatri
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
Dear Experts
I'm very new to R and after some days or reading and testing I tried to make
my first small application (and failed ...)
In general I would like to work with sqldf , ggplot2 to create some graphical
output.
At present I got stuck with this:
PROG
Hello everybody,
i have to save a 100 iteration computation in a file every 5 iterations
until the end.
I first give a vector A of 100 elements for the 100 iterations and i want
to update A every 5 iterations.
I use save but it doesn't work.
Someone has an idea, i need a help
Cheers.
--
R-help,
There are several packages that have package 'XML' as a dependency. I cannot
get this package (XML) to extract or install. I have tried manually
downloaded to local machines (client PC and Linux), unsuccessful extraction.
How can this be compiled manually? How can a new version
Dear all,
Sorry for bothering you but I am having a doubt for a long time and I cannot
find a solution for it. I don't know if you can help me but I am doing a
research study on Correspondence Analysis and I have been using the package
ANACOR. After performing the Correspondence Analysis, I
Dear Hermann,
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Hermann Norpois
Sent: Tuesday, October 01, 2013 8:07 AM
To: r-help
Subject: [R] leveneTest - output
Hello,
I have a question concerning the output of
Do you know about task views? Try
http://cran.r-project.org/web/views/NaturalLanguageProcessing.html
Best,
Ista
On Tue, Oct 1, 2013 at 6:06 AM, umesh khatri khatriumes...@gmail.com wrote:
Can anyone please guide me on any useful links or resource regarding text
mining in R?
--
Regards
I want to show over dispersion in my data by plotting the variance against
the mean for 8 groups. I used the following code to plot by type
(parasite).
scatterplot(variance~mean | parasite, reg.line=lm, smooth=FALSE,
spread=FALSE, boxplots=FALSE, span=0.5, by.groups=TRUE, data=Disp2)
I now want
Hi,
Not sure if this helps.
deviceSummary - data.frame(item = character(0) , value = numeric(0) , unit =
character(0) ,stringsAsFactors=FALSE)
newlst- list(primitiveSpace,1.1,mm2)
deviceSummary[nrow(deviceSummary)+1,]- newlst
newlst2- list(primitiveSpace,2.2,mm2)
Hello,
The main problem is in the way you form newRow. You can't mix data
classes in a vector created with c(), so all its elements become
characters, the least common denominator:
newRow - c(primitiveSpace, 1.1 , mm2)
newRow
[1] primitiveSpace 1.1mm2
Then when you rbind it
Hello,
Just to add that the value 0.9829 [Pr(F)] can be checked with
pf(0.0173, df1 = 2, df2 = 7, lower.tail = FALSE)
Rui Barradas
Em 01-10-2013 19:02, John Fox escreveu:
Dear Hermann,
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On
R running on my MacPro under OS X 10.6.8 no longer installs packages. I
would greatly appreciate help in sorting this out.
Typically I get the messages like:
install.packages(party)
Warning: unable to access index for repository
Thanks Jim,
The picture width is adjusted. That was great help! I really appreciate it.
Does lwd mean 'line width'? I'm not exactly sure how that will help. Can
you please explain the intuition behind it?
Can I also ask another question? My graph is made from 7 plots from
On Tue, Oct 1, 2013 at 2:05 PM, carlisle thacker
carlisle.thac...@gmail.com wrote:
R running on my MacPro under OS X 10.6.8 no longer installs packages. I
would greatly appreciate help in sorting this out.
Typically I get the messages like:
install.packages(party)
Warning: unable to access
On 01/10/2013 2:05 PM, carlisle thacker wrote:
R running on my MacPro under OS X 10.6.8 no longer installs packages. I
would greatly appreciate help in sorting this out.
Typically I get the messages like:
install.packages(party)
Warning: unable to access index for repository
On 10/02/13 00:22, Steven Dwayne Randolph wrote:
R-help,
There are several packages that have package 'XML' as a dependency. I
cannot get this package (XML) to extract or install. I have tried manually
downloaded to local machines (client PC and Linux), unsuccessful extraction.
How
I'm used to using ctrl-c to end operations without killing R. But I've used
xlsx in this session, which loads Java, which apparently intercepts the
ctrl-C. Accordingly, I hit ctrl-C, R died, and I lost a lot of work.
I did some looking, and found a
On 10/01/2013 07:32 PM, happyR wrote:
hey lobachevsky - thanks!
but: Error in barplot.default (matrix, col =
c(quot;yellowquot;,quot;purplequot;), beside = TRUE) :
#39;height#39; must be a vector or a matrix
what to do? :)
First, one must always avoid confusing the general case with a
I am running a large mixed model, 65k entries on 11 fixed effects and one
random. One of the fixed effects is Course a factor that takes on 14
different values
levels(Master.complete$Course)
[1] B101 B2A B2B B2C C118A C118B C118C C2A C2B
[10] C2C N101 P7A P7B P7C
and another
On 10/02/2013 04:43 AM, Xianwen Chen wrote:
Thanks Jim,
The picture width is adjusted. That was great help! I really appreciate it.
Does lwd mean 'line width'? I'm not exactly sure how that will help. Can
you please explain the intuition behind it?
Yes, lwd refers to line width. The
I'm running R version 3.0.2 (2013-09-25) on a Windows computer.
Here is a simplified version of my data:
a - expand.grid(SITE=c(A,B),TREAT=c(low,med,high), REP=c(1:5))
VAR - round(rep(c(10,12,14,15,17,18),5)+rnorm(30)*2,1)
dat - cbind(a,VAR)
I'm running a GAM on my data (package: mgcv, version
40 matches
Mail list logo
