Hi
I am not sure if you can do what you want. Segments are not points so your
pch option is (I believe) ignored. You could play with lmitre and lend
parameters, but it probably would not help much.
You cold try to look at
?symbols
but you probably need to change source code to suit your needs.
Hi
My aim is to read a large .csv file into R. I ran the following code and am
using R version 10.1 on Windows.
>library(ff)
> read.csv.ffdf(x=NULL,"file.csv",fileEncoding="",nrows=-1,first.rows=NULL,next.rows=NULL,levels=NULL,appendLevels=TRUE,FUN="read.table",transFUN=NULL,asffdf_args=list(),
Dear Gregory ,
Thnaks for your reply and help. I am explaining you my problems again,
below is my script for the same .
Dom <-c (195,568,559)
fkbp <- barplot (Dom, col="black", xlab="", border = NA, space = 7,
xlim=c(0,650), ylim =c(0, 87), las = 2, horiz = TRUE)
axis (1, at = seq(0,600,10), l
Thanks Jannis!
I feel so dumb! I didn't specify the par(mar) before I wrote the
code...thanks for checking me :)
beloitstudent
On Thu, Jun 10, 2010 at 4:28 AM, jannis-2 [via R] <
ml-node+2250158-1003270618-278...@n4.nabble.com
> wrote:
>
> > I am having a few problems with this. 1) My y
> >
Thanks for the advice! I'm not sure I fully understand how to get the
text() function to rotate the axis, but I'll fiddle and see if I can't get
it to work.
On Thu, Jun 10, 2010 at 7:50 AM, Marc Schwartz wrote:
> On Jun 10, 2010, at 7:18 AM, Jim Lemon wrote:
>
> > On 06/10/2010 01:44 PM, beloit
On Fri, 11 Jun 2010, Marcin Gomulka wrote:
AFAIK a script run through source() does not have any legit way to learn
about it's own location.
You mean like 'test.R' being able to figure out that that is what it is
called when it is invoked from source()?
cat("print(eval(sys.calls()[[1]][[
Hello.
I am trying to get my axis label to read as follows
(The symbol) Delta AUC blah blah...
then below it...(some other text)
The problem is the Delta symbol shows up beside the "(some other text)"
rather than the "AUC". Does any one know how I can get the Delta to remain
beside AUC?
H
Hey Jorge,
Don't know if you received this but, essentially I found something weird
with the outputs.
Edit:
There is something funky about the code. It definitely returns the right
column of the "distance" data, but returns an incorrect row.
Code:
NCols=250
NRows=829
myMat<-matrix(runif(NCols
Hi ,
you can try using cov2cor(vcov(lm(calorie~height))) to get the correlation
matrix of estimated coefficients
boltonboy999 wrote:
>
> Hi everyone,
>
> This is pretty urgent so if anyone can help that would be great.
>
> I have a table of information. The categories are weight, height
On Fri, Jun 11, 2010 at 6:03 AM, Marcin Gomulka wrote:
> AFAIK a script run through source() does not have any legit way to learn
> about it's own location.
>
> I need this to make sure that the script will find its datafiles after I
> move the whole directory. (The datafiles are in the same dire
Edit:
There is something funky about the code. It definitely returns the right
column of the "distance" data, but returns an incorrect row.
Code:
NCols=250
NRows=829
myMat<-matrix(runif(NCols*NRows), ncol=NCols)
d<-dist(myMat)
e<-sort.list(d)
e<-e[1:5] ##Retrieve minimum 5 distances
k <-
patt <- ".*Prostate.*cvs.*"
dir(pattern = patt)
On Thu, Jun 10, 2010 at 7:07 PM, wrote:
> Actually this pattern will select the files whose name contains "csv".
> Whereas I'd like to exscude them.
> Maura
>
> -Messaggio originale-
> Da: Jorge Ivan Velez [mailto:jorgeivanve...@gmail.com]
AFAIK a script run through source() does not have any legit way to learn
about it's own location.
I need this to make sure that the script will find its datafiles after I
move the whole directory. (The datafiles are in the same directory.)
Here is a hack I invented to work around it:
print(getwd
Ted
I'm not sure I fully understand the question, but you may want to consider
creating a temporary table with a join, which you can do with a query from your
R session, and then query that table to bring the data into R. Roughly, the
logic is to leave the data in the db if you are not doing an
Actually this pattern will select the files whose name contains "csv".
Whereas I'd like to exscude them.
Maura
-Messaggio originale-
Da: Jorge Ivan Velez [mailto:jorgeivanve...@gmail.com]
Inviato: gio 10/06/2010 17.30
A: mau...@alice.it
Cc: r-h...@stat.math.ethz.ch
Oggetto: Re: [R] select
Just as a side note, perhaps this will help someone else who comes
across this topic of packages no longer being on CRAN and built before
2.11.0 in that case they can't install the package...
I use package 'sma' which was last maintained sometime in the mid 2000's
and obviously it's not on CRA
Philip wrote:
Dear list,
Is it possible to replace the median in the following code so as to calculate 5
and 95% tiles?
with( raunheim, aggregate(no, by = list(wd,hr), FUN = "median", na.rm=TRUE) )
Please give reproducible examples that we can run by cutting and pasting!
You should just g
On Thu, Jun 10, 2010 at 10:37 AM, Patrick Burns
wrote:
> See 'The R Inferno' Circle 2 for why
> this takes so long, and what to do
> about it.
>
Also check this Rnews issue [1], page 46, on how to make loops faster.
Liviu
[1] http://www.r-project.org/doc/Rnews/Rnews_2008-1.pdf
__
Hi,
You could reorder the factor levels before plotting,
x$n = factor(x$n, levels=c("va","vp", letters[1:3]))
last_plot() %+% x
or you could avoid using factors in the first place,
x <- data.frame(cbind(n,p,pm,pn), stringsAsFactors=FALSE)
last_plot() %+% x
HTH,
baptiste
On Thu, Jun 10, 20
Thanks,
Actually, thanks to the info Henrique sent, I have made decent progress.
In actuality, I could have just submitted a SELECT * on the second table,
which would give me everything, just like Henrique's suggestion, and yours,
would give me. The problem is that that table is HUGE (I don't wa
Here's one possibility:
allfiles = list.files()
csvfiles = list.files(pattern='csv')
allfiles[!allfiles %in% csvfiles]
- Phil Spector
Statistical Computing Facility
Departmen
Hi Dr. Lumley,
Sorry about the succession of the 2 messages. It wasn't intended. This is my
first post to the forum and I first signed up to Nabble and submitted the
post and later got the message from Nabble that the post has not been
accepted by the R forum and that I need separate subscription
Hello Felipe,
Is this what you want?
format(as.Date("3/10/10", format="%m/%d/%y"), "%B %d, %Y")
Josh
On Thu, Jun 10, 2010 at 8:29 AM, Felipe Carrillo
wrote:
> Hi:
> Can't find a way to convert from shortDate to LongDate format. I got:
> 3/10/10 that I want to convert to March 10, 2010. I am us
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of baptiste auguie
> Sent: Thursday, June 10, 2010 11:50 AM
> To: Kim Jung Hwa
> Cc: r-help@r-project.org; ggplot2
> Subject: Re: [R] Order labels in qplot() - ggplot2 {help}
>
> Hi
Dear list,
Is it possible to replace the median in the following code so as to calculate 5
and 95% tiles?
with( raunheim, aggregate(no, by = list(wd,hr), FUN = "median", na.rm=TRUE) )
Regards,
Philip Sinclair
__
R-help@r-project.org mailing list
ht
I don't think that this would be considered a bug. The
reason for the discrepancy between use="complete.obs"
and use="pairwise.complete.obs" for the case of the
Spearman correlation of two vectors x, y is this:
"pairwise" does complete.cases(x,y) and then ranks;
this is also what's done in cor.te
Hello,
I want to arrage the label according to my preference eg.. (va, vp, a, b,
c) but don't know how to supress default ordering. Any
suggestions?
Please try the code below:
n <- c("va", "vp", "a", "b", "c")
p <- c(2, 2,1, 3,5)
pm<- c(3,4,2,5,4)
pn <- c(1,1,1,2,3)
x<-data.frame(cbind(n,p,pm,pn
On 2010-06-10 8:27, Graves, Gregory wrote:
What am I failing to understand here?
Several things; see below.
The script below works fine if the dataset being used is
DNase1<- DNase[ DNase$Run == 1, ] per the example given in
help(nlrob).
Obviously, I am trying to understand how to use nl
On Thu, Jun 10, 2010 at 9:30 PM, wrote:
> Hello list,
>
> Is there a way to add contour lines to a levelplot at different breakpoints
> than are used for the colors? For example:
>
>
> library(lattice)
>
> # colors good but too many contours
> levelplot(volcano, at=94:195, contour=TRUE)
>
> # I t
Your question is not really clear, do either of these examples do what you want?
with(anscombe, plot(x1, y2, ylim=range(y2,y3)) )
with(anscombe, points(x1, y3, col='blue', pch=2) )
with(anscombe, segments(x1, y2, x1, y3, col=ifelse( y2>y3, 'green','red') ) )
with(anscombe, plot(x1, y2, ylim=
> newmat <- cbind( oldmat, order=seq(nrow(oldmat)) )
See ?seq and/or ?rev for descending.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r
On Jun 10, 2010, at 12:19 PM, Marc Schwartz wrote:
> Felipe,
>
> I would not do the processing in TeX, but do it in R and then pass the
> results to the \Sexpr{}'s.
>
> If I am correctly understanding the process flow, put the following R code
> chunk before the point where you need to output
Marc:
Thanks for reinforcing that, I was just trying to go that route.
It seems to be simpler to import the dataset and just grab the first and
last date from it and then format it. Thanks again.
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish & Wildlife Serv
Felipe,
I would not do the processing in TeX, but do it in R and then pass the results
to the \Sexpr{}'s.
If I am correctly understanding the process flow, put the following R code
chunk before the point where you need to output the formatted dates:
<>
START <- format(as.Date(report[1, 1],
correction.
See ?daisy in the PACKAGE cluster. *slaps head*
cheers
Joris
On Thu, Jun 10, 2010 at 7:02 PM, Joris Meys wrote:
> see ?daisy in the library cluster
>
> Cheers
> Joris
>
> On Thu, Jun 10, 2010 at 6:12 PM, kapil mahant wrote:
>> All,
>>
>> How can we find a distance matrix for categori
see ?daisy in the library cluster
Cheers
Joris
On Thu, Jun 10, 2010 at 6:12 PM, kapil mahant wrote:
> All,
>
> How can we find a distance matrix for categorical data
>
> ie. given a csv below
>
> var1 var2 var3 var4
> element1-1 yes x a
The lm and aov functions can take a matrix response allowing you to
fit all of the responses for a single attribute simultaneously.
On Thu, Jun 10, 2010 at 8:47 AM, melissa wrote:
> Dear all R users,
> I want to realize 800 000 ANOVAS and to store Sum of Squares of the effects.
> Here is an ext
Marc:
My report is done every two weeks and is created automatically.
I click a command button on an Excel form and it runs a .rnw script
in R creating a latex dynamic report. Excel sends 15 days of data
to R, eg: 6/1/10 to 6/15/10. Right above my report I usually write the range
of the report m
I'm doing some coxph() analyses with a large and complex dataset. The
data was collected in different centers, so I am using
strata(centers) to stratify the analysis.
My main issue is, not all centers collected all the variables, so for
a model such as:
coxph(Surv(days, cancer) ~ varA + sex +
Hi everyone,
I'm doing some coxph() analyses with a large and complex dataset. The data was
collected in different centers, so I am using strata(centers) to stratify the
analysis.
My main issue is, not all centers collected all the variables, so for a model
such as:
coxph(Surv(days, cancer)
You don't say what you mean by validation and calibration of the model (these
terms have multiple meanings) but if you mean looking at how good the
prediction is, all the raw materials are available through svykm() and the
predict() method for svycoxph(). There aren't any predefined function
All,
How can we find a distance matrix for categorical data
ie. given a csv below
var1 var2var3var4
element1-1 yesx a k
element1-2 no y b l
element1-3 maybe y c m
how can i
Please let me know where i have to type below thing to RSQLite database get
installed.Please let me know the solution.Thanks in advance
RSQLite -- Embedding the SQLite engine in R
(The RSQLite package includes a recent copy of the SQLite
distribution from http://www.sqlite.org.)
Installat
Dear all R users,
I want to realize 800 000 ANOVAS and to store Sum of Squares of the effects.
Here is an extract of my table data
Product attribute subject rep t1 t2 t3 ⦠t101
P1 A1 S1 R1 1 0 0 ⦠1
I want to realize 1 ANOVA per timepoint and per attribute, there are 101
timepoints and 8 attr
Its best to use dput when displaying your data in r-help as in dput(z)
or dput(head(z)) if large.
Try this:
> library(zoo)
> # test data
> z <- zoo(1:2, factor(c("04/16/09","04/17/09")))
>
> # 1
> aggregate(z, function(x) as.Date(x, "%m/%d/%y"))
2009-04-16 2009-04-17
1 2
>
> # o
It gets my vote.
-Peter Ehlers
On 2010-06-10 9:30, Bert Gunter wrote:
Is the quote below from Deepayan a possible Fortune? It put a sardonic smile
on my lips.
Bert Gunter
Genentech Nonclinical Biostatistics
... You should really get into the habit of reading
documentation; it can often tu
Hi
r-help-boun...@r-project.org napsal dne 10.06.2010 15:56:06:
> Hello everyone,
>
> I have a matrix of over 4 line and about 30 columns.
Matrix or data.frame?
>
> For my analysis I would like to add another column with ascending
numbers
> (column header should be "order", and than 1,2,
Hello list,
Is there a way to add contour lines to a levelplot at different
breakpoints than are used for the colors? For example:
library(lattice)
# colors good but too many contours
levelplot(volcano, at=94:195, contour=TRUE)
# I thought something like this might work
levelplot(volcano,
Hello, I am trying to run some examples from the book of P.Congdon. If I run
the following script
# Program 7.2 Bayesian Statistical Modelling - Peter Congdon
#
library(R2WinBUGS)
setwd("c:/temp/R")
mo <- function() {
rho ~ dbeta(1,1)
th ~ dgamma(0.001,0.001)
Y[1] ~ dpois(th)
for (t in 2:
On Jun 10, 2010, at 10:21 AM, Felipe Carrillo wrote:
> Hi:
> Can't find a way to convert from shortDate to LongDate format. I got:
> 3/10/10 that I want to convert to March 10, 2010. I am using:
>
> \documentclass[11pt]{article}
> \usepackage{longtable,verbatim}
> \usepackage{ctable}
> \usepackag
Hi Maura,
Try list.files(pattern = 'csv')
HTH,
Jorge
On Thu, Jun 10, 2010 at 11:24 AM, <> wrote:
> I have the following files list:
>
> > list.files()
> [1] "Prostate-Cancer_cvs_Dir"
> [2] "Prostate_Cancer-miRNAs&Genes.Pathway.xml"
> [3] "Prostate_Cancer_Pathways-miRNAs-GeneTargets-Dir"
>
Is the quote below from Deepayan a possible Fortune? It put a sardonic smile
on my lips.
Bert Gunter
Genentech Nonclinical Biostatistics
... You should really get into the habit of reading
documentation; it can often turn out to be very useful.
-Deepayan Sarkar
_
Hi:
Can't find a way to convert from shortDate to LongDate format. I got:
3/10/10 that I want to convert to March 10, 2010. I am using:
\documentclass[11pt]{article}
\usepackage{longtable,verbatim}
\usepackage{ctable}
\usepackage{datetime}
\title{my title}
\begin{document}
% Convert date
\dd
I have the following files list:
> list.files()
[1] "Prostate-Cancer_cvs_Dir"
[2] "Prostate_Cancer-miRNAs&Genes.Pathway.xml"
[3] "Prostate_Cancer_Pathways-miRNAs-GeneTargets-Dir"
[4] "Prostate_Cancer_Pathways-miRNAs-GeneTargets-Dir.zip"
[5] "Prostate
Hi:
Can't find a way to convert from shortDate to LongDate format. I got:
3/10/10 that I want to convert to March 10, 2010. I am using:
\documentclass[11pt]{article}
\usepackage{longtable,verbatim}
\usepackage{ctable}
\usepackage{datetime}
\title{my title}
\begin{document}
% Convert date
\dd
> prop.test(c(271,433),c(6164113,5572041))
2-sample test for equality of proportions with continuity correction
data: c(271, 433) out of c(6164113, 5572041)
X-squared = 54.999, df = 1, p-value = 1.206e-13
alternative hypothesis: two.sided
95 percent confidence interval:
-4.291428e-05 -
The OP asked for a smooth curve,
foo = splinefun(x,y)
curve(foo, min(x), max(x))
# points(x,y)
I'm sure a R wizard could make it a one-liner.
HTH,
baptiste
On 10 June 2010 16:48, Mario Valle wrote:
> x<-c(1:6)
> y<-c(.01,.09,.08,.03,.001,.02)
> plot(x,y, type='l')
>
> Please try ?plot before
x<-c(1:6)
y<-c(.01,.09,.08,.03,.001,.02)
plot(x,y, type='l')
Please try ?plot before asking to the list.
Ciao!
mario
On 10-Jun-10 12:11, suman dhara wrote:
Sir,
I have a problem regarding drawing curve.I pose the problem as follows:
suppose I have two vectors:
x<-c(1:6)
y<-c(.0
or
> or # 3
> aggregate(z, as.Date(time(z), "%m/%d/%y"))
2009-04-16 2009-04-17
1 2
On Thu, Jun 10, 2010 at 9:58 AM, Gabor Grothendieck
wrote:
> Its best to use dput when displaying your data in r-help as in dput(z)
> or dput(head(z)) if large.
>
> Try this:
>
>> library(zoo)
>
Hi:
I downloaded the patch, how do I incorporate it to my current version of
ggplot2?
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish & Wildlife Service
California, USA
- Original Message
> From: Hadley Wickham
> To: Karsten Loesing
> Cc: "r-help
What am I failing to understand here?
The script below works fine if the dataset being used is
DNase1 <- DNase[ DNase$Run == 1, ] per the example given in
help(nlrob).
Obviously, I am trying to understand how to use nls and nlrob to fit
curves to data using R.
#package=DAAG
attac
Hello everyone,
I have a matrix of over 4 line and about 30 columns.
For my analysis I would like to add another column with ascending numbers
(column header should be "order", and than 1,2,3,4 the end of the
matrix).
During my analysis I reorder them ( due to merge commands by a differen
Hello, please excuse me if this is a repost, I think I got confused
about how to reply. Anyway, I posted the message below about a day and a
half ago through Nabble but have not been cleared until today. Here it
goes again. Please bear with me I will get better with time.
Start of message s
Hi,
How to mark points on x axis of a graph keeping x axis as constant and
changing y from y1 to y2 respectively. I want to highlight the area from y1
to y2.
Any suggestions
Thank you
Jeet
[[alternative HTML version deleted]]
__
R-help@r-pro
Thanks to both Prof. Ripley and Duncan Murdoch for correctly diagnosing my
problem as failure to redirect stderr (nothing to do with buffering).
Both pointed me to R CMD BATCH which would have done the job, if I hadn't
wanted to run interactively - in my real task (not the example below) I had
Dear R People:
I have a zoo object with its date index as a factor.
> xAle1.zoo$index
Error in xAle1.zoo$index : $ operator is invalid for atomic vectors
> str(xAle1.zoo)
Class 'zoo' atomic [1:32] 1253 1316 1038 1157 1710 1489 1159 1142 945 1245 ...
..- attr(*, "index")= Factor w/ 32 levels "0
Hi Florian,
Have you tried to replace each '\n' with '\r\n'. That did the trick for
me.
HTH,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek
team Biometrie & Kwaliteitszorg
Gaverstraat 4
9500 Geraa
Have a look at the ggplot2 package
x <- 1:6
y <- c(.01,.09,.08,.03,.001,.02)
dataset <- data.frame(x, y)
library(ggplot2)
ggplot(data = dataset, aes(x = x , y = y)) + geom_smooth()
ggplot(data = dataset, aes(x = x , y = y)) + geom_smooth() +
geom_point()
ggplot(data = dataset, aes(x = x , y = y))
Hi Jim and Hi Jannis,
Thanks very much to both of you for your help! Both methods work perfectly!
Always good to know that there is more than one way to skin a cat when it
comes to R! I will just need to get a grip on the regular expressions, it
would seem.
Many thanks again for you r help,
muc
Sir,
I have a problem regarding drawing curve.I pose the problem as follows:
suppose I have two vectors:
x<-c(1:6)
y<-c(.01,.09,.08,.03,.001,.02)
plot(x,y)
It gives me the plotted points.But I want to draw a smooth curve passing
througt these points.
How can I do this?
Thanks & Regards,
Suman D
Hi Bill,
No worries, always a million things to do! Thanks very much for the reply,
that has cleared that up and I'll look out for the update next week.
Many thanks,
Ross
--
View this message in context:
http://r.789695.n4.nabble.com/Intra-Class-correlation-psych-package-missing-data-tp17739
Hi,
I'm totally new to R so I apologise for the basic request. I am looking at
the incidence of a disease over two time periods 1990-1995 and 2003-2008. I
have counts for each year, subdivided into three disease categories and by
males/females.
I understand that I need to analyse the data using p
Hi,
I have trouble with Sweave (I think) cutting of my newlines.
As stated in the help of Sweave, I generate tex code straight from R for
dynamically computed reports.
If I do this in R:
for (i in 0:4) {cat("\n",i,"\n")};cat("\n 3")
0
1
2
3
4
3
The output looks correct.
However
Try this also:
rowsum(`mode<-`(!is.na(m[,-(1:2)]), 'numeric'), m[,2])
On Wed, Jun 9, 2010 at 10:03 PM, steven mosher wrote:
> # create a matrix with some random NAs in it
> > m<-matrix(NA,nrow=15,ncol=14)
> > m[,3:14]<-52
> > m[13,9]<-NA
> > m[4:7,8]<-NA
> > m[1:2,5]<-NA
> > m[,2]<-rep(1800:1804
Cool! Thanks Karsten. If you send me a github pull request I'll incorporate it.
Hadley
On Thursday, June 10, 2010, Karsten Loesing wrote:
> Hi everyone,
>
> here's the same patch as a new branch on GitHub.
>
> http://github.com/kloesing/ggplot2/commit/a25e4fbfa4017ed1
>
> Best,
> --Karsten
>
>
On Jun 10, 2010, at 7:18 AM, Jim Lemon wrote:
> On 06/10/2010 01:44 PM, beloitstudent wrote:
>>
>> Hello experts.
>>
>> Sorry this is such a novice question, but I have been trying, fruitlessly to
>> get my x axis to angle at 45 degrees. I have tried the text() call with srt
>> and adj and just
I think you can do this:
ids <- dbGetQuery(conn, "SELECT id FROM my_table")
other_table <- dbGetQuery(conn, sprintf("SELECT * FROM my_other_table WHERE
t1_id in (%s)", paste(ids, collapse = ",")))
On Wed, Jun 9, 2010 at 11:24 PM, Ted Byers wrote:
> I have not found anything about this except th
This is a typical problem we give as a homework to students. If you
can't solve this yourself, you really need to brush up your
statistical knowledge or look for a statistician close by to cooperate
with.
Take a look at these :
http://www.ats.ucla.edu/stat/r/seminars/repeated_measures/repeated_me
On Thu, 10 Jun 2010, Keith Jewell wrote:
In MS Windows I
a) invoke Rterm from a batch file (test.bat)
b) to execute commands from a script (m:\test.rsc)
c) capturing output in a log file (m:\test.log)
BUT if the script results in an error the error message is NOT written to
the log file, lea
On 6/10/2010 8:26 AM, Samuel Okoye wrote:
> Hello,
>
> Is there any R function which does power calculation for unbalanced groups
> (n1 neq n2)? Since power.t.test has n
>
> Number of observations (per group).
>
> Many thanks,
> Samuel
See pwr.t2n.test() in the pwr package.
http://finzi.psy
On Thu, Jun 10, 2010 at 2:04 PM, Gabor Grothendieck
wrote:
> On Thu, Jun 10, 2010 at 5:33 AM, Joris Meys wrote:
>> This is only valid in case your X matrix is exactly the same, thus
>
> The poster did not give a full explanation so the best that can be
> done without getting into an extended ques
Hello,
Is there any R function which does power calculation for unbalanced groups (n1
neq n2)? Since power.t.test has n
Number of observations (per group).
Many thanks,
Samuel
[[alternative HTML version deleted]]
__
R-help@r-project
On 06/10/2010 01:44 PM, beloitstudent wrote:
Hello experts.
Sorry this is such a novice question, but I have been trying, fruitlessly to
get my x axis to angle at 45 degrees. I have tried the text() call with srt
and adj and just cannot get my labels to tilt. Does anyone have any other
sugges
On Thu, Jun 10, 2010 at 5:33 AM, Joris Meys wrote:
> This is only valid in case your X matrix is exactly the same, thus
The poster did not give a full explanation so the best that can be
done without getting into an extended question and answer is to make
some assumptions, show the result and hop
On 06/09/2010 10:14 PM, Mustafa Tufan Turp wrote:
Dear R users,
I need to learn plotting Taylor diagram. I dowloaded the "plotrix" that
includes "taylor.diagram" but I dont know how to study with real datasets.
Anyone help me?
Hi Tufan,
Try this:
www-pcmdi.llnl.gov/about/staff/Taylor/CV/Tayl
On Mon, 2010-06-07 at 10:25 -0700, Dipa Hari wrote:
>
> Hello Sir,
> I am using mgcv package for my data.
> My model isy~x1+f(x2),I want to find out the function f(x2).
> Following is the code.
>
> sm1=gam(y~x1+s(x2),family=binomial, f)
> summary(sm1)
> plot(sm1,residuals=TRUE, xlab="AGE",pch
Thank you very much for your answers.
The setup is not as you described - all participants have 4 stimuli - 2 for
punishment and 2 for reward. all participants see them in mixed blocks.
I try to measure the correlation between some personality traits (measured
with personality questionnaires) and t
First of all, read the help file ?data.frame. What you do is adding a
variable row.names, column.names and title to your dataframe. What you
want to do is set the rownames and colnames.
> results <- data.frame(matrix(c(1,2,3,4),nrow=2,ncol=2))
> rownames(results) <- c("a","b")
> colnames(results
Dear Chris,
Just define the following function:
transf.iasin <- function(x) {
z <- sin(x)^2
return(z)
}
to give the inverse of the arcsine transformation. Then specify this function
under the atransf argument.
I will add the transf.iasin() function to the package.
Best,
--
Wolfgang Vie
As I explained, you cannot just test two models with different
dependent variables. You can model the difference you calculated, if
you think you can give a sensible interpretation to it. But be aware
of the fact that your model will tell you something about the relation
between your predictors and
Hi:
I had Harold's idea (matrix indexing), but I was curious to see which of
these ran fastest. I simulated
1000 rows and three columns of binary data, along with a fourth column that
sampled the values 1:3
1000 times. Here are the timings:
> f <- as.data.frame(matrix(rbinom(3000, 1, 0.4), nrow =
Keith Jewell wrote:
In MS Windows I
a) invoke Rterm from a batch file (test.bat)
b) to execute commands from a script (m:\test.rsc)
c) capturing output in a log file (m:\test.log)
BUT if the script results in an error the error message is NOT written to
the log file, leaving me problems wh
markw...@afrihost.co.za wrote:
Package arulesSequences isn't on CRAN for automatic package install.
I downloaded a *.tar.qz version because no *.zip for Windows offered.
Why is this?
I expanded *.tar.qz in ~R/R-2.11.0/library/arulesSequences
I then assumed that R would list it in the packages lis
On 10.06.2010 12:28, markw...@afrihost.co.za wrote:
Package arulesSequences isn't on CRAN for automatic package install.
I downloaded a *.tar.qz version because no *.zip for Windows offered.
Why is this?
Because the DESCRIPTION file and hence the web page of the package at
http://cran.r-proje
Hi all,
I'm having great trouble working with the Cforest (from the party package)
and Random forest functions. Large data set seem to create very large model
objects which means I cannot work with the number of observations I need to,
despite running on a large 8GB 64-bit box. I would like the
In MS Windows I
a) invoke Rterm from a batch file (test.bat)
b) to execute commands from a script (m:\test.rsc)
c) capturing output in a log file (m:\test.log)
BUT if the script results in an error the error message is NOT written to
the log file, leaving me problems when the error is from a
On 10.06.2010 11:15, Petr PIKAL wrote:
Hi
Uwe Ligges napsal dne 10.06.2010
10:37:05:
On 10.06.2010 10:19, Petr PIKAL wrote:
Hi
Peter Ehlers napsal dne 09.06.2010 19:05:24:
Soapbox:
Well, if you're just starting out with R it would be
a VERY good idea to learn right away that T is no
Package arulesSequences isn't on CRAN for automatic package install.
I downloaded a *.tar.qz version because no *.zip for Windows offered.
Why is this?
I expanded *.tar.qz in ~R/R-2.11.0/library/arulesSequences
I then assumed that R would list it in the packages list but it doesn't.
What am I doing
Hello,
Or ?aggregate
stdate<-c(rep(1,5),rep(2,5))
domaindesc<-c(1,1,1,2,2,2,3,3,3,3)
logins<-sample(1:10, 10)
dummy.data<-as.data.frame(cbind(stdate,domaindesc,logins))
aggregate(dummy.data[,3], by=list(dummy.data[,1],dummy.data[,2]), FUN=sum)
# same result as tapply, but with only levels that
I'll try to add some more information regarding my experiment - maybe that
would help clear things out.
Instead of actually measuring the learning curve (i.e. number of correct
responses per block) I created a variable that substract the number of
correct answers from the last block with that of th
1 - 100 of 127 matches
Mail list logo
