Hi,
I'm just learning about poison links for the glm function.
One of the data sets I'm playing with has several of the variables as
factors (i.e. month, group, etc.)
When I call the glm function with a formula that has a factor variable,
R automatically converts the variable to a series of
On Tue, 2 Mar 2010 08:58:25 +1300 Peter Alspach
wrote:
> This brings up another confusion for new users. Simply typing the
> object name at the command line gives just one view of the object (that
> provided by print()).
Good point. Any good introduction to R should include a brief discussion
Quoting Uwe Ligges :
R is not able to re-encode the file to the native encoding. But if you
keep it in UTF-8, what is the problem to grep for the specific
characters (as grep and friends support the argument useBytes these
days)?
The Problem with UTF-8 is that I'm not able to cat a valid xml-
I have been reading binary files, and parsing the output, for some
time now. I have tried to develop a technique that is as robust as
possible to all the strange things that appear in text fields, not to
mention different global/regional encodings. I have no control over
the data generated by use
On 02-Mar-10 08:02:27, Noah Silverman wrote:
> Hi,
> I'm just learning about poison links for the glm function.
>
> One of the data sets I'm playing with has several of the
> variables as factors (i.e. month, group, etc.)
>
> When I call the glm function with a formula that has a factor
> variabl
Emil Davtyan wrote:
Hello -
I need to do simple linear autoregressive model with R software for my
thesis. I looked into all your documentation and I am not able to find
anything too helpful. Can someone help me with the codes?
Thanks
Emil
[[alternative HTML version deleted]]
__
On Mon, 1 Mar 2010, Emil Davtyan wrote:
Hello -
I need to do simple linear autoregressive model with R software for my
thesis. I looked into all your documentation and I am not able to find
anything too helpful. Can someone help me with the codes?
By "all documentation" you mean that you ha
On Mon, 2010-03-01 at 16:55 -0500, Diana Connett wrote:
> Hello.
> I just downloaded R onto a new computer, and after entering library(MASS), I
> still get the message "Error: could not find function "predict.lda"" when I
> try
> to use the predict.lda function (even just "predict.lda()")
>
> Can
On Mon, 1 Mar 2010, FMH wrote:
Dear All,
I'm trying to find changepoints in a data series which only consist of 11
measurements of altitude(x) and temperature(y), respectively, in which the data
are as followed:
y = 16.3, 16.2, 16.1, 15.6, 14.2, 10, 8.2, 8.0, 7.5, 7.3, 7.2
x = 1, 2, 5,
Thanks, Berwin. That works just great!
--sundar
On Tue, Mar 2, 2010 at 12:57 AM, Berwin A Turlach
wrote:
> G'day Sundar,
>
> On Mon, 1 Mar 2010 23:46:55 -0800
> Sundar Dorai-Raj wrote:
>
> > Thanks for the input, but I don't want "try" in the Sweave output. I
> > want the output to look just li
Hello John,
As you said, I could also take a means model and test linear hypothesis for
the desired effects - would this also be the case for the repeated measure i
did in the first place.
I copied the model from the car model where you first call:
> modx<-lm(cbind(div_h, div_l) ~ site, divrep
On Mon, 1 Mar 2010 14:23:04 -0500 Liaw, Andy
wrote:
> Indeed this is one of the (few, I believe) traps of R,
Oh, no; there are many more:
http://www.burns-stat.com/pages/Tutor/R_inferno.pdf :-)
--
Karl Ove Hufthammer
__
R-help@r-project.org mailing
On Mon, 01 Mar 2010 16:26:37 - (GMT) ted.hard...@manchester.ac.uk
wrote:
> In vim (to which I'm wedded
> for life) it will pick up matching (), {} and [].
You can also easily move between matching delimeter by typing '%'. A
similar feature should be available in all *good* text editors.
--
On Mon, 01 Mar 2010 10:00:07 -0500 Duncan Murdoch
wrote:
> Suppose X is a dataframe or a matrix. What would you expect to get from
> X[1]? What about as.vector(X), or as.numeric(X)?
All this of course depends on type of object one is speaking of. There
are plenty of surprises available, and
Dear R users
Today I discovered that function lm.influence() stops when applied to
glm objects with the following error message
Error in if (NROW(e) != n) stop("non-NA residual length does not match
cases used in fitting") :
argument is of length zero
After inspecting lm.influence.R (both into
Ted,
Brilliant explanation (as usual)
I'm back in school, just starting on a post-graduate degree in stats so
the help is really appreciated.
Now, I have a slightly trickier question about the same model.
I've seen more than one way to get "values" out of the glm model.
i.e. If we're looki
The dsread output is little-endian, as that's the native format for
floats on the Wintel platform. The byte order should stay the same if
converting directly to a float, using a data structure like (C/C++):
union {
char bytes[8];
double value;
}
If reading the values with a SAS HEX i
G'day Sundar,
On Mon, 1 Mar 2010 23:46:55 -0800
Sundar Dorai-Raj wrote:
> Thanks for the input, but I don't want "try" in the Sweave output. I
> want the output to look just like it does in the console, as if an
> uncaptured error really did occur.
I don't think that you will get around using "
>From CRAN:
"2.8 What's the best way to upgrade?
That's a matter of taste. For most people the best thing to do is to
uninstall R (see the previous Q), install the new version, copy any
installed packages to the library folder in the new installation, run
update.packages(checkBuilt=TRUE, ask=FAL
I am trying to load package fpc in order to use the 'plotcluster' function
however everytime I attempt to do so I get the following warning message:
> library(fpc)
Loading required package: MASS
Error: package 'MASS' could not be loaded
In addition: Warning messages:
1: pac
On Mon, Mar 1, 2010 at 11:49 PM, Liviu Andronic wrote:
> On 3/1/10, Keo Ormsby wrote:
>> Perhaps my biggest problem was that I couldn't (and still haven't) seen
>> *absolute beginners* documents.
>>
> there was once a link posted on r-sig-teaching that would probably fit
> your needs, but I cann
On Mon, Mar 1, 2010 at 10:49 AM, Duncan Mackay wrote:
> Dear All
>
> Below is a toy example of a modified standard bwplot.
>
> require(lattice)
> DF <-
> data.frame(site = rep(1:5, each = 20),
> height = rnorm(100))
>
> bwplot(site ~ height,DF,
> pch = "|",
> par.settings = list(strip.ba
Hi,
I have some vector v1,v2,...,vk, with different lengths. I want to
consider these vectors as a matrix with k rows.
Can you please guide me how I can do it?
Regards
khazaei
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r
What were your biggest misconceptions or
stumbling blocks to getting up and running
with R?
Easy. I terms of materials I have been unable to find good books that
introduce users to R from the perspective of someone familiar only
with packages like SPSS or STATA, or not familiar with statist
Please take what follows not as an ad hominem statement, but rather as an
attempt to improve what is already an excellent program, that has been built as
a result of many, many hours of dedicated work by many, many unpaid, unsung
volunteers.
It troubles me a bit that when a confusing aspect of
Hi,
I hope somebody can give me an idea where can I can find the code for empirical
copula.
I have a bivariate data.
Thank you so much for your help.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https:
Brandon Zicha wrote:
What were your biggest misconceptions or
stumbling blocks to getting up and running
with R?
Easy. I terms of materials I have been unable to find good books that
introduce users to R from the perspective of someone familiar only
with packages like SPSS or STATA, or not f
Hi,
How could i generate random real numbers between 0 en 2*pi?
Thanks,
Frederik
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://ww
Greetings, fellow travelers.
When I paste a series of commands into R, they execute serially, and when I
go back through commands by hitting the up key, they show up as a block,
rather than as individual lines. Is there any way to change this behavior?
I'm running the 32-bit build of R 2.10.1 on M
On Tue, 2 Mar 2010 12:31:45 +0100 Brandon Zicha
wrote:
> Easy. I terms of materials I have been unable to find good books that
> introduce users to R from the perspective of someone familiar only
> with packages like SPSS or STATA,
Have you read these books:
R for SAS and SPSS Users
http:/
On Tue, 2 Mar 2010 12:29:28 +0100 (CET) khaz...@ceremade.dauphine.fr
wrote:
> I have some vector v1,v2,...,vk, with different lengths. I want to
> consider these vectors as a matrix with k rows.
> Can you please guide me how I can do it?
What do you want to do with the missing elements?
--
Ka
John Sorkin wrote:
Please take what follows not as an ad hominem statement, but rather as an
attempt to improve what is already an excellent program, that has been built as
a result of many, many hours of dedicated work by many, many unpaid, unsung
volunteers.
It troubles me a bit that when a
Hi Brandon,
I just read this book, which I am sure you will be interested in:
http://www.amazon.com/SAS-SPSS-Users-Statistics-Computing/dp/0387094172
Cheers!!
Albert-Jan
~~
In the face of ambiguity, refuse the temptation to gue
Hi all,
I have a (hopefully) simple newbie-level question.
# I have data like this:
dtf <- data.frame(read.table(textConnection("var value
company 9887.1
company 91117.0
blaah 91.1
etc 11
etc 97111"), header=TRUE))
# I would like to have output like this (the index number may
On Tue, Mar 2, 2010 at 7:27 AM, Duncan Murdoch wrote:
> John Sorkin wrote:
>>
>> Please take what follows not as an ad hominem statement, but rather as an
>> attempt to improve what is already an excellent program, that has been built
>> as a result of many, many hours of dedicated work by many, m
On Tue, 2 Mar 2010 11:51:39 +0100 frederik vanhaelst
wrote:
> How could i generate random real numbers between 0 en 2*pi?
Ten such numbers from the uniform distribution:
2*pi*runif(10)
--
Karl Ove Hufthammer
__
R-help@r-project.org mailing list
ht
> runif(20,0,2*pi)
[1] 1.29417642 1.10933879 4.31669186 2.41339484 4.83705630 3.12713657
4.50893007 6.23232980 2.38783146 4.88483239 5.87292617
[12] 1.33293077 4.09458703 0.7593 1.67899698 2.42602639 0.08413394
2.40261439 5.46442874 2.13847582
>
On Tue, Mar 2, 2010 at 5:51 AM, frederik vanha
Hello,
is there an elegant way, how I can convert each row of a data frame into
distinct elements of a list?
In essence, what I'm looking for is something like
rows.to.lists <- function( df ) {
ll <- NULL
for( i in 1:nrow(df) )
ll <- append( ll, list(df[i,]) )
This is what you would use a 'list' for.
On Tue, Mar 2, 2010 at 6:29 AM, wrote:
> Hi,
> I have some vector v1,v2,...,vk, with different lengths. I want to
> consider these vectors as a matrix with k rows.
> Can you please guide me how I can do it?
>
> Regards
> khazaei
>
> _
>From what distribution?
If the uniform,
runif(100,0,2*pi)
If another, install package Runuran, and do this
?Runuran
Vignette("Runuran")
HTH
Dr. Rubén Roa-Ureta
AZTI - Tecnalia / Marine Research Unit
Txatxarram
Dear All,
Consider the following trivial code snippet
rm(list=ls())
name_vec <- c("color1", "color2")
pdf("test_color.pdf")
plot(seq(5), seq(5), main=paste(name_vec[1]," and ",name_vec[2], sep=""))
dev.off()
What I would like to achieve is rather simple to explain, but it is
giving me a head
> as.data.frame(t(df))
For example
> x <- as.data.frame(t(mtcars))
> typeof(x)
[1] "list"
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Sebastian Bauer
Sent: Tuesday, March 02, 2010 8:12 AM
To: r-help@r-project.org
Subject: [R]
Lower case vignette("Runuran")
Sorry, it's this clever email program.
R.
Dr. Rubén Roa-Ureta
AZTI - Tecnalia / Marine Research Unit
Txatxarramendi Ugartea z/g
48395 Sukarrieta (Bizkaia)
SPAIN
-Mensaje or
See this thread :
http://finzi.psych.upenn.edu/Rhelp10/2009-January/185693.html
On 03/02/2010 02:18 PM, Lorenzo Isella wrote:
Dear All,
Consider the following trivial code snippet
rm(list=ls())
name_vec <- c("color1", "color2")
pdf("test_color.pdf")
plot(seq(5), seq(5), main=paste(name_vec[1
Hi!
On 03/02/2010 02:22 PM, Nutter, Benjamin wrote:
as.data.frame(t(df))
For example
x<- as.data.frame(t(mtcars))
typeof(x)
[1] "list"
Thanks for the quick reply!
I would never have guessed that as.data.frame() works that way!
BTW
This one seems also to do the trick:
rows.to.list <- f
If memory serves, Bill Venables said in the paper cited several times
here, that there's only one type of sums of squares. So there's only
one type of "ANOVA" (if I understand what you mean by ANOVA).
Just forget about the different types of tests, and simply ask yourself
this (hopefully simple a
Here is an example using proto based on converting Duncan's example:
library(proto)
Stack <- proto(new = function(.) proto(Stack,
stack = NULL,
push = function(., el) .$stack <- c(list(el), .$stack),
pop = function(.) { stopifnot(length(.$stack) > 0)
out <-
Try this:
reshape(cbind(id = as.numeric(dtf$var), dtf, time = with(dtf,
ave(value, var, FUN = seq))), timevar="time", direction="wide")
Or:
xtabs(value ~ var + ave(value, var, FUN = seq), data = dtf)
On Tue, Mar 2, 2010 at 9:40 AM, Albert-Jan Roskam wrote:
> Hi all,
> I have a (hopefully) si
On Mar 2, 2010, at 8:11 AM, Sebastian Bauer wrote:
Hello,
is there an elegant way, how I can convert each row of a data frame
into distinct elements of a list?
split(dfrm, rownames(dfrm))
In essence, what I'm looking for is something like
rows.to.lists <- function( df ) {
ll <-
Ravi Kulkarni wrote:
Hello,
I think I am beginning to understand what is involved in the so-called
"Type-I, II, ..." ANOVAS (thanks to all the replies I got for yesterday's
post). I have a question that will help me (and others?) understand it
better (or remove a misunderstanding):
I know tha
Brandon Zicha wrote:
Hey Paul,
Hey Brandon, (adding R-help in the cc)
I agree with you that the documentation of R could be better, especially
with more examples in code showing not only the common cases, but also
more esoteric cases. It would be great if everyone invested a lot of
time to w
Sorry there were 2 typos in my note:
John Fox's posting provided a lot of information. Briefly, the "Types"
refer to whether effects are adjusted for all other effects in the model
(Types II, III, IV) or no (Type I or sequential tests only adjust for
no -> not
EARLIER terms in the model).
Hi Henrique,
*Thank you!* The reshape code does precisely what I want.
Cheers!!
Albert-Jan
~~
In the face of ambiguity, refuse the temptation to guess.
~~
--
Hello !
I`m using SVMs for multi-class classification problems. Therefore I`m using the
svm() function in the package "e1071".
If I use svm(...type="C-classification") everything works fine. But if I want
to use nu-SVM with svm(..., type="nu-classification", nu=0.5) R crashes
immediately. No er
On Tue, 2010-03-02 at 00:58 -0800, Noah Silverman wrote:
> Ted,
>
> Brilliant explanation (as usual)
>
> I'm back in school, just starting on a post-graduate degree in stats so
> the help is really appreciated.
>
> Now, I have a slightly trickier question about the same model.
>
> I've seen mo
On 02.03.2010 15:41, Häring, Tim (LWF) wrote:
Hello !
I`m using SVMs for multi-class classification problems. Therefore I`m using the svm()
function in the package "e1071".
If I use svm(...type="C-classification") everything works fine. But if I want to use
nu-SVM with svm(..., type="nu-clas
-- Forwarded message --
From: Duncan Murdoch
Date: 25 February 2010 22:35
Subject: Re: [R] Building R packages in Windows 7
To: RHelp , Eric Ferreira
On 25/02/2010 7:56 PM, Eric Ferreira wrote:
> Thank you, Sir, but how can I demand it to create HTML files?
>
The tools::Rd2HTM
> > I`m using SVMs for multi-class classification problems. Therefore I`m
> using the svm() function in the package "e1071".
> > If I use svm(...type="C-classification") everything works fine. But
> if I want to use nu-SVM with svm(..., type="nu-classification", nu=0.5)
> R crashes immediately. No
frederik vanhaelst wrote:
Hi,
How could i generate random real numbers between 0 en 2*pi?
Thanks,
Frederik
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Hello R wizards,
What is the best way to read a data file containing both fixed-width and
tab-delimited files? (More detail follows.)
_*Details:*_
The U.S. Bureau of Labor Statistics provides local area unemployment
statistics at ftp://ftp.bls.gov/pub/time.series/la/, and the data are
document
The indexing in xts is very nice; it may do what you want.
library(xts)
x.xts <- as.xts(x)
plot(x.xts)
plot(x.xts['2005::2006-10'])
HTH,
David Reiner
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Erin Hodgess
Sent: Tuesday, March
On 02.03.2010 17:05, Häring, Tim (LWF) wrote:
I`m using SVMs for multi-class classification problems. Therefore I`m
using the svm() function in the package "e1071".
If I use svm(...type="C-classification") everything works fine. But
if I want to use nu-SVM with svm(..., type="nu-classifica
Try this:
plot(window(x, end = c(2006, 10)))
On Tue, Mar 2, 2010 at 2:19 AM, Erin Hodgess wrote:
> Dear R People:
>
> I have the following time series and plot:
>> x <- ts(rnorm(50),start=2005,freq=12)
>> plot(x)
>>
>
> which works fine.
>
> I would like to plot a subset of that time series, whi
On Tue, Mar 2, 2010 at 4:05 PM, Paul Hiemstra wrote:
> frederik vanhaelst wrote:
>>
>> Hi,
>>
>> How could i generate random real numbers between 0 en 2*pi?
>>
>> Thanks,
>
> Googeling for "R generate random number" gave this as a second hit (on my
> machine):
>
> http://blog.revolution-computing
hi all.
If someone have the same problems this is the answer:
-vertical legend :
> legend.krige(x.leg=c(X,X),y.leg=c(X,X),kr$pred,vert=TRUE,col="gray"(seq(.7,0,l=10)))
-sample's positions on the map:
###coords.dat=table$coords### like in
>image(kr,col="gray"(seq(.7,0,l=10)),xlim=c(-1
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of John Sorkin
> Sent: Tuesday, March 02, 2010 3:46 AM
> To: Karl Ove Hufthammer; r-h...@stat.math.ethz.ch
> Subject: Re: [R] two questions for R beginners
>
> Please take what foll
G'day Sundar,
On Tue, 2 Mar 2010 01:03:54 -0800
Sundar Dorai-Raj wrote:
> Thanks, Berwin. That works just great!
You are welcome.
I noticed by now that "cat(tmp)" is sufficient; the "tmp[1]" in
"cat(tmp[1])" was a left over from earlier attempts to get the output
to look correct.
Cheers,
Hi,
On Tue, Mar 2, 2010 at 11:05 AM, Häring, Tim (LWF)
wrote:
>
>> > I`m using SVMs for multi-class classification problems. Therefore I`m
>> using the svm() function in the package "e1071".
>> > If I use svm(...type="C-classification") everything works fine. But
>> if I want to use nu-SVM with s
I have recently upgraded to R 2.10.1 on Windows XP and am using
scripts that I've used in previous versions successfully. I'm having a
problem with choose.files. The lines read:
fura_scan_file<-choose.files(caption="Select log file (*.log) for fura-2 scans")
PI_scan_file<-choose.files(caption="Sel
I tried to shoehorn the read.* functions and match both the fixed width and
the variable width fields
in the data but it doesn't seem evident to me. (read.fwf reads fixed width
data properly but the rest
of the fields must be processed separately -- maybe insert NULL stubs in the
remaining fields a
Dear R-community,
I'm currently trying to assemble a forest plot using the "forest" function
from package "metaphor".
Works well. Even the regular "main"-argument works for adding a title to the
graph.
However, I would like to add one top row which explains the nature of the
columns. Very
Ah, I should have mentioned this. Personally I work on Macs (Leopard)
and PC's (XP Pro and XP Pro x64). Even though the PC's do have Cygwin,
I'm trying to make this code portable. So I want to avoid such things as
sed, perl, etc.
I want to do this in R, even if processing is a bit slower. Event
Hello R community,
I am hoping for some help with the following problem.
I have a data frame containing various groups. These groups are identified
by a grouping variable. I would like to add a sequential ID number to each
group to later sort these individuals within each group by this ID number.
On 02/03/2010 11:53 AM, William Dunlap wrote:
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of John Sorkin
> Sent: Tuesday, March 02, 2010 3:46 AM
> To: Karl Ove Hufthammer; r-h...@stat.math.ethz.ch
> Subject: Re: [R] two que
Try this:
data$ID <- with(data, ave(group, group, FUN = seq))
On Tue, Mar 2, 2010 at 2:53 PM, Alexander Schwall
wrote:
> Hello R community,
>
> I am hoping for some help with the following problem.
>
> I have a data frame containing various groups. These groups are identified
> by a grouping va
William,
I agree that changing syntax can lead to problems. I don't, however think
extending the language will break existing code. Providing a common syntax for
accessing matrices and dataframes will not change the way things have been done
to date, but rather how things will be done in the fut
Like this?
group<- c(1,1,1,1,2,2,2,2,2,3,3,3,3,3,3)
var2<- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15)
data<-data.frame(group, var2)
data
ddply(data,"group",transform,ID=1:length(group))
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish & Wildlife Service
California,
Dear R users,
I would like to draw map and import it in maptools/spatstat packages.
The 'raw data' I have come from a web page (...) and are
basically a list of coordinates of a polygon.
I would like to know how to import them in R; I checked the maptools
packages, but all the examples use exist
-- Forwarded message --
From: Duncan Murdoch
Date: 25 February 2010 22:30
Subject: Re: [R] Building R packages in Windows 7
To: Eric Ferreira
On 25/02/2010 7:56 PM, Eric Ferreira wrote:
> Thank you, Sir, but how can I demand it to create HTML files?
>
The tools::Rd2HTML functi
Hi all,
I have a table in database that is very long and when simplified it has only
two columns in it (id, text). id is the row, and text is the column.
Technically the text is a term and and id is the document.
If simplifying this and assuming there is only one occurrence of the term
per the doc
Liviu Andronic escribió:
On Mon, Mar 1, 2010 at 11:49 PM, Liviu Andronic wrote:
On 3/1/10, Keo Ormsby wrote:
Perhaps my biggest problem was that I couldn't (and still haven't) seen
*absolute beginners* documents.
there was once a link posted on r-sig-teaching that would prob
Thank you gentlemen,
all three solutions are working and very insightful.
Your help and time is very much appreciated.
Alexander
On Tue, Mar 2, 2010 at 1:08 PM, Felipe Carrillo wrote:
> Like this?
>
> group<- c(1,1,1,1,2,2,2,2,2,3,3,3,3,3,3)
> var2<- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15)
> da
I am trying to plot fitted lme values as a smooth line of a graph
showing the exponential relationship between temperature and soil
respiration.
In the plot, the x-axis has temperature, and the y-axis has soil
respiration. When I try to add a line showing temperature versus the
fitted va
Hello together,
Is there a tool to test the statistical differences between parameter estimates
of an nlsList fit, with more than two groups? I am able to complete the nlme
function for two groups after getting starting paramaters in nlsList, as seen
below.
>fit.nlme <- nlme(rate ~ SSmicmen(con
Try this:
DF <- read.table(textConnection("1 this
1 is
1 the
1 first
1 row
2 this
2 is
2 the
2 send
2 row"))
reshape(DF, v.names = 'V2', idvar = 'V1', timevar = 'V2', direction = 'wide')
On Tue, Mar 2, 2010 at 3:35 PM, Jan Hornych wrote:
> Hi all,
>
> I have a table in database that is very long
Hi everyone,
I am trying to find a way to calculate an adjacency matrix from a polygon
shapefile to make a CAR model in winbugs, but I would like to vary the order
of neihborhood. I used the the spdep package to create an adjacency matrix
(in list format) of order 2 and it works
Jan -
Here's one way:
tbl = data.frame(id=c(1,1,1,1,1,2,2,2,2,2),
text=c('this','is','the','first','row','this','is','the','second','row'))
xtabs(~id+text,tbl)
text
id first is row second the this
1 1 1 1 0 11
2 0 1 1 1 11
It'
Or better:
reshape(cbind(DF, value = 1), v.names = 'value', idvar = 'V1', timevar
= 'V2', direction = 'wide')
On Tue, Mar 2, 2010 at 3:49 PM, Henrique Dallazuanna wrote:
> Try this:
>
> DF <- read.table(textConnection("1 this
> 1 is
> 1 the
> 1 first
> 1 row
> 2 this
> 2 is
> 2 the
> 2 send
> 2
What I ended up using was:
cat(unclass(tmp))
--sundar
On Tue, Mar 2, 2010 at 8:58 AM, Berwin A Turlach wrote:
> G'day Sundar,
>
> On Tue, 2 Mar 2010 01:03:54 -0800
> Sundar Dorai-Raj wrote:
>
> > Thanks, Berwin. That works just great!
>
> You are welcome.
>
> I noticed by now that "cat(tmp)" i
On Mon, 01-Mar-2010 at 12:01PM -0500, Max Kuhn wrote:
|> In theory, the choice between two perfectly correlated predictors is
|> random. Therefore, the importance should be "diluted" by half.
|> However, this is implementation dependent.
|>
|> For example, run this:
|>
|> set.seed(1)
|> n <-
In most implementations of boosting, and for that matter, single tree,
the first variable wins when there are ties. In randomForest the
variables are sampled, and thus not tested in the same order from one
node to the next, thus the variables are more likely to "share the
glory".
Best,
Andy
Fro
On 2/03/2010, at 9:02 PM, Noah Silverman wrote:
> Hi,
>
> I'm just learning about poison links for the glm function.
>
> One of the data sets I'm playing with has several of the variables as
> factors (i.e. month, group, etc.)
>
> When I call the glm function with a formula that has a factor
Hi R users,
I am using the following model to analyze data from a factorial experiment
(randomized complete block design with no replication within blocks):
model <- glm(survival ~ density * vegetation + (1|block), data=sal2005,
family="binomial")
Does R use a binomial distribution in this for
Hi R-fans,
I would like put out a question to all R users on this list and hope
it will create some feedback and discussion.
1) What are your 3 most useful R package? and
2) What R package do you still miss and why do you think it would make
a useful addition?
Pulling answers together for these
On Tue, Mar 2, 2010 at 2:43 PM, Liaw, Andy wrote:
> In most implementations of boosting, and for that matter, single tree,
> the first variable wins when there are ties.
They must be in a union :-)
>> What happens if there's a third?
If they were P perfectly correlated predictors, the importanc
Hi Sebastian,
Here is an example showing a forest plot with some column headings:
library(metafor)
data(dat.bcg)
dat <- dat.bcg
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat, measure="RR")
windows(width=6.5, he
Try this; specify where you want the second one to start:
files.a <- choose.files()
# now change to the directory of the first file name to continue search
files.b <- choose.files(paste(dirname(files.a[1L]), "*", sep='/'))
On Tue, Mar 2, 2010 at 12:17 PM, Caleb Rounds wrote:
> I have recently up
Friends
Seems I've run into another snag. More of the nitty-gritty r-details I don't
understand.
So, as I mentioned below, dataset[[var_sub]] seems to be understood well by the
functions I previously used and I was able to run my loop successfully with the
[[var_sub]] as a variable-substitutio
Welcome to R/ R-PLUS Webminar Series. R-PLUS 3.3 Rocks: Interactive,
Comprehensible and Highly Visual.
http://www.xlsolutions-corp.com/webminar.asp.
March 12th @ 12PM ET (USA Time)
Increase your productivity with R-PLUS 3.3 by attending the webminar and
learning how to:
1. Interactively clicki
Hello I have 2000 univariate timeSeries of about 20 observations each, as
the following, I would like to store all of them in one object, sort of a
data frame, and to be able to recall each by its column name, which by the
way is the same as the first date. Do you know how can I do this. Thank you
1 - 100 of 143 matches
Mail list logo