Dear Greg,
thanks for your answer, but that is not the point.
My question is another: from the data in the R prompt how
can I compare the models?
I don´t understand the output of anova(fit1, fit2, fit3)
What I have to look to understand something about the comparison of the models?
Look the
Dear Robert,
thanks so much!!! Now I understand!
So you also think that I have to check only the residuals and not the data
directly.
Now just for curiosity I did the the shapiro test on the residuals. The problem
is that on fit3 I don´t get from the test
that the data are normally
When introduced to R, I learned how to use *apply whenever I could to avoid
for-loops and all. And, getting the habit, I think I somehow got the
mis-conception that it is a magic source, always an optimal way of coding in
R.
Thanks a lot for all of your helpful advice and comment!
Young
On Wed,
Dear Greg,
thanks so much. Now I understand!!
Thanks!!
From: Greg Snow greg.s...@imail.org
r-help@r-project.org
Sent: Wed, January 5, 2011 11:27:53 PM
Subject: RE: [R] Comparing fitting models
The output with all three fits gives you 2 comparisons, fit1
On Wed, Jan 05, 2011 at 09:19:49PM +, Anthony Staines wrote:
Dear colleagues,
This may be a question with a really obvious answer, but I
can't find it. I have access to a large file with real
medical record identifiers (mixed strings of characters and
numbers) in it. These represent
Dear R users,
I used to be able to use the package bipartite to calculate network metrics
using the function networklevel, but now (with the new version 1.14) I keep
getting the error messages below. I don't understand why I am getting these
error messages given that colSums is not 2, and the
Dear All,
I have a problem in understanding how the interactions of 2 ways ANOVA work,
because I get conflicting results
from a t-test and an anova. For most of you my problem is very simple I am sure.
I need an help with an example, looking at one table I am analyzing. The table
is in
Hello,
I am new to R and need some help.
I have data in following format
DATA matrix
DateTime oh l c
2009-01-01 07:30:00 2 3 4 5
2009-01-01 07:33:00 4 2 5 7
I am able to fill the gap using combination of seq and chron
On Wed, Jan 5, 2011 at 10:49 PM, Young Cho young.s...@gmail.com wrote:
When introduced to R, I learned how to use *apply whenever I could to avoid
for-loops and all. And, getting the habit, I think I somehow got the
mis-conception that it is a magic source, always an optimal way of coding in
On 05.01.2011 22:49, Young Cho wrote:
When introduced to R, I learned how to use *apply whenever I could to avoid
for-loops and all. And, getting the habit, I think I somehow got the
mis-conception that it is a magic source, always an optimal way of coding in
R.
That is right, but your apply
Hello,
I have applied all tips (except moving to different DB lib) :
move to R-2.12.1
try R CMD INSTALL instead of install.packages('ROracle');
run as root
checked that I have full 32 bit env
It still fails with same error at installation when trying to load ROracle
lib
-
Dear help-list,
I have a problem with timeSequence {timeDate}.
When I use it like
timeSequence(from = 2008-01-01, to = 2010-12-13, by = 1 month)
GMT
[1] [2008-01-01] [2008-02-01] [2008-03-01] [2008-04-01] [2008-05-01]
[2008-06-01] [2008-07-01] [2008-08-01] [2008-09-01] [2008-10-01]
Hi:
Look into the na.approx() function in package zoo. The discussion below may
be of help:
http://r.789695.n4.nabble.com/Filling-in-missing-time-samples-with-na-approx-td3063682.html
HTH,
Dennis
On Wed, Jan 5, 2011 at 10:17 PM, Rustamali Manesiya rmanes...@gmail.comwrote:
Hello,
I
Hi,
I have several matrix in a list, for example:
e
[[1]]
[,1] [,2]
[1,]13
[2,]24
[[2]]
[,1] [,2]
[1,]14
[2,]25
[3,]36
[[3]]
[,1] [,2]
[1,]21
I would like to join them by column i.e.
[,1] [,2] [,3] [,4][,5] [,6]
[1,]1
try this:
matLis - list(matrix(1:4, 2, 2), matrix(1:6, 3, 2),
matrix(2:1, 1, 2))
n - max(sapply(matLis, nrow))
do.call(cbind, lapply(matLis, function (x)
rbind(x, matrix(, n-nrow(x), ncol(x)
I hope it helps.
Best,
Dimitris
On 1/6/2011 11:56 AM, emj83 wrote:
Hi,
I have
Excellent- that is just what I need. Thank you so much for your prompt help,
Emma
--
View this message in context:
http://r.789695.n4.nabble.com/How-to-join-matrices-of-different-row-length-from-a-list-tp3177212p3177252.html
Sent from the R help mailing list archive at Nabble.com.
Hello everyone and season's greetings.
I have an array that looks like that R-c(1,2,3,4,5,6).
Is it possible to select all the elements but except one? For example to not
select the third element and get back
(1,2,4,5,6)?
How can I do that?
I would like to thank you in advance for your help
Alaios alaios at yahoo.com writes:
Hello everyone and season's greetings.
I have an array that looks like that R-c(1,2,3,4,5,6).
Is it possible to select all the elements but except one?
For example to not select the third element and get back
(1,2,4,5,6)?
How can I do that?
I would
Please report bugs in contributed packages to the corresponding package
maintainer.
Uwe Ligges
On 06.01.2011 10:30, Benjamin B. wrote:
Dear help-list,
I have a problem with timeSequence {timeDate}.
When I use it like
timeSequence(from = 2008-01-01, to = 2010-12-13, by = 1 month)
GMT
Hi all,
I read in a text book, that you can examine a variable that is colinear
with others, and giving different ANOVA output and explanatory power
when ordered differently in the model forula, by modelling that
explanatory variable, against the others colinear with it. Then, using
that
Hopefully a quick question. My package has a NEWS.Rd file that is not
being found by news.
The news function calls tools:::.build_news_db which has this line:
nfile - file.path(dir, inst, NEWS.Rd)
So it appears that the news function is searching for
mypackage/inst/NEWS.Rd.
However, Writing R
Dear Benjamin,
timeSequence() ultimately is relying on seq.POSIXt(). If you look at
the last paragraph of the Details section in ?seq.POSIXt it seems to
basically indicate that using by = month just sequences through the
months and the day is only changed if it is invalid for a particular
month
If you look at tools:::.build_news_db, the plain text NEWS file is
searched for in pkg/NEWS and pkg/inst/NEWS, but NEWS.Rd in only
searched for in pkg/inst/NEWS.Rd.
Looks like a bug to me.
I *think*.
Thanks,
Kevin
On Thu, Jan 6, 2011 at 7:09 AM, Kevin Wright kw.s...@gmail.com wrote:
Peter,
thank you, that's what I was looking for!
David, I forgot to tell you my OS. Sorry... it's Win7. I'm running a RKWard
session.
And this is strange:
Cairo(example.pdf, type=pdf,width=12,height=12,units=cm,dpi=300)
Error: could not find function Cairo
... maybe you're not using the
Hi there
I want to test the expression of a subset of genes for correlation with
patient survival. I found out that the coxph function is appropriate for
doing this since it works with continuous variables such as Affy mRNA
expression values.
I applied the following code:
cp -
No one replied to this, so I'll try again, with a simple example. I
calculate a set of log odds ratios, and turn them into a data frame as
follows:
library(vcdExtra)
(lor.CM - loddsratio(CoalMiners))
log odds ratios for Wheeze and Breathlessness by Age
25-2930-3435-3940-44
On Thu, Jan 6, 2011 at 5:27 AM, Joshua Wiley jwiley.ps...@gmail.com wrote:
timeSequence() ultimately is relying on seq.POSIXt(). If you look at
My apologies, I spoke nonsense---timeSequence() does NOT rely on
seq.POSIXt(). The timeDate package has its own method defined for
seq() which is what
Dear list,
This must be an easy one. I have a data frame like this one:
test.df - data.frame(x1=c(2,3,5), x2=c(5, 3, 4), w=c(0.8, 0.3, 0.5))
and I want to construct a weighted mean of the first two columns using the
third column for weighting; i.e.
y[1] = x1[1]*w[1] + x2[1]*(1-w[1])
y[2] =
Dear List,
I have a data frame called trait with roughly 800 species in, each species have
15 columns of information:
Species 1 2 3 etc..
a t y h
b f j u
c r y u
Dear List,
I have a data frame called trait with roughly 800 species in, each species have
15 columns of information:
Species 1 2 3 etc..
a t y h
b f j u
c r y u
Hi Vassilis,
Try
test.df$y - with(test.df, x1*w + x2*(1-w))
test.df
HTH,
Jorge
On Thu, Jan 6, 2011 at 8:33 AM, Vassilis wrote:
Dear list,
This must be an easy one. I have a data frame like this one:
test.df - data.frame(x1=c(2,3,5), x2=c(5, 3, 4), w=c(0.8, 0.3, 0.5))
and I want to
Hi list,
Please tell me why sometimes RSymphony.dll crashes on
Windows 7 (may be on other Windows too).
I saw that it depends on
input data. For example this code bring to crash:
library('Rsymphony')
mat = c(0, 1, -100.37967, 0, 0, 1, 0, 0, -200, 0,
0, 1, -0.4, 0, 0, 1, 0, 0, -0.4, 0, 0.5,
Dear R users,
I've installed Rserve for R version 2.11.0 on x64 Windows 7. I've added the
Rserve_d and Rserve files to the /bin/ folder where the R.exe is installed
in the program files. I have also created an Rserv.cfg file that contains
the following text:
remote enable
auth disable
ear ALL,
The last part of my thesis analysis is the cross validation. Right now I am
having difficulty using the cross validation of gstat. Below are my commands
with the tsport_ace as the variable:
nfold - 3
part - sample(1:nfold, 69, replace = TRUE)
sel - (part != 1)
m.model - x2[sel, ]
Em 6/1/2011 11:33, Vassilis escreveu:
Dear list,
This must be an easy one. I have a data frame like this one:
test.df - data.frame(x1=c(2,3,5), x2=c(5, 3, 4), w=c(0.8, 0.3, 0.5))
and I want to construct a weighted mean of the first two columns using the
third column for weighting; i.e.
y[1]
Dear all,
I have a data frame that is created like that
data.frame(x=CRX[-1],y=CRY[-1],z=CRagent[[1]]$sr)
the output looks like
45 116 162 -30.89105988567164
46 128 79 -42.66296679571184
47 180 195 -30.45626175641315
48 114 83 -45.26843476475688
49 118 73 -46.85389245327003
How can I select
Is there any common variable? From your description, I don't see how you would
link a species to a community. I mean if you select species a in df1 how would
you know what community it is in?
--- On Thu, 1/6/11, Chris Mcowen chrismco...@gmail.com wrote:
From: Chris Mcowen
With dataframe xx (naming a data.frame as data.frame is a bit dicey
subset(xx, xx[,4] -45)
--- On Thu, 1/6/11, Alaios ala...@yahoo.com wrote:
From: Alaios ala...@yahoo.com
Subject: [R] Find and remove elemnts of a data frame
To: r-help@r-project.org
Received: Thursday, January 6, 2011,
On Jan 6, 2011, at 6:23 AM, Dimitris Rizopoulos wrote:
try this:
matLis - list(matrix(1:4, 2, 2), matrix(1:6, 3, 2),
matrix(2:1, 1, 2))
n - max(sapply(matLis, nrow))
do.call(cbind, lapply(matLis, function (x)
rbind(x, matrix(, n-nrow(x), ncol(x)
It's good that you solved the OP's
On 06.01.2011 15:04, Vladyslav Kolbasin wrote:
Hi list,
Please tell me why sometimes RSymphony.dll crashes on
Windows 7 (may be on other Windows too).
I saw that it depends on
input data. For example this code bring to crash:
library('Rsymphony')
mat = c(0, 1, -100.37967, 0, 0, 1, 0, 0,
You're definitely doing something I don't understand.
Here's a quick mock-up of the sample data sets you provided. How do you know
if species a in trait is from NA1102 or AT1302 in com? I may be blind but I
just don't see what you're matching when you use which. Perhaps you might want
to
On Jan 6, 2011, at 6:36 AM, Chris Mcowen wrote:
Dear List,
I have a data frame called trait with roughly 800 species in, each
species have 15 columns of information:
Species 1 2 3 etc..
a t y h
b f j
Hi,
Sorry the formatting has messed up, the common variable is species ( spelt
different in both but i have corrected this now)
I am having some luck with the which function but am struggling to automate it
for all communities
On 6 Jan 2011, at 15:09, John Kane wrote:
Is there any common
Thanks a lot
worked fine
Regards
--- On Thu, 1/6/11, John Kane jrkrid...@yahoo.ca wrote:
From: John Kane jrkrid...@yahoo.ca
Subject: Re: [R] Find and remove elemnts of a data frame
To: r-help@r-project.org, Alaios ala...@yahoo.com
Date: Thursday, January 6, 2011, 3:27 PM
With dataframe xx
Hello,
Can the x and y axis limits be specified in a density plot with the
mixtools package for a finite mixture model? Uncommenting the xlim2/
ylim2 lines in the plot command below generates 'not a graphical
parameter' warnings (and does not change the axis settings), and
uncommenting the
On Jan 6, 2011, at 3:11 AM, thomas.car...@bnpparibas.com wrote:
Hello,
I have applied all tips (except moving to different DB lib) :
move to R-2.12.1
try R CMD INSTALL instead of install.packages('ROracle');
run as root
checked that I have full 32 bit env
It still fails with same
Hello, after calculating a multinomial logit regression on my data, I
compared the output to an output retrieved with SPSS 18 (Mac). The
coefficients appear to be the same, but the logLik (and therefore fit)
values differ widely. Why?
The regression in R:
set.seed(1234)
df - data.frame(
On Thu, Jan 6, 2011 at 5:56 AM, emj83 stp08...@shef.ac.uk wrote:
Hi,
I have several matrix in a list, for example:
e
[[1]]
[,1] [,2]
[1,] 1 3
[2,] 2 4
[[2]]
[,1] [,2]
[1,] 1 4
[2,] 2 5
[3,] 3 6
[[3]]
[,1] [,2]
[1,] 2 1
I would
On Jan 6, 2011, at 11:06 AM, Sören Vogel wrote:
Hello, after calculating a multinomial logit regression on my data, I
compared the output to an output retrieved with SPSS 18 (Mac). The
coefficients appear to be the same, but the logLik (and therefore fit)
values differ widely. Why?
The
Sören Vogel sovo0815 at gmail.com writes:
Hello, after calculating a multinomial logit regression on my data, I
compared the output to an output retrieved with SPSS 18 (Mac). The
coefficients appear to be the same, but the logLik (and therefore fit)
values differ widely. Why?
Since
Thanks for your replies. I am no mathematician or statistician by far,
however, it appears to me that the actual value of any of the two LLs
is indeed important when it comes to calculation of
Pseudo-R-Squared-s. If Rnagel devides by (some transformation of) the
actiual value of llnull then any
Look forward to it.
Thanks.
Shige
On Sat, Jan 1, 2011 at 8:45 AM, Michael Lawrence
lawrence.mich...@gene.com wrote:
Please watch for 2.20.5 and let me know if it helps. Not really sure what is
going on here, but someone else has reported the same issue.
Thanks,
Michael
On Wed, Dec 29,
You need RGtk2 2.20.7 which is now on CRAN. Others have seen this,
but it has taken a while to track down the exact cause.
The diagnosis was that ML used a recent GNU tar which created a
tarball with hard links that R's untar was not prepared to deal with.
We consider that is a bug in GNU
On Jan 6, 2011, at 11:23 AM, Sören Vogel wrote:
Thanks for your replies. I am no mathematician or statistician by far,
however, it appears to me that the actual value of any of the two LLs
is indeed important when it comes to calculation of
Pseudo-R-Squared-s. If Rnagel devides by (some
Dear David,
Thats great, thanks very much for the help, much appreciated.
On 6 Jan 2011, at 15:53, David Winsemius wrote:
On Jan 6, 2011, at 6:36 AM, Chris Mcowen wrote:
Dear List,
I have a data frame called trait with roughly 800 species in, each species
have 15 columns of
This question of mine is now solved, thanks to a suggestion by Homer Strong,
the organizer of the R user group in Portland, Oregon.
The unmet dependencies as reported by install was caused by an incorrect
entry in /etc/apt/sources.list. Previously I had
deb
Santanu,
I second Phil's suggestion. sas.get is actually quite nice.
Another current option is using a command-line utility called dsread
(http://www.oview.co.uk/dsread/) to convert the sas7bdat file to a csv or tsv
format, which can then easily be read into R using read.table and its
'merge' comes in handy:
spec - read.table(textConnection(Species 1 2 3
+ a t y h
+ b f j u
+ c r y u), header=TRUE)
comm - read.table(textConnection(community species
You really need to spend more time with a good aov textbook and probably a
consultant that can explain things to you face to face. But here is a basic
explanation to get you pointed in the right direction:
Consider a simple 2x2 example with factors A and B each with 2 levels (1 and
2). Draw
Has anyone succeeded in porting any PDF compression tools to R so far?
Regards,
Yihui
--
Yihui Xie xieyi...@gmail.com
Phone: 515-294-2465 Web: http://yihui.name
Department of Statistics, Iowa State University
2215 Snedecor Hall, Ames, IA
On Wed, Jan 5, 2011 at 7:34 PM, Andrew Miles
Remember that an non-significant result (especially one that is still near
alpha like yours) does not give evidence that the null is true. The reason
that the 1st 2 tests below don't show significance is more due to lack of power
than some of the residuals being normal. The only test that I
Dear R forum helpers,I have following datatrans -
data.frame(currency_transacted = c(EURO, USD, USD, GBP, USD, AUD),
position_amt = c(1, 25000, 2, 15000, 22000, 3))date -
c(12/31/2010, 12/30/2010, 12/29/2010, 12/28/2010, 12/27/2010,
12/24/2010, 12/23/2010, 12/22/2010, 12/21/2010,
I think that you are looking for the 'resid' and 'fitted' functions, these will
give you the residuals and fitted values from an lm object (that added together
gives the original response but are orthogonal to each other). Those values
can then be assigned to a data frame or used by
I'm guessing this page will answer
your question:
http://www.portfolioprobe.com/2010/10/04/a-tale-of-two-returns/
If not, then you need to be more
specific.
On 06/01/2011 18:07, Amelia Vettori wrote:
Dear R forum helpers,I have following datatrans- data.frame(currency_transacted = c(EURO, USD,
Hi,
I wrote a simple script to retrieve an n number of followers for a given
user in Twitter. I used a sample of n=10 to test my script and worked
perfectly but once I started to changes n I started to get the following
error:
Error in list_to_dataframe(res, attr(.data, split_labels)) :
I've got a dataset which looks like this in the beginning:
cbr dust smoking expo
1 0 0.20 15
2 0 0.25 14
3 0 0.25 18
4 0 0.25 14
5 0 0.25 14
(till no. 1240, anyway, a huge set)
I have to analyse cbr and
Folks:
The following has NOTHING (obvious) to do with R. But I believe that
all on this list would find it relevant and, I hope, informative. It
is LONG. I apologize in advance to those who feel I have wasted their
time.
http://www.newyorker.com/reporting/2010/12/13/101213fa_fact_lehrer
Best
Hi Everyone,
Thanks a lot for your guys help, I finally got it running, now I can call my
R function from my java code.
But there is one problem left.
In my R code , it open a window and plot a curve.
Once running in Java, the window did successfully opened, but nothing
inside, instead
Michael Friendly frien...@yorku.ca
on Thu, 06 Jan 2011 09:33:25 -0500 writes:
No one replied to this, so I'll try again, with a simple example. I
calculate a set of log odds ratios, and turn them into a data frame as
follows:
library(vcdExtra)
(lor.CM -
Yes, exactly. But the problem is with NEWS.Rd, not NEWS.
pkg/inst/NEWS.Rd is moved to pkg/NEWS.Rd at build time, but for
installed packages, news tried to load pkg/inst/NEWS.Rd.
I'm going to file a bug report.
Kevin
On Thu, Jan 6, 2011 at 7:29 AM, Kevin Wright kw.s...@gmail.com wrote:
If
Given the data structure below and the call to ggplot2, how can I increase the
size of the axis scale points, the line weight, and the size of the legend?
ddata -structure(list(year = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c(2003, 2007), class = factor),
On Jan 6, 2011, at 1:23 PM, kiotoqq wrote:
I've got a dataset which looks like this in the beginning:
cbr dust smoking expo
1 0 0.20 15
2 0 0.25 14
3 0 0.25 18
4 0 0.25 14
5 0 0.25 14
(till no. 1240,
On 06/01/2011 2:19 PM, Kevin Wright wrote:
Yes, exactly. But the problem is with NEWS.Rd, not NEWS.
I'm not sure who you are arguing with, but if you do file a bug report,
please also put together a simple reproducible example, e.g. a small
package containing NEWS.Rd in the inst directory
Dear Greg,
thanks so much, I think that now I have understood. Please confirm me this
reading what follows ;-)
To summarize from the beginning, the table I analyzed is the result of a simple
experiment. Subjects where exposed to some stimuli
and they where asked to evaluate the degree of
I used chisq.test(read.table(C:/Users/Maggy/Downloads/dust.asc,
header=TRUE))
and got this
Pearson's Chi-squared test
data: read.table(C:/Users/Maggy/Downloads/dust.asc, header = TRUE)
X-squared = 5226.164, df = 3735, p-value 2.2e-16
and I think it should be right for the whole
# Can anyone suggest why this works
datafilename -
http://personality-project.org/r/datasets/maps.mixx.epi.bfi.data;
person.data - read.table(datafilename,header=TRUE)
# but this does not?
dd - https://sites.google.com/site/jrkrideau/home/general-stores/trees.txt;
treedata - read.table(dd,
Ok,
I see ;-)
Let´s put in this way then. When do I have to use the kruskal wallis test? I
mean, when I am very sure that I have
to use it instead of ANOVA?
Thanks
Best regards
P.S. In addition, which is the non parametric methods corresponding to a 2
ways
anova?..or have I to
repeat
I was communicating with Kevin off-list.
The problem seems to be run time, not install time. News() calls
tools:::.build_news_db(), and the 2nd line of that function is:
nfile - file.path(dir, inst, NEWS.Rd)
and that's the problem: an installed package shouldn't have an inst/
subdirectory,
On Jan 6, 2011, at 2:34 PM, kiotoqq wrote:
I used chisq.test(read.table(C:/Users/Maggy/Downloads/dust.asc,
header=TRUE))
So, where did you download this data and when is your homework due?
and got this
Pearson's Chi-squared test
data:
Some would argue to always use the kruskal wallis test since we never know for
sure if we have normality. Personally I am not sure that I understand what
exactly that test is really testing. Plus in your case you are doing a two-way
anova and kruskal.test does one-way, so it will not work for
Andy, thanks for providing a clear way of saying it. I thought I was
clear in the first place, but oh well).
Here is the structure of my source files:
hwpkg/DESCRIPTION
hwpkg/R/hw.R
hwpkg/inst/NEWS.Rd
I'm using Windows XP. When I install this package, I do this:
Rcmd INSTALL hwpkg
Which
See inline
From: Frodo Jedi [mailto:frodo.j...@yahoo.com]
Sent: Thursday, January 06, 2011 12:37 PM
To: Greg Snow; r-help@r-project.org
Subject: Re: [R] Problem with 2-ways ANOVA interactions
Dear Greg,
thanks so much, I think that now I have understood. Please confirm me this
reading
David,
I think the poster wants to use one of the columns as x and the other as y,
ignoring the remaining columns. If that is the case then he/she needs to read
the section in Introduction to R on subsetting data frames.
I agree that the output so far is meaningless, from the degrees of
I don't know. I can access it from within OOo Calc so it might be an authority
problem but it seems unlikely or would Google sites have some strange
restrictions ?
I did discover that the trees.txt file was rather messed up so I tried again
with a clean csv file and I'm still getting the same
On 06/01/2011 20:29, Greg Snow wrote:
Some would argue to always use the kruskal wallis test since we never know
for sure if we have normality. Personally I am not sure that I understand
what exactly that test is really testing. Plus in your case you are doing a
two-way anova and
Hi:
On Thu, Jan 6, 2011 at 5:36 AM, Eduardo de Oliveira Horta
eduardo.oliveiraho...@gmail.com wrote:
Peter,
thank you, that's what I was looking for!
David, I forgot to tell you my OS. Sorry... it's Win7. I'm running a RKWard
session.
And this is strange:
Cairo(example.pdf,
--- On Thu, 1/6/11, Alberto Negron albertoneg...@gmail.com wrote:
From: Alberto Negron albertoneg...@gmail.com
Subject: Re: [R] Accessing data via url
To: John Kane jrkrid...@yahoo.ca
Cc: R R-help r-h...@stat.math.ethz.ch
Received: Thursday, January 6, 2011, 4:20 PM
there is an option in
Hi:
I'm not sure what you mean by line weight - I interpreted it to mean size,
but you already have have it in your geom_line() call. If you want thicker
lines, increase size. As for the others, try this:
ggplot(ddata, aes(x = area, y = disc, colour = year)) +
geom_point() +
Well, the goal is to include a reference to f2 in the error message
returned by f2('char'); sys.call(1) appears to do the trick. You
mentioned this function could be unreliable, could you please provide an
example?
Sebastien
William Dunlap wrote:
-Original Message-
From:
Dear R-users,
Is there a way I can prevent global variables to be visible within my
functions?
Sebastien
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
On 06/01/2011 4:45 PM, Sebastien Bihorel wrote:
Dear R-users,
Is there a way I can prevent global variables to be visible within my
functions?
Yes, but you probably shouldn't. You would do it by setting the
environment of the function to something that doesn't have the global
environment
Dear all,
I´m specifying the fisher.exact test for use with summary.formula as follows:
u-function(a,b){
j-fisher.test(a)
p-list(P=j$p.value,stat=NA,df=NA,testname=j$method,statname=)
return(p)
}
However I´m also required to specify stat df.
Hi, I am seeking help on designing an algorithm to identify the locations of
stretches of 1s in a vector of 0s and 1s. Below is an simple example:
dat-as.data.frame(cbind(a=c(F,F,T,T,T,T,F,F,T,T,F,T,T,T,T,F,F,F,F,T)
,b=c(4,12,13,16,18,20,28,30,34,46,47,49,61,73,77,84,87,90,95,97)))
dat
a
On Thu, Jan 6, 2011 at 4:59 PM, Duncan Murdoch murdoch.dun...@gmail.com wrote:
On 06/01/2011 4:45 PM, Sebastien Bihorel wrote:
Dear R-users,
Is there a way I can prevent global variables to be visible within my
functions?
Yes, but you probably shouldn't. You would do it by setting the
The closest I get is
u-function(a,b){
j-fisher.test(a)
p-list(P=j$p.value,stat=1,df=1,testname=j$method,statname=)
return(p)
}
However then I manually have to edit the output. Is there a smart way of doing
this?
You're not explaining what
sys.call(1) will work now, but as my code evolves it tends
towards having more layers of function calls and so the count
given in sys.call(count) has to be altered. This is a maintainance
problem. Sometimes, as when using generic functions, you are
not even aware that the number of calls between
Does the prtest argument help when you actually use the 'print' function
around your summary.formula object? I think that's how I
solve it.
I.e.,
sf1 - summary(trt~sex+ascites,data=ex,test=T,method=reverse,catTest=u)
print(sf1, prtest = P)
Descriptive Statistics by trt
Allright..Works like a charm. However I do believe that the prtest vector
should have been mentioned in the catTest or conTest option. Appreciate your
time and effort.
Best,
//M
On 6. jan. 2011, at 23.24, Erik Iverson wrote:
Does the prtest argument help when you actually use the 'print'
The next week's New Yorker has some decent rebuttal letters. The case
is hardly as clear-cut as the author would like to believe.
Carl
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
try this:
?rle
Carl
**
From: array chip arrayprofile_at_yahoo.com
Date: Thu, 06 Jan 2011 14:16:38 -0800 (PST)
Hi, I am seeking help on designing an algorithm to identify the
locations of stretches of 1s in a vector of 0s and 1s. Below is an
simple example:
1 - 100 of 144 matches
Mail list logo