Rajesh,
A quick reply to your questions concerning the diagram package:
To use different colors per arrow: just define arr.col as a
Matrix or a data.frame, and give the color number or colorname to each arrow:
To use numbers:
M <- matrix(nrow=4,ncol=4,data=0)
M[2,1]<-1 ;M[4,2]<-2;M[3,4]<-3
Hi All,
I am using fleiss kappa for inter rater agreement. Are there any know
issues with Fleiss kappa calculation in R? Even when I supply mock data
with total agreement among the raters I do not get a kappa value of 1.
instead I am getting negative values.
I am using the irr package version 0
On Mon, Jul 13, 2009 at 5:55 AM, Weiwei Shi wrote:
> Hi, there:
>
> Assume I have a dataframe with rownames like A with rownames like a to e,
>
>> A
> [,1] [,2]
> a 1 6
> b 2 7
> c 3 8
> d 4 9
> e 5 10
>
> when I use A[1,], I lost the rowname for it, like below. How co
Try A[1,,drop=FALSE] - see help("\[")
--- On Mon, 13/7/09, Weiwei Shi wrote:
> From: Weiwei Shi
> Subject: [R] how to keep row name if there is only one row selected from a
> data frame
> To: "r-h...@stat.math.ethz.ch"
> Received: Monday, 13 July, 2009, 1:55 PM
> Hi, there:
>
> Assume I hav
> In this simple example, it took less than half a second to generate the
> result. That is on a 2.93 Ghz MacBook Pro.
>
>
> So, for your data, the code would look something like this:
>
>
> system.time(DF.new <- do.call(rbind,
> lapply(split(patch_summary,
> patch_summ
Hi megh,
Unfortunately there's a known bug which prevents you from setting the
y axis title in this way. Instead you can do :
qplot(...) + scale_y_continuous("")
I think I've figured out a work around and it will be fixed in a future release.
Hadley
On Mon, Jul 13, 2009 at 2:15 AM, megh wrote
Hello Everyone
I am calculating Fleiss Kappa, I have 28 raters, 5 Subjects and 5 ratings.
The problem is that there are 2 missing values in the data.
Would it better to replace those with "0" or should those be omitted? By
omission I will be left wit only 3 subjects.
and my second problem is t
Hi, there:
Assume I have a dataframe with rownames like A with rownames like a to e,
> A
[,1] [,2]
a16
b27
c38
d49
e5 10
when I use A[1,], I lost the rowname for it, like below. How could I keep
it? Is there an easy way instead that I have to modify by mys
Here is what I do:
R> library(fBasics)
R> ts<-dummyDailySeries(x = rnorm(365), units = NULL, zone = "",
FinCenter ="")
R> class(ts)
[1] "timeSeries"
attr(,"package")
[1] "timeSeries"
R> (s <- as.character(time(ts)[1]))
[1] "1970-01-01"
R> class(s)
[1] "character"
tradenet wrote:
I added
On Jul 11, 2009, at 8:42 PM, Jonathan Greenberg wrote:
I'm a bit confused why the following command is taking an
extrodinarily long time (> 2-3 hours) to run on an 3.06ghz iMac
(brand new). I'm trying to do a stratified random sample, drawing
only a single value per UniqueID from the patc
I have not used 'sna'.
Have you tried using "debug" to walk through the code line by
line, examining and even changing things at will?
For example, how big is "rho", passed as starting values to
"optim"? If that matches the size of your adjacency matrix, it could
expose a
Hi all, whenever I try to plot a histogram using qplot() function of
"ggplot2" library, I get error like this :
> qplot(rnorm(1000), geom="histogram", binwidth=0.2, main = "", xlab="",
> ylab="")
Error in scale[[1]] : subscript out of bounds
However if I remove ylab="" argument, then it is worki
Hi Saurav!
On Sun, Jul 12, 2009 at 6:06 PM, Pathak,
Saurav wrote:
> I am new to R, I have to do a 2 step Heckman model, my selection equation is
> below which I was successful in running but I am unable to proceed further,
>
>
>
> I have so far used the following command
>
> glm(formula = s ~ age
On Jul 12, 2009, at 5:06 PM, Kayce Anderson wrote:
Hi,
I am using mgcv:gam and have developed a model with 5 smoothed
predictors
and one factor.
gam1 <- gam(log.sp~ s(Spr.precip,bs="ts") + s(Win.precip,bs="ts") +
s(
Spr.Tmin,bs="ts") + s(P.sum.Tmin,bs="ts") + s( Win.Tmax,bs="ts")
+fact
You can do something like this:
con <- file("MyFle.txt","r")
repeat {
line <- readLines(con,n=1)
if (length(line) == 0) break
ParseLine(line)
}
On Sun, Jul 12, 2009 at 6:20 PM, wrote:
> I have a long text file with uneven record length and variable structure.
> Therefore I have to read
Dear R users,
Thanks in advance.
I am using R 2.9.1 on Windows XP.
I am trying to install above, but faced problem.
I have described the steps below.
1. I like to install locally. That means I like to install after downloading
the required files into my hard drive.
2. Java Runtime E
Dear Saurav,
If you don't wish to calculate it by hand, you can extract the Inverse Mill's
Ratio using the function "invMillsRatio" from the sampleSelection package.
(http://cran.r-project.org/web/packages/sampleSelection/index.html).
You can also use the "heckit" function in the same packag
Neotropical bat risk assessments-2 wrote:
>
> Hi all,
>
> Where can I find detail of use of seewave?
> Seems no help included in the package using ??seewave does not provide
> anything.
>
help(package=seewave)
which includes information on a vignette.
Please try to remember to use an in
Hi,
My question might be a little general.
I have a number of values to select for the complexity parameters in some
classifier, e.g. the C and gamma in SVM with RBF kernel. The selection is based
on which values give the smallest cross validation error.
I wonder if the randomized splitting of
I have a long text file with uneven record length and variable structure.
Therefore I have to read it line-by-line.
I found out I can open a connection to the file and read in one line at a time.
Something like:
con <- file("MyFle.txt","r")
while (End-Of-File) {
line <- readLines(con,n=1)
Hi,
I am using mgcv:gam and have developed a model with 5 smoothed predictors
and one factor.
gam1 <- gam(log.sp~ s(Spr.precip,bs="ts") + s(Win.precip,bs="ts") + s(
Spr.Tmin,bs="ts") + s(P.sum.Tmin,bs="ts") + s( Win.Tmax,bs="ts")
+factor(site),data=dat3)
The total deviance explained = 70.4%.
On Jul 12, 2009, at 2:56 PM, David Winsemius wrote:
On Jul 12, 2009, at 3:21 PM, maram salem wrote:
Dear group,
Thank u so much 4 ur help. I've tried the link,
http://finzi.psych.upenn.edu/R/library/quantreg/html/akj.html
for adaptive kernel density estimation.
But since I'm an R beginer an
On 13/07/2009, at 8:05 AM, David Winsemius wrote:
Appears I am wrong about this. I was basing my assumption on this
interaction with the R interpreter:
?break
Error in genericForPrimitive(f) :
methods may not be defined for primitive function "break" in this
version of R
Since ``break''
On Sun, Jul 12, 2009 at 1:05 PM, David Winsemius wrote:
>
> On Jul 12, 2009, at 3:35 PM, David Winsemius wrote:
>
>>
>> On Jul 12, 2009, at 2:53 PM, Mark Knecht wrote:
>>
>>>
>>> As a test I tried to print down to the string "(all)" and then
>>> break but this code and everything I've tried so far
On Fri, 10-Jul-2009 at 09:41AM -0700, Michael wrote:
|> Here is my code:
|> mygbm<-gbm.fit(y=mytraindata[, 1], x=mytraindata[, -1],
|> interaction.depth=4, shrinkage=0.001, n.trees=2, bag.fraction=1,
|> distribution="bernoulli")
|>
|> Here is the error:
|> Error in gbm.fit(y = mytraindata[, 1
On Jul 12, 2009, at 3:35 PM, David Winsemius wrote:
On Jul 12, 2009, at 2:53 PM, Mark Knecht wrote:
As a test I tried to print down to the string "(all)" and then
break but this code and everything I've tried so far is terribly
wrong. Every attempt prints lots of error messages. I'm not gr
On Jul 12, 2009, at 3:21 PM, maram salem wrote:
Dear group,
Thank u so much 4 ur help. I've tried the link,
http://finzi.psych.upenn.edu/R/library/quantreg/html/akj.html
for adaptive kernel density estimation.
But since I'm an R beginer and the topic of adaptive estimation is
new for me, i s
On Jul 12, 2009, at 2:53 PM, Mark Knecht wrote:
Hi,
Newbie alert on for loops...
I have a bunch of data.frames built using rbind that have repeated
values in the EnTime column. I want to read the value in the EnTime
column and use it as an input to a function, but only down to the
first oc
Dear group,
Thank u so much 4 ur help. I've tried the link,
http://finzi.psych.upenn.edu/R/library/quantreg/html/akj.html
for adaptive kernel density estimation.
But since I'm an R beginer and the topic of adaptive estimation is new for me,
i still can't figure out some of the arguments of
ak
On Sun, Jul 12, 2009 at 11:53 AM, Mark Knecht wrote:
So this gets better in terms of error messages but still has problems
for(n in SystemResults$EnTime) {
if ( SystemResults$EnTime[n] == "(all)") break else X =
SystemResults$EnTime[n]
print(X)
}
> for(n in SystemResults$EnTi
Hi,
Newbie alert on for loops...
I have a bunch of data.frames built using rbind that have repeated
values in the EnTime column. I want to read the value in the EnTime
column and use it as an input to a function, but only down to the
first occurrence of the string "(all)" where I want to bre
I have no idea to be honest. I have never used the package, I simply did a
search for it. Hopefully a more experienced user can help
--- On Sun, 7/12/09, maram salem wrote:
> From: maram salem
> Subject: Re: [R] (no subject)
> To: "John Kane"
> Received: Sunday, July 12, 2009, 1:01 PM
> Dea
On Jul 12, 2009, at 1:05 PM, David Winsemius wrote:
On Jul 12, 2009, at 8:10 AM, tradenet wrote:
I added a reproducible example to my question...
ts<-dummyDailySeries(x = rnorm(365), units = NULL, zone = "",
FinCenter =
"")
"Heads up" is an English colloquial warning.
It may be repr
Hi
I have been trying so many different things to get my Inverse Mills Ratio
going for a Two stage Heckman Model, I have tried the following so far (the
commands are listed below till teh point where I get an error), I get an
error in the last sentence (marked in bold below), if this were successfu
On Jul 12, 2009, at 8:10 AM, tradenet wrote:
I added a reproducible example to my question...
ts<-dummyDailySeries(x = rnorm(365), units = NULL, zone = "",
FinCenter =
"")
"Heads up" is an English colloquial warning.
It may be reproducible on a machine that has all of your particular
You might also want to look at the doBy package - one function is summaryBy:
summaryBY(var1 + var2 ~ patient_type, data=d, FUN=summary)
david freedman
Hayes, Rachel M wrote:
>
> Hi All,
>
>
>
> I'm trying to automate a data summary using summary or describe from the
> HMisc package. I wan
On Jul 9, 2009, at 6:52 PM, Hayes, Rachel M wrote:
I'm trying to automate a data summary using summary or describe from
the
HMisc package. I want to stratify my data set by patient_type. I was
hoping to do something like:
Describe(myDataFrame ~ patient_type)
by(myDataFrame, patient_ty
On Sat, 11 Jul 2009, saurav pathak wrote:
I have so far used the following command
glm(formula = s ~ age + gender + gemedu + gemhinc + es_gdppc +
imf_pop + estbbo_m, family = binomial(link = "probit"))
My question is
1. How do i discard the non significant selection variables (one out of th
I added a reproducible example to my question...
ts<-dummyDailySeries(x = rnorm(365), units = NULL, zone = "", FinCenter =
"")
> (ts[1,0]) #returns first date in return series
GMT
1970-01-01
> ttt<-(sprintf("%s",ts[1,0]))
> print(ttt)
character(0)
> ttt<-(ts[1,0])
> print(ttt)
GMT
1970-01
Hi all,
Where can I find detail of use of seewave?
Seems no help included in the package using ??seewave does not provide
anything.
Tnx
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Hi Benoit,
Have a look at http://had.co.nz/ggplot2/stat_function.html. Does that help?
Hadley
On Tue, Jul 7, 2009 at 11:15 AM, Benoit
Boulinguiez wrote:
> Hi all,
>
> I'm smoothly transferring my lattice graphs to ggplot2 graphs, but I'm stuck
> on representing a curve from a formula.
> I'm loo
Hi Benoit,
What do you expect width to do? You are already setting the left and
right extents with xmin and xmax.
Hadley
On Thu, Jul 9, 2009 at 10:37 AM, Benoit
Boulinguiez wrote:
> Hi all,
>
> quick question: is the optional command "width" effective in the
> geom_errorbarh() layer of ggplot?
threeftmetered
Sent from my Windows Mobile® phone.
-Original Message-
From: Muenchen, Robert A (Bob)
Sent: Saturday, July 11, 2009 4:43 PM
To: R-help@r-project.org
Subject: [R] Reading data entered within an R program
Dear R-helpers,
I know of two ways to reading data within an R prog
Hi, I am trying to use the nls() function to closely approximate a vector of
values, colC and I'm running into trouble. I am not sure how if I am asking
the program to do what I think its doing, because the same minimization in
Excel's Solver does not run into problems. If anyone can tell me wha
Hi all,
I've searched threads and read up on some ways of doing this but I'm having
a hard time to get it to work. Here's my basic problem. I have the
following linear mixed model
y = Xb+Zu+e
where u~N(0,s^2*K) where K is a matrix.
I read a thread that basically suggested to decompose Zu into
45 matches
Mail list logo
