hi everyone,
i have a problem with a function that i defined:
the function needs to use a function which is defined outside the function
but
i realised that this is not working. a friend told me that this must be a
problem with
hidden parameters.
a workaround works when i just define all
It rather depends on how you have your data stored. Here is one possibility
you might want to look at:
> con <- textConnection("
+ id_name x1 x2 x3
+ aa101 1,4,5 2 1
+ aa102 1,2,5 1 2
+ aa103 1,2,5 1 1
+ aa104 1,2,3 1 2
+ aa105 1,5 2 2
+ aa106 1,2,5 2 2
+ aa107 1,2,5 2 1
Dale Steele wrote:
I would like to randomly shuffle a distance object, such as the one
created by ade4{dist.binary} below. My first attempt, using
sample(jc.dist) creates a shuffled vector, losing the lower triangular
structure of the distance object. How can I Ishuffle the lower
triangular part
Does anybody know how to solve this error?
postscript(file= "file.ps", family = "Helvetica", font = "Helvetica")
plot.igraph(g, layout=layout.circle, vertex.label.font=2)
dev.off()
error in text.default(x, y, labels = labels, col = label.color, family
= label.family, :
family 'serif' not inc
> i have a problem with a function that i defined:
>
> the function needs to use a function which is defined outside the
function
> but
> i realised that this is not working. a friend told me that this must be
a
> problem with
> hidden parameters.
>
> a workaround works when i just define
Hello,
I'm trying to use adabag to make bagging and boosting with bagging() and
adabost.M1(), respectively, but in both cases it produces an abnormal
termination of R.
My code is:
bagging(I.NOSOCO~EDAD+SEXO+ESTANCIA+ADMISIÓN+T.CIRUGÍ+DURACIÓN+CONTAMIN
+PROFILAX+E.PREOPE+V.PERIFE+V.CENTRA+S.VES
Knut, see this thread here:
http://lists.gnu.org/archive/html/igraph-help/2007-07/msg00010.html
Best,
Gabor
On Thu, Apr 16, 2009 at 10:03 AM, Knut Krueger wrote:
> Does anybody know how to solve this error?
>
> postscript(file= "file.ps", family = "Helvetica", font = "Helvetica")
> plot.igraph(g
Dear all, I have a database x,y,value imported in R with read.table:
dati<-
read.table("dati.dat")
value is a categorical data (land use) and i want to
plot in the same colour the same land use. It is possible with R. Thanks a lot
__
R-help@r-project.
At 11:10 16.04.2009, giuseppef...@libero.it wrote:
Dear all, I have a database x,y,value imported in R with read.table:
dati<-
read.table("dati.dat")
value is a categorical data (land use) and i want to
plot in the same colour the same land use. It is possible with R. Thanks a lot
Dear R users,
Is there an R function to compute the pooled mean?
Many Thanks
Eleni Rapsomaniki
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.h
Dear karin,
Try rescaler from the reshape package:
install.packages("reshape")
library("reshape")
newDF <- apply(oldDF, 2, rescaler, type = "sd")
HTH,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderz
On Wed, 15 Apr 2009, j.k wrote:
Hello alltogheter,
I have a little problem regarding merging to zoo series.
I want to merge two zoo series to reduce the timegaps between the stamps.
I use the following code:
data.test <-
as.POSIXct(seq(data.input01[1,1],data.input01[nrow(data.input01),1],900),
David Hajage writes:
> Hi Karin,
>
> I'm not sure I understand... Is this what you want ?
>
> d$y - mean(d$y)/sd(d$y)
Yes, and also a bit no.
Each column in my data frame represents one data set. For every
element in this data set I want to know the z value for that
element. I.e: I want to c
Hi,
I need to draw a line segment between two points on different plots in the
same multigraph.I've tried looking at the zoominplot function in plotrix but
havent understood much.any help is appreciated
~Aks
[[alternative HTML version deleted]]
__
Hi,
You should give a minimal reproducible example so that we know more
precisely what you want to do (what's a "multigraph"?).
Perhaps you can get inspiration from Paul Murrell's R graphics book,
in particular Figure 5.22,
http://www.stat.auckland.ac.nz/~paul/RGraphics/customgrid-lattice
As Michael notes, scale is what you want. Also, your request has an
incorrect definition of z scores:
d$y - mean(d$y)/sd(d$y) #incorrect
(d$y - mean(d$y) ) / sd(d$y) #correct
On Thu, Apr 16, 2009 at 7:40 AM, Michael Conklin
wrote:
> newDF<-as.data.frame(scale(oldDF))
>
> see ?scale
>
> Hope tha
newDF<-as.data.frame(scale(oldDF))
see ?scale
Hope that helps.
Michael Conklin
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Karin Lagesen
Sent: Thursday, April 16, 2009 5:29 AM
To: r-help@r-project.org
Subject: Re: [R] perfor
Here's an example adapted from ?grid.arrows
pushViewport(viewport(layout=grid.layout(1, 2)))
pushViewport(viewport(layout.pos.col=1,
layout.pos.row=1))
grid.rect(gp=gpar(col="grey"))
grid.points()
grid.points(x=0.5, y=0.5, gp=gpar(col="red"))
popViewport()
pushViewport(vie
Dieter,
You could always try the "Classification, clustering, and phylogeny
estimation" list which often includes posts regarding MDS:
http://lists.sunysb.edu/index.cgi?A0=CLASS-L
regards
Bob
__
R-help@r-project.org mailing list
https://stat.eth
Hi there,
I'm travelling right now so I can't really check this but it seems that
the problem is that cluster.stats needs a partition as input. hclust
doesn't give you a partition but you can generate one from it using
cutree.
BTW, rather use "<-" than "=".
Best wishes,
Christian
On Wed, 1
You can look in MASS 4 for this formula on page 308 . Go to the
source and ask the horse he'll give you an answer that you endorse.
On Thu, Apr 16, 2009 at 8:13 AM, Bob Green wrote:
>
> Dieter,
>
> You could always try the "Classification, clustering, and phylogeny
> estimation" list which ofte
dear R experts:
does R have "static" variables that are local to functions? I know
that they are usually better avoided (although they are better than
globals).
However, I would like to have a function print how often it was
invoked when it is invoked, or at least print its name only once to
STD
On 4/16/2009 8:46 AM, ivo welch wrote:
dear R experts:
does R have "static" variables that are local to functions? I know
that they are usually better avoided (although they are better than
globals).
You put such things in the environment of the function. If you don't
want them to be global
You can use 'local'
> counter()
> # create a function with a 'local static' variable
> counter <- local({
+ i <- 0
+ function() {
+ i <<- i + 1
+ i
+ }
+ })
>
> counter()
[1] 1
> counter()
[1] 2
> i <- 42 # change 'i' at this level
> counter() # has not effect on the
Hello
I have a data frame looking like this;
> head(SNP2)
AnimalMarker All1 All2 x
3213 194073197 P1001222
1295 194073197 P1002121
0915 194073197 P1004222
2833 194073197 P1005110
1487 194073197 P1006222
1885 19407
Worked perfect...
Thanks a lot!
Johannes
j.k wrote:
>
> Hello alltogheter,
> I have a little problem regarding merging to zoo series.
> I want to merge two zoo series to reduce the timegaps between the stamps.
> I use the following code:
>
> data.test <-
> as.POSIXct(seq(data.input01[1,1],dat
Hi,
I would like to draw a graph as follows:
A simplified example is that on the X axis are different countries, I have
several temperature measurements taken from each country and would like to plot
these linearly above each country. So one would imagine that cold countries
would have lots o
Is there a way to change the default settings for Sweave figure output?
Sweave produces an .eps and a .pdf file as output (7x7 inch). Changing pdf
defaults changes the pdf() function output size but does not affect the pdf
size produced by Sweave.
pdf.options(height=3)
pdf("myfile.pdf")
p
Hi,
How can I find the p-value for the F test for the interaction terms in a
regression linear model lm ?
I appreciate your help
--
View this message in context:
http://www.nabble.com/F-test-tp23078122p23078122.html
Sent from the R help mailing list archive at Nabble.com.
Hi R,
I have explored R archives a lot but couldn't find an efficient way of
doing the following:
I want to split a vector into sets of equal sizes. Is there any inbuilt
function of doing so with the option of specifying how to treat the
remaining observations. For example: suppose I want to
Hi
Using R 2.8.1. I have list object called "AuxData". Inside a browser(),
get("AuxData")
succeeds, while
getAnywhere("AuxData")
fails with the error "no object named âAuxDataâ was found". I'm
curious to know if this could be a bug. If yes, I'll try to come up
with a reproducible examp
Ok, I'll try to explain my issue. I have a monthly series (CPI index) and I
want to interpolate it using a specific lagged harmonized formula to get the
corresponding daily series. The formula is the following
CPI^=CPI(t-3)+(d-1)/D*(CPI(t-2)-CPI(t-3))
where
CPI^ is the CPI for the day we are ca
Sorry if I am cluttering up this list with too many rank beginner
questions. But I *am* a rank beginner, and I am at my wit's end with
this one:
I have a data frame df, with components x,y,z. This command to make a
contour plot with lattice graphics works just fine:
> contourplot(z~x*y,data=df,
Hi,
Using R 2.8.0 on Windows,
I have a mathematical model written in C++, the model writes to file (.txt) a
set of numbers I want to plot in R.
The model iterates over about 10,000 runs, each time overwriting the old with
the new set of results to the output file.
I want to be able to continu
On 4/16/2009 9:52 AM, Benjamin Tyner wrote:
Hi
Using R 2.8.1. I have list object called "AuxData". Inside a browser(),
get("AuxData")
succeeds, while
getAnywhere("AuxData")
fails with the error "no object named âAuxDataâ was found". I'm
curious to know if this could be a bug. If yes,
Liz Webb schrieb:
Hi,
I would like to draw a graph as follows:
A simplified example is that on the X axis are different countries, I have
several temperature measurements taken from each country and would like to plot
these linearly above each country. So one would imagine that cold countries
Call off the dogs. The problem was that because of the scale of my data
(<<1) contourplot somehow can't handle labeling the contours if it
"paints" the plot region as well. I can easily get around this by
re-scaling the data. But, if anyone can explain to me why throwing the
region=T flag makes
Hi,
as a R-newcomer I would like to create some new data frames from a given
data frame. The first new data frame should content all pairs of the
columns of the original data frame. The second new data frame should
content all tripels of of the columns of the original data frame and the
last the q
I want write article by russian language using Sweave. For cyrillic text
LaTeX use T2A encoding
\usepackage[T2A]{fontenc}
But in Sweave.sty we find:
\RequirePackage[T1]{fontenc}
It is source of critical problem.
For example Rnw file
$ cat estimation.Rnw
\documentclass[A4paper]{article}
\us
good morning
This question is not a stats question per say but a data management and
lattice plotting problem. I apologize now if I'm asking an inappropriate
question to this gracious group.
I'm need to bring in approximately 100 *.dbf files into R but I'm having
difficultly understanding sever
I have csv files imported in r each with 2 columns and many many rows. I have
sorted the data in them but want to extract some values.
The first column is an ID
The second is a p-value ( now sorted in increasing order with NA's last)
I want to extract the rows with a p-value of less than 0.05)
hello, I am trying to merge a large number of datasets into one using
"merge.rec".
> merge.rec(db, by.x = "X", by.y = "Y")
>
where "db" is my list of datasets. My problem is trying to create a code to
create the list "db" for all the files without missing data: Currently my
code looks like:
i
I used points graphic (geoR) in this way:
points(zn,pt.divide=c("quintiles"),
main="Location map of Zn",lambda=1,col="gray",yl="coord y",xl="coord x",x.leg=0.
3, y.leg=5.5,dig.leg=2,cex.min=0.5,cex.max=1,)
Can I reduce legend characters
(a sort of cex)?
_
Looking at ?contourplot, it would appear it should be ok to set region=TRUE.
So this is a case of "time to explore".
Have you tried levelplot instead of contourplot? The default for
levelplot is region=TRUE, vs region=FALSE for contourplot. If
levelplot also fails, that tells you something dif
On Apr 16, 2009, at 10:14 AM, Juergen Rose wrote:
Hi,
as a R-newcomer I would like to create some new data frames from a
given
data frame. The first new data frame should content all pairs of the
columns of the original data frame. The second new data frame should
content all tripels of of
Hi
r-help-boun...@r-project.org napsal dne 16.04.2009 15:23:15:
> Hi R,
>
> I have explored R archives a lot but couldn't find an efficient way of
> doing the following:
>
> I want to split a vector into sets of equal sizes. Is there any inbuilt
> function of doing so with the option of speci
On Apr 16, 2009, at 10:41 AM, Andy Barenberg wrote:
hello, I am trying to merge a large number of datasets into one using
"merge.rec".
merge.rec(db, by.x = "X", by.y = "Y")
where "db" is my list of datasets. My problem is trying to create a
code to
create the list "db" for all the file
Petr Pikal
petr.pi...@precheza.cz
724008364, 581252140, 581252257
r-help-boun...@r-project.org napsal dne 16.04.2009 16:45:15:
>
> I have csv files imported in r each with 2 columns and many many rows. I
have
> sorted the data in them but want to extract some values.
>
> The first column is
On Apr 16, 2009, at 10:45 AM, Amit Patel wrote:
I have csv files imported in r each with 2 columns and many many
rows. I have sorted the data in them but want to extract some values.
The first column is an ID
The second is a p-value ( now sorted in increasing order with NA's
last)
I wan
On Apr 16, 2009, at 11:10 AM, David Winsemius wrote:
On Apr 16, 2009, at 10:45 AM, Amit Patel wrote:
I have csv files imported in r each with 2 columns and many many
rows. I have sorted the data in them but want to extract some values.
The first column is an ID
The second is a p-value (
Hi,
I was trying to get error bars in my matplot. I looked at an earlier thread,
and the sample code that I made is:
#--
library(plotrix)
mat1 <- matrix(sample(1:30,10),nrow=5,ncol=2)
ses <- matrix(sample(1:3,10,replace=T),nrow=5,ncol=2)
vect <- seq(20,100,20)
rownames(mat1) <-
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
On 4/16/2009 10:29 AM, Угодай n/a wrote:
I want write article by russian language using Sweave. For cyrillic text
LaTeX use T2A encoding
\usepackage[T2A]{fontenc}
But in Sweave.sty we find:
\RequirePackage[T1]{fontenc}
It is source of critical problem.
For example Rnw file
$ cat estimation
A process like the following is how I would do it:
inputData <- lapply(listOfFiles, function(.file){
input <- read.table(.file, whatever other parameters...)
# now do the modifications that you need
input # return the updated dataframe
})
# combine into one dataframe
inp
http://www.stata.com/help.cgi?egen -- it creates new variables dealing
with some special relatively non-standard tasks that don't boil down
to a one-line arithmetic expressions. For that reason, there will be
no equivalent to -egen- in general, as it has so many functions that
are so different. -ro
[Env: R 2.8.1, Win XP]
For a package I'm working on, I need two small helper functions to
manipulate the
data used in an lm or mlm object, given the *name* of a term, which will
always be
a character string representing
a factor ("A") or an interaction of two or more factors ("A:B", "A:B:C",
Hi R,
Wanted to check if there are any packages available for getting the
(x,y) data points of a graph or a plot, which is in the image format.
Say, the plot could be a published report, and I want to get the points
of the curve plotted. (I am speaking something related the subject
discussed in
Hi all,
I have a question about linear model with interaction:
I created a data frame df like this:
>df
V1 V2 V3 V4 V5
1 6.414094 c t a g
2 6.117286 t a g t
3 5.756922 a g t g
4 6.090402 g t g t
...
which holds the response in the first column and letters (a,c,g,t) in th
Dear R-listers,
I know that there have been many, many posts on the output from
Survreg. To summarise what I have read, Scale is 1/shape of the
Weibull which is also the standard deviation of the normal
distribution which is also the standard deviation of the log survival
time and Intercept is lo
Am Donnerstag, den 16.04.2009, 10:59 -0400 schrieb David Winsemius:
Thanks David,
is there also a shorter way to get the columns names of the new data
frames?
Juergen
> On Apr 16, 2009, at 10:14 AM, Juergen Rose wrote:
>
> > Hi,
> >
> > as a R-newcomer I would like to create some new data fram
Perhaps,
apply(combn(letters[1:4],2), 2, paste,collapse="")
Hope this helps,
baptiste
On 16 Apr 2009, at 17:33, Juergen Rose wrote:
Am Donnerstag, den 16.04.2009, 10:59 -0400 schrieb David Winsemius:
Thanks David,
is there also a shorter way to get the columns names of the new data
frames?
Hi!
> Wanted to check if there are any packages available for getting the
> (x,y) data points of a graph or a plot, which is in the image format.
> Say, the plot could be a published report, and I want to get the points
> of the curve plotted. (I am speaking something related the subject
The grImport package seems to provide such possibility for vector
graphics,
http://www.stat.auckland.ac.nz/~paul/Talks/gddg.pdf
imageJ is another open-source option.
baptiste
On 16 Apr 2009, at 16:44, Shubha Vishwanath Karanth wrote:
Hi R,
Wanted to check if there are any packages avail
Am Donnerstag, den 16.04.2009, 17:41 +0100 schrieb baptiste auguie:
> Perhaps,
>
> apply(combn(letters[1:4],2), 2, paste,collapse="")
>
> Hope this helps,
Thanks Babtiste,
I use now:
Lines <- "abcd
13 015 16
2324250
3334 0 36
0
summary(my_lm) will give you t-values, anova(my_lm) will give you
(equivalent) F-values. summary() might be preferred because it also
provides the estimates & SE.
> a=data.frame(dv=rnorm(10),iv1=rnorm(10),iv2=rnorm(10))
> my_lm=lm(dv~iv1*iv2,a)
> summary(my_lm)
Call:
lm(formula = dv ~ iv1 * iv2,
The following approach works for both of your examples:
xx <- rep(x, length.out=n)
xx[m:n] <- NA
Thus:
> n <- 2
> aa <- rep(a, length.out=n)
> aa[(length(a)+1):n] <- NA
> aa
[1] "2008-01-01" NA
> bb <- rep(b, length.out=n)
> bb[(length(b)+1):n] <- NA
> bb
[1] a
Levels: a
>
R. Raub
Those are not actually dataframes. They are matrices. If you want to
make them into dataframes, use a coercive function. The names can be
generated from the original column names using the same construction
as the column creation:
> apply(combn(colnames(DF),2), 2, paste, collapse="*")
[1] "
The help page for the cnvrt.coords function in the TeachingDemos package shows
some examples of drawing lines/rectangles between multiple graphs created using
par(mfrow=c(...
This function and the examples were written before the grconvertX and
grconvertY functions, so use the example for the g
I apologize in advance that this question is not specific to R, but I thought
some R users may be using this in their work process flow.
I would like to be able to have a tool (prefer "scriptable") that will take two
images and some pre-written text and put it on an simple webpage.
That is,
I am using the package hdrcde to get the highest density region. I have the
data from an unknown distribution. And I used the subroutine hdr from the
package to get the highest density region. But I always got a error message. I
do not know why. Who can help?Thanks!
The codes look like this:
Terse is OK by me as long as I get told what goes in (allowable data
types, argument names and effects) and what comes out. What seemed to
be lacking in that Stata doc for egen was a description of the purpose
or behavior and then could find no description of the values produced.
Perhaps it
Hi Jason,
I apologize in advance that this question is not specific to R, but I thought some R users may be using this in their work process flow.
I would like to be able to have a tool (prefer "scriptable") that will take two images and some pre-written text and put it on an simple webpage.
Great idea - that's a little faster than my previous approach of
setting length() and then re-adding the attributes. Thanks!
Hadley
On Thu, Apr 16, 2009 at 12:16 PM, Raubertas, Richard
wrote:
> The following approach works for both of your examples:
>
> xx <- rep(x, length.out=n)
> xx[m:n] <- N
Hello,
I use the map package in conjunction with the mapdata package,
and I want to draw the map of China without interior regions.
After reading the doc and searched in this mailing list I've tried
map("china",interior=FALSE)
and
map("china",fill=TRUE)
but I allways see the interior regions...
Mike,
I kind of have the same question. What if for a mixed effect model, say
using lme(), how to specify the interaction effect (between a fixed effect
and a random effect)? and where to find the result of the interaction?
Thanks.
Jun
On Thu, Apr 16, 2009 at 12:08 PM, Mike Lawrence wrote:
> su
On Thu, Apr 16, 2009 at 2:01 PM, Peter Dalgaard
wrote:
> Gabor Grothendieck wrote:
>>
>> Using string from another responder's post here are two
>> solutions:
>>
>
> Er, a rather more obvious solution could be to point to the definition of
> substr...
Only if there are a fixed number of digits.
Gabor Grothendieck wrote:
Using string from another responder's post here are two
solutions:
Er, a rather more obvious solution could be to point to the definition
of substr...
1. The first converts to numeric and manipulates that:
[snip]
[37] "0.30" " NA" "0.30" "0.30" "0.30" "0.30" "
On Apr 16, 2009, at 1:39 PM, CheXiaohong wrote:
I am using the package hdrcde to get the highest density region. I
have the data from an unknown distribution. And I used the
subroutine hdr from the package to get the highest density region.
But I always got a error message. I do not know
I'm new to LME myself, so it would be best for others to advise on this.
On Thu, Apr 16, 2009 at 3:00 PM, Jun Shen wrote:
> Mike,
>
> I kind of have the same question. What if for a mixed effect model, say
> using lme(), how to specify the interaction effect (between a fixed effect
> and a random
How would one print the information in a table without having to view it as
a table? I have a dataframe with about 30 columns and 50 rows. About 7 rows
contain human subjects where something is just not right and I need to
manually work out what is going on with them and maybe even call them to
ens
andrew@lshtm.ac.uk wrote:
I want to be able to continuously plot the output from the model in R
each time a new run generates data.
From the C++ program, run the R script that plots the data. Something
like this:
system("Rscript myplotter.R");
That assumes Rscript is in the PA
Folks:
Not to be picky, but depending on exactly what's meant by "attributes," I
think it's impossible:
> x <- 1:6
> attr(x,"length") <- length(x)
> xx <- rep(x,length=10)
> xx[7:10] <- NA
> attr(xx,"length")
NULL
> attr(x,"length")
[1] 6
So "attributes" aren't preserved. The whole point of obj
Hi,
I think this question is best explained using the following
self-contained toy example:
## cut code here and paste to R window
z <- cbind(rep(c("BIC", "hist"), each = 150), rep(rep(c(5, 10, 30), each =
50),2))
z <- as.data.frame(z)
z <- cbind(z, runif(300))
names(z) <- c("Method", "si
Hi Tim,
there are a couple of problems in your example.
(1) The most important is that your 'x' values for the matplot
are 1:5 (that is row numbers of your mat1 matrix)
and are seq(20,100,20) (that is, your vect vector) for your error bars.
Error bars are thus plotted outside the plotting are
Le jeudi 16 avril 2009 à 14:08 -0300, Mike Lawrence a écrit :
> summary(my_lm) will give you t-values, anova(my_lm) will give you
> (equivalent) F-values.
Ahem. "Equivalent", my tired foot...
In simple terms (the "real" real story may be more intricate) :
The "F values" stated by anova are s
The second beginner question. I want to create a new dataframe, where
each element of the original dataframe is translated to 1 if it was "+",
to 0 if it was "-" to -1 otherwise. I could do with:
Lines <- "abcd
+-+ +
+++ -
+1- '+ '
-
Hi,
I'm trying to figure out how to have R estimate partial derivatives for
logit models. As an example, I'm providing a (fake) scored observation in a
MNL with 3 categories of Y and 2 predictors (x01 and x02), and show the
right way to calculate it, but am looking for how to use an R function, s
> Ahem. "Equivalent", my tired foot...
My bad, I wasn't paying attention.
> May I suggest consulting a textbook *before* flunking ANOVA 101 ?
Harsh but warranted given my carelessness.
On Thu, Apr 16, 2009 at 3:47 PM, Emmanuel Charpentier
wrote:
> Le jeudi 16 avril 2009 à 14:08 -0300, Mike La
On 4/16/2009 2:58 PM, Juergen Rose wrote:
> The second beginner question. I want to create a new dataframe, where
> each element of the original dataframe is translated to 1 if it was "+",
> to 0 if it was "-" to -1 otherwise. I could do with:
>
> Lines <- "abcd
> +-+ +
>
Am Donnerstag, den 16.04.2009, 15:14 -0400 schrieb Chuck Cleland:
> On 4/16/2009 2:58 PM, Juergen Rose wrote:
> > The second beginner question. I want to create a new dataframe, where
> > each element of the original dataframe is translated to 1 if it was "+",
> > to 0 if it was "-" to -1 otherwise
Hi,
I think this question is best explained using the following
self-contained toy example:
## cut code here and paste to R window
z <- cbind(rep(c("BIC", "hist"), each = 150), rep(rep(c(5, 10, 30),
each = 50),2))
z <- as.data.frame(z)
z <- cbind(z, runif(300))
names(z) <- c("Method", "sig
I don't know of anything that does the automated process like the example
software you link to, but here is one procedure that you can use to "do it by
hand".
1. Brind the graph into R and plot it (as an image) using rimage or EBImage (or
other) packages.
2. use the locator function to find th
Now that we know what egen is, the answers are one-liners in R:
# Make up some data
vasdat <- matrix ( sample ( 1:100, 3000, replace = TRUE ), ncol = 3 )
# Use apply for each ( MARGIN = 1 means rows, 2 means columns )
anycountresult <- apply ( vasdat, MARGIN = 1, FUN = function ( x ) sum ( x %in%
See, we just jave different expectations of what is to be seen in the
help system, and are used to different formats. Yes, Stata thinks of
data as a rectangular array (although it stores it in memory, unlike
SAS). The inputs to -egen-, as well as the values produced, depend on
the particular functi
Try:
z <- cbind(rep(c("BIC", "hist"), each = 150), rep(rep(c(5, 10, 30),
each = 50),2))
z <- as.data.frame(z)
z <- cbind(z, runif(300))
names(z) <- c("Method", "sigma", "Error")
z$sigma <- factor(z$sigma, c("5", "10", "30"))
library(lattice)
sigma <- as.numeric(levels(z$sigma))
sigmaExprList <-
On Apr 16, 2009, at 3:58 PM, Stas Kolenikov wrote:
See, we just jave different expectations of what is to be seen in the
help system, and are used to different formats. Yes, Stata thinks of
data as a rectangular array (although it stores it in memory, unlike
SAS). The inputs to -egen-, as well
Sorry, can't be done at the moment with the maps package.
The "china" map database is based on data which was provided in
latitude/longitude
rectangles, and has not been properly processed to generate the line segments
and
polygons that the maps package uses to choose individual provinces or t
Duncan Murdoch wrote:
On 4/16/2009 10:29 AM, Угодай n/a wrote:
I want write article by russian language using Sweave. For cyrillic text
LaTeX use T2A encoding
\usepackage[T2A]{fontenc}
But in Sweave.sty we find:
\RequirePackage[T1]{fontenc}
It is source of critical problem.
For example Rnw
Sorry, that should be:
sigma <- as.numeric(levels(z$sigma))
sigmaExprList <- lapply(sigma, function(s) bquote(sigma == .(s)))
sigmaExpr <- as.expression(sigmaExprList)
bwplot(Error~Method | sigma, data = z,
horiz = F, xlab = "Method",
strip = function(which.given, which.panel, var.na
Le jeudi 16 avril 2009 à 13:00 -0500, Jun Shen a écrit :
> Mike,
>
> I kind of have the same question. What if for a mixed effect model, say
> using lme(), how to specify the interaction effect (between a fixed effect
> and a random effect)?
With lme, you have to specify a *list* of random effec
1 - 100 of 126 matches
Mail list logo