##Q1. confint.glm(...) fails for an example of HSAUR
data(womensrole, package = HSAUR);
## summary(womensrole);
womensrole_glm_2 - glm(fm2, data = womensrole,family = binomial())
## summary(womensrole_glm_2);
confint(womensrole_glm_2);
## ---Fail-
# Waiting for profiling to be done...
[EMAIL PROTECTED] wrote:
Hello, I was just wondering if you had found a solution? I am having
the same difficulty of converting pdf's into plain text documents in
R. I originally thought I could use the readLines() function, but as
you can see below that did not work.
Why the hell should it?
Dear R Users,
I am having a weird problem. I have three zoo time series, foo, bar and
baz. I run a simple linear regression with foo as the dependent and
bar+baz as independents. Even though the regression runs fine, summary
seems to fail.The code is below. I am happy to send the data along.
Hello,
* I read in a server weblog with read.table.
- OK.
* I look for the downloaded-size-values (filesize of the download)
- OK
* I found - and wanted to substitute them with 0 and
used: weblog$V8[ weblog$V8 == - ] - 0
- OK
* checked the contents on - vs. 0 and found all -
* when then looking at str(weblog),
the - will stay in the levels, mentioned for the variable weblog$V8
- BAD!
Is this snormal behaviour?
Yes, it is. The idea is that a factor has a given set of levels
independent of how often you find them in your data - including
the case that a
Zitat von Weiss, Bernd [EMAIL PROTECTED]:
Philipp Pagel schrieb:
* when then looking at str(weblog),
the - will stay in the levels, mentioned for the variable
weblog$V8
- BAD!
Is this snormal behaviour?
Yes, it is. The idea is that a factor has a given set of levels
Tom,
On Fri, Nov 14, 2008 at 4:35 PM, Tom Backer Johnsen [EMAIL PROTECTED] wrote:
Gabor:
Actually, I meant a script GUI management of ugraph I managed to locate on
the net.
Oh, ok, it will be part of igraph from the next version. But its
capabilities are very limited. We are considering
On Sun, Nov 16, 2008 at 02:52:10PM +0100, Oliver Bandel wrote:
OK, but I thought, when touching the data, it will
recalculate the levels. Now I see, it does not.
No it doesn't - for the reasons given in my explanation.
x - factor(c('A','B','C','A','C'))
y - x[x!='C']
y
[1] A B A
On Sun, Nov 16, 2008 at 3:01 PM, Gábor Csárdi [EMAIL PROTECTED] wrote:
Tom,
On Fri, Nov 14, 2008 at 4:35 PM, Tom Backer Johnsen [EMAIL PROTECTED] wrote:
[...]
suspicion that if the relations between two vertices is non-symmetric, only
one of the is shown. Is it possible to have a pair of,
Hi
I would like to implement the Earth Movers Distance into R and found
the code in C on the following page:
http://www.cs.duke.edu/~tomasi/software/emd.htm
Unfortunately, I don't know much (nearly nothing) about C and how to
call it from R.
As far as I have seen, the calling of the code from R
On Nov 16, 2008, at 9:25 AM, Philipp Pagel wrote:
snip
As you are reading the data from a file anyway, the simplest
solution would probably be to use the colClasses argument ot
read.table in order to get numeric avlues in the first place.
Or use stringsAsFactors = FALSE,
If you have a
See the manual we wrote to help you, 'Writing R Extensions'. If that is
not enough, the topic is covered in several books on R, including 'S
Programming'. You can also find on-line tutorials of various qualities
(often obsolete or inaccurate).
On Sun, 16 Nov 2008, Rainer M Krug wrote:
Hi
Many thanks to Henrik and Spencer, Professor Ripley, and others.
The trick that I had missed (found in 'Writing R Extensions') is that the
CITATION file must be located in the 'inst' subdirectory of the package
sources.
Again, my gratitude to all.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone:
On Sun, 16 Nov 2008, Tolga Uzuner wrote:
I believe there is a recently introduced inconsistency between timeSeries and
zoo which is causing a problem with rbind. I had previously reported that I
was having problems with rbind in the following code:
library(zoo)
Philipp Pagel p.pagel at wzw.tum.de writes:
[...]
As you are reading the data from a file anyway, the simplest
solution would probably be to use the colClasses argument ot
read.table in order to get numeric avlues in the first place.
[...]
Hey, I tried this colClasses-option.
It's really
Try upgrading to R 2.8.0 patched. This works for me
using R 2.8.0 patched from Nov 10th:
library(zoo)
z - 1:10
x - z*z
y - x*z
lm(z ~ x + y)
summary(lm(z ~ x + y))
packageDescription(zoo)$Version
[1] 1.5-4
R.version.string # Vista
[1] R version 2.8.0 Patched (2008-11-10 r46884)
On Sun, Nov
Hello,
I got a hint to use the gmane-interface...
...nice tool. :)
But the gmane-server which is mentioned ont he healp-pages of
http://www.r-project.org/
http://dir.gmane.org/gmane.comp.lang.r.general
and the links on that page also,
just only bring me to an gamne-archive, which is outdated.
Philipp Pagel p.pagel at wzw.tum.de writes:
[...]
foo$bar - factor(foo$bar)
This was my first attemot, before posting here,
and it somehow did not worked...
...now it works so I maybe was too tired,
when trying it and messed something up. :(
[...]
x - factor(c(0,1,3,4,5,7))
x
Dear R users,
From today, you can find a cran task view on probability
distributions I maintain cf. http://cran.r-project.org/web/views/Distributions.html
. The Devel version can be found on R-forge along
with other CTVs (cf.
On 16 November 2008 at 15:51, Prof Brian Ripley wrote:
| See the manual we wrote to help you, 'Writing R Extensions'. If that is
| not enough, the topic is covered in several books on R, including 'S
| Programming'. You can also find on-line tutorials of various qualities
| (often obsolete
I never said it *should* work.
I was simply trying something out that works on other types of files
I've needed in the past (eg: html, csv, dat, etc.). I don't know the
details of the pdf format, but I thought it was worth a try, certainly
no harm in experimenting, and hence I learned that pdfs
R has introduced a new function xtfrm and in order for zoo to
work with it there must be an xtfrm zoo method. The development
version of zoo has such a method but its not yet released. Try this:
xtfrm.zoo - coredata
and then run your code.
On Sun, Nov 16, 2008 at 12:20 PM, Tolga Uzuner
Hei R Users,
i have the following dataframe:
Datetime Temperature and many more
collumns
12008-6-1 00:00:00 5
22008-6-1 02:00:00 5
32008-6-1 03:00:00 6
42008-6-1 03:00:00 0
52008-6-1 04:00:00
Gábor Csárdi wrote:
Tom,
On Fri, Nov 14, 2008 at 4:35 PM, Tom Backer Johnsen [EMAIL PROTECTED] wrote:
Gabor:
Actually, I meant a script GUI management of ugraph I managed to locate on
the net.
Oh, ok, it will be part of igraph from the next version. But its
capabilities are very limited. We
Gábor Csárdi wrote:
On Sun, Nov 16, 2008 at 3:01 PM, Gábor Csárdi [EMAIL PROTECTED] wrote:
Tom,
On Fri, Nov 14, 2008 at 4:35 PM, Tom Backer Johnsen [EMAIL PROTECTED] wrote:
[...]
suspicion that if the relations between two vertices is non-symmetric, only
one of the is shown. Is it possible
Gábor Csárdi wrote:
On Sun, Nov 16, 2008 at 3:01 PM, Gábor Csárdi [EMAIL PROTECTED] wrote:
Tom,
On Fri, Nov 14, 2008 at 4:35 PM, Tom Backer Johnsen [EMAIL PROTECTED] wrote:
[...]
suspicion that if the relations between two vertices is non-symmetric, only
one of the is shown. Is it possible
When I look at that page and choose the all messages with excerpted
text link, I get what appears reasonably up-to-date entries compared
to what I get in my mail subscription. Have you considered the
possibility that it is not Gmane which is without update but rather
your knowledge of
Antje -
I may be missing something, but I usually do this with the negation of
duplicated instead of unique.
So, as an example:
test - data.frame(a = rep(1:5, each = 2),
b = rep(1:5, each = 2), c = rnorm(10))
test[!duplicated(test[c(a, b)]), ]
Hope that helps!
Erik
Antje Nöthlich
This should do it for you:
x - read.table(textConnection( Date time
Temperature
+ 12008-6-1 00:00:00 5
+ 22008-6-1 02:00:00 5
+ 32008-6-1 03:00:00 6
+ 42008-6-1 03:00:00 0
+ 52008-6-1 04:00:00 6
+ 6
On Sun, Nov 16, 2008 at 7:14 PM, Tom Backer Johnsen [EMAIL PROTECTED] wrote:
Gábor Csárdi wrote:
On Sun, Nov 16, 2008 at 3:01 PM, Gábor Csárdi [EMAIL PROTECTED]
wrote:
Tom,
On Fri, Nov 14, 2008 at 4:35 PM, Tom Backer Johnsen [EMAIL PROTECTED]
wrote:
[...]
suspicion that if the
Hi
Martin Maechler wrote:
GF == Guillaume Filteau [EMAIL PROTECTED]
on Thu, 13 Nov 2008 22:38:04 -0500 writes:
GF Hello all, I'm trying to use text symbols in R, such are
GF greek letters and mathematical operators. Do you know
GF where I can find a list telling me how to
Hi R users,
Imagine X as a data frame:
X=read.table(textConnection(
Var1 Var2
1 A H
2 B K
3 A H
4 B H
5 C L
),header=TRUE)
closeAllConnections()
XX is a data frame
Var1 Var2
1AH
2BK
3AH
4BH
5CL
Try:
as.data.frame.matrix(table(X))
On Sun, Nov 16, 2008 at 1:28 PM, CE.KA [EMAIL PROTECTED] wrote:
Hi R users,
Imagine X as a data frame:
X=read.table(textConnection(
Var1 Var2
1 A H
2 B K
3 A H
4 B H
5 C L
),header=TRUE)
closeAllConnections()
X
Hi,
Uwe -- ahh, thank you kindly, I was able to do a web search after
reading your post above in order to find a guide on how to set the
path in windows (i wasn't aware that this is how a file is made
avaiable to the system). I haven't got it to work yet, but at least
i'm on the right track! also
Try putting PATH under name, and the directory path (not the file)
under value. That looks more appropriate to me...
Kind regards
Joris Meys
On Sun, Nov 16, 2008 at 8:41 PM, Tony Breyal [EMAIL PROTECTED] wrote:
Hi,
Uwe -- ahh, thank you kindly, I was able to do a web search after
reading
There is information on http;//batchfiles.googlecode.com
on setting your PATH.
On Sun, Nov 16, 2008 at 2:41 PM, Tony Breyal [EMAIL PROTECTED] wrote:
Hi,
Uwe -- ahh, thank you kindly, I was able to do a web search after
reading your post above in order to find a guide on how to set the
path
On Mon, 17 Nov 2008, Paul Murrell wrote:
Martin Maechler wrote:
GF == Guillaume Filteau [EMAIL PROTECTED]
on Thu, 13 Nov 2008 22:38:04 -0500 writes:
GF Hello all, I'm trying to use text symbols in R, such are
GF greek letters and mathematical operators. Do you know
GF where I
Tony Breyal wrote:
Hi,
Uwe -- ahh, thank you kindly, I was able to do a web search after
reading your post above in order to find a guide on how to set the
path in windows (i wasn't aware that this is how a file is made
avaiable to the system). I haven't got it to work yet, but at least
i'm
Hi Gabor, yes, i used the link from that website to figure out the
steps to setting my path: http://www.computerhope.com/issues/ch000549.htm
but i'm still doing something wrong it seems (see my last post in
response to Joris).
Cheers,
Tony Breyal
On 16 Nov, 20:07, Gabor Grothendieck [EMAIL
Hi Joris, there is already a variable called 'Path', therefore i
appended the directory path to the other strings already in the value
section:
Name: Path
Value: %SystemRoot%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem;
%SystemRoot%\system32\nls;%SystemRoot%\system32\nls\ENGLISH;C:\Program
Hi Uwe, that is what i have now done, following the sugestion from
Joris, but unfortuantly that hasn't worked.
Cheers,
Tony Breyal
Ps. i would like to appologise, i am replying to these posts through
google.groups.co.uk and there is a delay beween me posting a reply and
that post appearing in
The PATH is a list of directories, not filenames. It will look
for files in those directories.
The easiest method discussed on the batchfiles site is to
use the Redmond utility.
You may or may not need to exit R and start it up again once you
have set your path.
If your programs are self
Zitat von David Winsemius [EMAIL PROTECTED]:
When I look at that page and choose the all messages with excerpted
text link, I get what appears reasonably up-to-date entries compared
to what I get in my mail subscription. Have you considered the
possibility that it is not Gmane which is
Many thanks Gabor, as always, much appreciated.
Regards,
Tolga
Gabor Grothendieck wrote:
R has introduced a new function xtfrm and in order for zoo to
work with it there must be an xtfrm zoo method. The development
version of zoo has such a method but its not yet released. Try this:
Oliver Bandel wrote:
Hello,
I got a hint to use the gmane-interface...
...nice tool. :)
But the gmane-server which is mentioned ont he healp-pages of
http://www.r-project.org/
http://dir.gmane.org/gmane.comp.lang.r.general
and the links on that page also,
just only bring me to an
Hi Tony,
You can name several variables 'Path' without problems. So you best
restore the original variable PATH to its original value (or it ain't
going to work any more) and just add a new one, call that PATH as
well, and add the directory C:\Program Files\xpdf , like Uwe already
suggested.
Zitat von Peter Dalgaard [EMAIL PROTECTED]:
Oliver Bandel wrote:
Hello,
I got a hint to use the gmane-interface...
...nice tool. :)
But the gmane-server which is mentioned ont he healp-pages of
http://www.r-project.org/
http://dir.gmane.org/gmane.comp.lang.r.general
and the
On Nov 16, 2008, at 5:05 PM, Oliver Bandel wrote:
Zitat von David Winsemius [EMAIL PROTECTED]:
When I look at that page and choose the all messages with excerpted
text link, I get what appears reasonably up-to-date entries compared
to what I get in my mail subscription. Have you considered
Deepayan,
Thanks again for your explanations. I must admit I am still
struggling with some of these concepts (though I have in fact read
more than just help pages!), and I would be delighted to find a
holistic introduction to the concepts of formulae and of the lattice
system -- in fact I have
Dear R users...
I made this by help of one of R users.
_
X=matrix(seq(1,4), 2 , 2)
B=matrix(c(0.6,1.0,2.5,1.5) , 2 , 2)
func - function(i,y0,j) { y0*exp(X[i,]%*%B[,j]) }
list1 - expand.grid( i=c(1,2) , y0=c(1,2) , j=c(1,2) )
Sampling with and without replacement
I seem unable to use replace = F when I want to sample without
replacement. I would think
that it comes down to F is not a legitimate abbreviation for FALSE.
except that
Dalgaard (p. 118) uses F for FALSE and it works
pairwise.t.test(folate, ventilation,
Dear Kathie,
Try this:
X=matrix(seq(1,4), 2 , 2)
B=matrix(c(0.6,1.0,2.5,1.5) , 2 , 2)
list1 - expand.grid( i=c(1,2) , y0=c(1,2) , j=c(1,2) )
apply(list1,1,function(x) {
i=x[1]
y0=x[2]
j=x[3]
y0*exp(X[i,]%*%B[,j])
On 17/11/2008, at 1:56 PM, David C. Howell wrote:
Sampling with and without replacement
I seem unable to use replace = F when I want to sample without
replacement. I would think
that it comes down to F is not a legitimate abbreviation for FALSE.
except that
Dalgaard (p. 118) uses F for FALSE
It is better programming practice to use FALSE for false and TRUE for
true, and not F and T. This is because it is quite legal to do this:
T - FALSE
F - TRUE
or any other assignment. If you re-assign T or F (which are set to TRUE
and FALSE at the beginning of a session), you run into the sort of
Deemed colleagues
I would appreciate your help with a sentence to calculate the linear
distance between 2 geographical points (coordinates in UTM).
In advance thnks for your attention,
--
Ricardo Bandin Llanos
[EMAIL PROTECTED]
Estudiante - Magíster Cs. m. Pesquerías
Universidad de Concepción,
Here are a couple of hits I got from searching the archives. Any of them help:
http://finzi.psych.upenn.edu/R/library/spectralGP/html/rdist.earth.html
http://finzi.psych.upenn.edu/R/library/RFOC/html/GreatDist.html
On Sun, Nov 16, 2008 at 8:40 PM, Ricardo Bandin [EMAIL PROTECTED] wrote:
On Sun, 16 Nov 2008, kathie wrote:
func needs to return a result of the same length as the input vectors
try this:
func2 - function(i,y0,j) { y0* exp( rowSums(X[i,] * t(B[,j])) ) }
do.call(func2, list1 )
Chuck
Dear R users...
I made this by help of one of R users.
Hi,
I have a data set similar to the following
State Gender Quantity
TX Male1
NY Female 2
TX Male3
NY Female 4
I need to calculate cumulative sum of the quantity by State and Gender. The
expected output is
State Gender QuantityCumQuantity
TX
I have a dataset which has categories A-G (7 categories) with corresponding
x and y variables. I want to delete one category from the dataset and run a
logit model on only 6 categories instead of 7. How do I delete a category? I
have tried the omit command but it just returns TRUE and False
On 11/16/08, Stavros Macrakis [EMAIL PROTECTED] wrote:
Deepayan,
Thanks again for your explanations. I must admit I am still
struggling with some of these concepts (though I have in fact read
more than just help pages!), and I would be delighted to find a
holistic introduction to the
Deepayan --- you are too patient by half! :-)
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
Thank you for the pointer to your paper Getting Started with Lattice
Graphics -- it looks very useful.
I think we're talking past each other on the question of the semantics
of formula operators, and it's probably not productive to continue. I
would guess that the underlying issue is that we
on 11/16/2008 08:12 PM jeffc wrote:
Hi,
I have a data set similar to the following
State Gender Quantity
TXMale1
NYFemale 2
TXMale3
NYFemale 4
I need to calculate cumulative sum of the quantity by State and Gender. The
expected output is
State Gender
Hi Stavros,
I think we're talking past each other on the question of the semantics
of formula operators, and it's probably not productive to continue. I
would guess that the underlying issue is that we come from different
communities, and so have different background assumptions. Perhaps
Hello R-Community,
I am pretty new to R and I am fascinated what R can do! I am doing
phylogenetic analysis in R, and my current project includes two problems
that I am unable to solve, unfortunately. I am hoping that someone knows
an easy, and R-typical solution, any help would be
Try this. The first line appends the cumulative sum column and the second
displays it in sorted fashion:
DF$cumQuantity - ave(DF$Quantity, DF$State, DF$Gender, FUN = cumsum)
DF[order(DF$State, DF$Gender), ]
On Sun, Nov 16, 2008 at 9:12 PM, jeffc [EMAIL PROTECTED] wrote:
Hi,
I have a data
On Sun, Nov 16, 2008 at 9:56 PM, Christian Arnold [EMAIL PROTECTED] wrote:
Hello R-Community,
I am pretty new to R and I am fascinated what R can do! I am doing
phylogenetic analysis in R, and my current project includes two problems
that I am unable to solve, unfortunately. I am hoping that
Hello R-Community,
I am a rookie in R and I am fascinated with the power of bio computing
by R. I am analysing gene expression data from Real time PCR. I have
used absolute gene quantitation to measure gene copy number in all my
transcripts. All my data has been normalised them to a housekeeping
Dear Jeff,
Try also
df=df[order(df$State, df$Gender), ]
df$cQuantity-unlist(tapply(df[,3],df[,-3],cumsum))
df
State Gender Quantity cQuantity
2NY Female2 2
4NY Female4 6
1TX Male1 1
3TX Male3 4
HTH,
Jorge
Thanks Gabor and Baptiste,
Both your suggestions give me the results I was looking for!
Robert
CC: r-help@r-project.org
From: [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Subject: Re: [R] Lattice: groups and no groups with panel.superpose
Date: Sat, 15 Nov 2008 16:59:31 +
Hi,
I
70 matches
Mail list logo