Hi Ann,
>> I want to delete one category from the dataset ... I have tried the omit
>> command but it just returns TRUE and False values.
You are leaving the list to guess at what you have tried, and which
functions you are using.
There are several different ways of omitting or dropping catego
new R package FactoClass to combine factorial methods and cluster
analysis is uploaded to CRAN. This package is implemented in order to perform a
multivariate exploration of a data table according to Lebart et al. (1995). We
use some ade4 functions (Chessel et al. 2004) to perform the factorial a
Rolf Turner wrote:
> (Note: You *cannot* have spurious objects name "FALSE" (and not equal
> to FALSE) hanging around; R won't let you. That's why you use FALSE and
> not F.)
yes you can, r will let you:
assign("FALSE", TRUE)
ls()
# see "FALSE"
get("FALSE")
# see TRUE
this does not matter much
Wacek Kusnierczyk wrote:
> Rolf Turner wrote:
>
>> (Note: You *cannot* have spurious objects name "FALSE" (and not equal
>> to FALSE) hanging around; R won't let you. That's why you use FALSE and
>> not F.)
>>
>
> yes you can, r will let you:
>
> assign("FALSE", TRUE)
> ls()
> # see "FALS
David,
Are you using the normalized residuals from the $lme part of the gam object
(i.e. something like residuals(foo$lme,type="normalized"))? Without
standardization the raw residuals will look pretty much as bad for the gamm
as they did for the gam (actually they might even lookl a little wo
I get following result, while I write this code :
> lapply(1:5, function(i) c(1,2,3)^i)
[[1]]
[1] 1 2 3
[[2]]
[1] 1 4 9
[[3]]
[1] 1 8 27
[[4]]
[1] 1 16 81
[[5]]
[1] 1 32 243
This is fine. However my goal is : each element of this list should depend
on previous element like :
lis # Lis
Dear all,
does anyone of you know where to get some simple data from dragster racing
(the race where you go for 1/4 mile straight)?
Best would be some data where you have two columns with Col 1 "Time since
start of race in miliseconds" and Col 2 "Distance traveled".
Thank you very much,
Roland
Zitat von David Winsemius <[EMAIL PROTECTED]>:
[...]
> When I used the URL you offered to the Gmane website, both the RSS
> feeds and the web interface links were as current as Gmane ever is.
> Furthermore when I also look at the link entitled "Original list
> archive", it too is also up-to-date.
On Thursday 06 November 2008 05:56:23, Yang Wan wrote:
> I am using the function {aidsEst} in package [micEcon] to do an AIDS
> model now. So far, everything is good. But I want to test the auto
> correlation and heteroskedasticity of the individual equation from AIDS
> demand system. How can I
Sally,
Is Kept.CPUE discrete valued? If it isn't then the poisson family won't like
it, and it would be better to use quasipoisson, if the variance proportional
to mean is the right assumption. (mgcv 1.5 includes a Tweedie family, which
would be an alternative, but it's not ready for release ye
The layout.svd-function in the igraph-package seems to support only
two-dimensional layouts currently. Is anybody aware of a workaround to obtain
the coordinate-values for the third dimension? (And by the way, is it possible
to change the background-color in rglplots?)
I'm using the R-2.6.2 vers
Dear David,
You could put the numbers on the graph with geom_text()
HTH,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and Forest
Cel biometrie, methodologie en
Rolf Turner wrote:
>
> On 17/11/2008, at 1:56 PM, David C. Howell wrote:
>
>> Sampling with and without replacement
>>
>> I seem unable to use "replace = F" when I want to sample without
>> replacement. I would think
>> that it comes down to "F is not a legitimate abbreviation for FALSE."
>> exce
Dear all,
On Windows XP with R 2.8.0, I get the error message below when I try install a
package from the command line.
(Installing as a local zip-file from the menu in the GUI works fine.) Could
anybody point me towards the source (and solution) to the problem?
Regards
Søren
-
I'm trying to plot multiple lines using different colours/symbols to
distinguish them. If I try to plot more than 6 lines, I get an error:
> ggplot( dat, aes(x=time,y=value,group=variable,shape=variable) ) +
geom_line()
Error: scale_shape can deal with a maximum of 6 discrete values, but
yo
Hi, could somebody tell me how to perform non-parametric pairwise
comparisons between related (= paired) series, but only to compare different
series to a control (= not to perform all pairwise comparisons) ?
Many thanks
gilles
__
R-help@r-project.org
Hi Dave,
On Mon, Nov 17, 2008 at 6:35 AM, Dave Murray-Rust
<[EMAIL PROTECTED]> wrote:
> I'm trying to plot multiple lines using different colours/symbols to
> distinguish them. If I try to plot more than 6 lines, I get an error:
>
>> ggplot( dat, aes(x=time,y=value,group=variable,shape=variable) )
On Sun, Nov 16, 2008 at 7:41 PM, Simon Blomberg <[EMAIL PROTECTED]> wrote:
> It is better programming practice to use FALSE for false and TRUE for
> true, and not F and T. This is because it is quite legal to do this:
>
> T <- FALSE
> F <- TRUE
It may be better programming practice, but is it bett
I have used both the trigonometric identity (a^2 + b^2) = c^2, the
haversine and the vincenty equations
(http://www.movable-type.co.uk/scripts/latlong.html
), it really depends on how many calculations you need to make, how
big of an error you can tolerate, and how far apart your measurements
On 17/11/2008 8:03 AM, hadley wickham wrote:
On Sun, Nov 16, 2008 at 7:41 PM, Simon Blomberg <[EMAIL PROTECTED]> wrote:
It is better programming practice to use FALSE for false and TRUE for
true, and not F and T. This is because it is quite legal to do this:
T <- FALSE
F <- TRUE
It may be bet
Hi Everyone,
I am pretty new to R and so far have mostly been using R interactively
through the Windows console.
I'm starting to write some scripts, and have been executing them using the
source() command, i.e. source(myRfile.R).
My questions is how can I pass command line arguments to R. My fi
Brigid Mooney wrote:
Hi Everyone,
I am pretty new to R and so far have mostly been using R interactively
through the Windows console.
I'm starting to write some scripts, and have been executing them using the
source() command, i.e. source(myRfile.R).
My questions is how can I pass command line
Duncan Murdoch wrote:
>
> paramValue <- 15
> source("myRfile.R")
>
> The quotes are necessary, because source(myRfile.R) would go looking
> for a variable named myRfile.R, rather than using "myRfile.R" as the
> filename.
why?
vQ
__
R-help@r-project.org
Is there a better command to use rather than source which would take command
arguments?
I ask because I currently have 6 parameters, will likely have additional
paramaters later, and would like to be able to have default values for each,
if I do not specify new values.
Thanks so much!
On Mon, No
On 17/11/2008 9:14 AM, Brigid Mooney wrote:
Is there a better command to use rather than source which would take
command arguments?
I ask because I currently have 6 parameters, will likely have additional
paramaters later, and would like to be able to have default values for
each, if I do no
Dear all,
I am for the first time trying to work with R, but I have bumped into a
problem.
I have four vectors:
r1 (49 values)
r21 (49 values)
r22 (49 values)
r3 (49 values)
s1 (300 values)
s21 (300 values)
s22 (300 values)
s3 (300 values)
And I would like to calculate the distances from all re
try this (presented only for two pairs):
r1 <- rnorm(49)
r2 <- rnorm(49)
s1 <- rnorm(300)
s2 <- rnorm(300)
dd <- sqrt(0.723523 * outer(r1, s1, "-")^2 + 0.215518 * outer(r2, s2,
"-")^2)
dd
I hope it helps.
Best,
Dimitris
P.Branco wrote:
Dear all,
I am for the first time trying to work with
Rainer,
here is a quick workaround:
layout.svd3 <- function (graph, d = shortest.paths(graph), ...)
{
if (!is.igraph(graph)) {
stop("Not a graph object")
}
l <- svd(d, 3)$u
l[, 1] <- l[, 1]/dist(range(l[, 1]))
l[, 2] <- l[, 2]/dist(range(l[, 2]))
l[, 3] <- l[, 3]/dist(range(l[, 3]
Roland,
I don't know of the specific data you are asking for but the official
website of the 'National Hot Rod Assoc' (?) is www.nhra.com.
DaveT.
*
Silviculture Data Analyst
Ontario Forest Research Institute
Ontario Ministry of Natural Resources
[EMAIL PROTECTE
The TRUE/FALSE vs. T/F issue brings up a related one. Can one assign a
variable a value during an R session that cannot be re-assigned any new
value during the session? That is, the variable is `protected' from
change during the session. `Session' here is not precisely defined, so
that will i
Sorry, it does not work.
If I do a rnorm I lose the original values of my vectors, and the equation
result must be attained by the use of the original values.
Thanks,
P.Branco
Dimitris Rizopoulos-4 wrote:
>
> try this (presented only for two pairs):
>
> r1 <- rnorm(49)
> r2 <- rnorm(49)
> s1
Hi,all
I am doing a nonlinear regression and I am going to use nls.
To Minimize:sum{(Y-BS(s,x,r,t,v)) ^2}
a<-nls(Y~BS(s,x,r,t,v))
BS is the black-scholes model.
BS<-function(s,x,r,t,v){
d1=(ln(s/x)+v^2*t/2)/(v*t^(0.5))
d2=(ln(s/x)-v^2*t/2)/(v*t^(0.5))
c=e^(rt)*(s*N(d1)-x*N(d2))
}
where s,x,
Thanks to you both
--
View this message in context:
http://www.nabble.com/Calculate-Specificity-and-Sensitivity-for-a-given-threshold-value-tp20481633p20541110.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org maili
Dear all,
I want to do some time plots and actually the dates are in the format of
dd/mm/. So first I input my dataframe in R in a csv form. What I do is
DF<-read.csv("C:/Documents and Settings/DF.csv")
DATE<-as.Date(DATE, "%d/%m/%Y") # to tell R that DATE column is indeed dates
with(DF,
"paul murima" <[EMAIL PROTECTED]> writes:
> Hello R-Community,
>
> I am a rookie in R and I am fascinated with the power of bio computing
> by R. I am analysing gene expression data from Real time PCR. I have
> used absolute gene quantitation to measure gene copy number in all my
> transcripts. Al
"Lucke, Joseph F" <[EMAIL PROTECTED]> writes:
>
> The TRUE/FALSE vs. T/F issue brings up a related one. Can one assign a
> variable a value during an R session that cannot be re-assigned any new
> value during the session? That is, the variable is `protected' from
> change during the session.
Dear David,
Thank you very much for your help! I think I have solved the issue.
Lucho
--- On Fri, 11/14/08, David Winsemius <[EMAIL PROTECTED]> wrote:
From: David Winsemius <[EMAIL PROTECTED]>
Subject: Re: [R] Epicalc package
To: "Luciano La Sala" <[EMAIL PROTECTED]>
Cc: "r help"
Date: Fri
Hello,
I have written a mail to gmane
so that they can correct their page (eventually throwing out that link).
Btw, the framed version of the link to the gmane-page (from r-cran project
pages)
appears like being part of the r-cran-pages; that was the reason, why
I thought the gmane-archive-serve
try
plot(DATE, pH, type='l')
On Mon, Nov 17, 2008 at 10:30 AM, Lathouri, Maria
<[EMAIL PROTECTED]> wrote:
> Dear all,
> I want to do some time plots and actually the dates are in the format of
> dd/mm/. So first I input my dataframe in R in a csv form. What I do is
>
> DF<-read.csv("C:/Docum
heatmap.2(m, margins= c(9,9), col = bluered(64), trace=c("none"),
breaks=c(seq(-60,0,60/20)), symkey=TRUE, density.info="histogram",
cexRow=1,)
Error in image.default(1:nc, 1:nr, x, xlim = 0.5 + c(0, nc), ylim = 0.5 + :
must have one more break than colour
Hie all.
I am working on scaling my d
On Mon, Nov 17, 2008 at 07:20:22AM -0800, P.Branco wrote:
>
> Sorry, it does not work.
>
> If I do a rnorm I lose the original values of my vectors, and the equation
> result must be attained by the use of the original values.
He didn't want you to actually do the rnorm part - he just used
that
On Mon, 17 Nov 2008 15:30:14 - "Lathouri, Maria"
<[EMAIL PROTECTED]> wrote:
LM> I want to do some time plots and actually the dates are in the
LM> format of dd/mm/. So first I input my dataframe in R in a csv
LM> form. What I do is DF<-read.csv("C:/Documents and Settings/DF.csv")
LM> DATE<
Have a look at ?lockBinding
Best,
luke
On Mon, 17 Nov 2008, Lucke, Joseph F wrote:
The TRUE/FALSE vs. T/F issue brings up a related one. Can one assign a
variable a value during an R session that cannot be re-assigned any new
value during the session? That is, the variable is `protected' f
Also if this is a time series you may wish to represent it as such
to facilitate other computations as well. This assumes that the
DATE column is first and the remaining columns are numeric:
library(zoo)
mypath <- "/whatever/myfile.csv"
z <- read.zoo(mypath, sep = ",", header = TRUE, format = "%d
Lathouri, Maria imperial.ac.uk> writes:
[...]
> Until here it works fine, but I have the graph plotting only the points. What
I want is to have a line (join
> these points) so to have a time plot. I have tried different commands such as
lines(DATE,pH) or with(DF,
> lines(DATE,pH) but nothing work
On Mon, 2008-11-17 at 15:12 +0100, Wacek Kusnierczyk wrote:
> Duncan Murdoch wrote:
> >
> > paramValue <- 15
> > source("myRfile.R")
> >
> > The quotes are necessary, because source(myRfile.R) would go looking
> > for a variable named myRfile.R, rather than using "myRfile.R" as the
> > filename.
>
On Mon, 2008-11-17 at 07:20 -0800, P.Branco wrote:
> Sorry, it does not work.
>
> If I do a rnorm I lose the original values of my vectors, and the equation
> result must be attained by the use of the original values.
Dimitris was generating some dummy data to test that the function
worked. He ha
On 17 Nov 2008, at 13:00, hadley wickham wrote:
Hi Dave,
On Mon, Nov 17, 2008 at 6:35 AM, Dave Murray-Rust
<[EMAIL PROTECTED]> wrote:
I'm trying to plot multiple lines using different colours/symbols to
distinguish them. If I try to plot more than 6 lines, I get an error:
ggplot( dat, aes(x
Hi,
Does anyone know of an R ORM (Object Relational Mapper)? I'm thinking of
something similar to sqlalchemy (http://www.sqlalchemy.org/).
Alternatively or additionally, can people offer suggestions about managing
relational databases from R. I'm currently using postgresql, but would
like a
On Nov 17, 2008, at 11:01 AM, paul murima wrote:
heatmap.2(m, margins= c(9,9), col = bluered(64), trace=c("none"),
breaks=c(seq(-60,0,60/20)), symkey=TRUE, density.info="histogram",
cexRow=1,)
Error in image.default(1:nc, 1:nr, x, xlim = 0.5 + c(0, nc), ylim =
0.5 + :
must have one more br
On Mon, Nov 17, 2008 at 5:00 AM, megh <[EMAIL PROTECTED]> wrote:
>> lapply(1:5, function(i) c(1,2,3)^i)
> [[1]]
> [1] 1 2 3
...
> This is fine. However my goal is : each element of this list should depend
> on previous element like :
>
> lis # List name
> then,
>
> lis[[i]] = lis[[i-1]] + c(1,2,3)^
On Nov 17, 2008, at 12:30 PM, David Winsemius wrote:
On Nov 17, 2008, at 11:01 AM, paul murima wrote:
heatmap.2(m, margins= c(9,9), col = bluered(64), trace=c("none"),
breaks=c(seq(-60,0,60/20)), symkey=TRUE, density.info="histogram",
cexRow=1,)
Error in image.default(1:nc, 1:nr, x, xlim = 0
To save my fingers and still being on the safe side, I always do:
> !0
[1] TRUE
> !1
[1] FALSE
;) ...still hackable though.
/Henrik
On Mon, Nov 17, 2008 at 5:25 AM, Duncan Murdoch <[EMAIL PROTECTED]> wrote:
> On 17/11/2008 8:03 AM, hadley wickham wrote:
>>
>> On Sun, Nov 16, 2008 at 7:41 PM, S
Hi All,
I am making a R package which is using a boost c++ library. Can some
suggest how can I compile the dll in windows using boost. It makes a
package but during installation in fails to link boost and give
errors
Any help will be useful
Thanks
Vidhu
_
Hi. I have built a levelplot with 3 variables, X, Y, and Z where X and Y are
the two axes and Z represents the intensity (i.e. Z~X*Y).
Now I want to adjust the size of the grid (like a mosaic plot) where the size
of each grid is weighted by a variable, W.
Just wonder how can I do that with lev
Specifically it is for a MNL I just want to drop temporarily to run a
test. I have tried upData (dataset, g=as.numeric (g) but then I am not sure
how I could apply na.omit to omit the variable. How are you definng
"response" and "predictor"?
On Mon, Nov 17, 2008 at 3:42 AM, Mark Difford <[EMAIL
Susan Ofner
Biostatistician II
Division of Biostatistics
Indiana University School of Medicine
410 West 10th St., Suite 3045
Indianapolis, IN 46202
(317) 278 - phone
(317) 274 - 2678 FAX
http://www.biostat.iupui.edu/
[[alternative HTML version deleted]]
___
Using the data set fgl in MASS the following code
layout(matrix(1:9,3,3))
for(i in 1:9){
boxplot(fgl[,i] ~ type, data = fgl,main=dimnames(fgl)[[2]][i])}
produces a 3 by 3 array of plots, each one of which consists of six
boxplots.
Is it possible to do this in lattice?
Steve
"R version 2.7.
Hi Søren,
http://jekyll.math.byuh.edu/other/howto/R/Rcmdr.shtml
Søren Højsgaard wrote:
>
> Dear all,
>
> On Windows XP with R 2.8.0, I get the error message below when I try
> install a package from the command line.
> (Installing as a local zip-file from the menu in the GUI works fine.)
I need help generating all possible combination of a vector of randomly
assigned values of 1 and -1. For example, a vector of 20 randomly placed 1s
and -1s, and all possibilities (which would amount to 2^20 total vectors). I
am able to generate one such vector (via sample(c(-1,1),20,replace=TR
Hello;
I have a function with 3 variables. I would like to take integration with
respect to just one variable. How can I write this function? You can find my
function above. (
I am trying to find \int_ZL^ZU
rho.bisquare(rho*z+u*sqrt(1-rho^2)*dnorm(u)du.Integrand has three variables, I
know the
On 11/17/2008 1:50 PM, steve wrote:
> Using the data set fgl in MASS the following code
>
> layout(matrix(1:9,3,3))
> for(i in 1:9){
> boxplot(fgl[,i] ~ type, data = fgl,main=dimnames(fgl)[[2]][i])}
>
> produces a 3 by 3 array of plots, each one of which consists of six
> boxplots.
>
> Is it pos
How about
!!0 # FALSE
!!1 # TRUE
Although its one more char, its very easy to press ! twice and when you
look at it its more obvious since 0 is associated with FALSE and
1 with TRUE.
On Mon, Nov 17, 2008 at 12:42 PM, Henrik Bengtsson <[EMAIL PROTECTED]> wrote:
> To save my fingers and still bein
Hi, does anybody know if it is possible to use an image (rather than
text) to label the leaves of a dendrogram? I realize that this will
not always lead to a useful on-screen display, but ideally I'd like
to directly go from the dendrogram to a PDF and on-screen viewing is
not to too import
## Question1: How to define IV with interaction alone, without main effects?
## Question2: Should Type III ANOVA in package car be independent of
the factor level order?
## data from http://www.otago.ac.nz/sas/stat/chap30/sect52.htm
drug <- c(t(t(rep(1,3)))%*%t(1:4));
disease <- c(t(t(1:3)) %*% t
This gives a data frame rather than a matrix, but you can convert it:
x <- do.call( expand.grid, rep( list( c(-1,1) ), 20 ) )
hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
801.408.8111
> -Original Message-
> From: [E
On Mon, Nov 17, 2008 at 11:15 AM, Chuck Cleland <[EMAIL PROTECTED]> wrote:
> On 11/17/2008 1:50 PM, steve wrote:
>> Using the data set fgl in MASS the following code
>>
>> layout(matrix(1:9,3,3))
>> for(i in 1:9){
>> boxplot(fgl[,i] ~ type, data = fgl,main=dimnames(fgl)[[2]][i])}
>>
>> produces a 3
Hi Susan,
This is one of the things you are required to do yourself. Head over to:
http://stat.ethz.ch/mailman/listinfo/r-help
and look towards the bottom of the page for unsubscribe options.
If that doesn't work for you then contact the list admins, email
addresses for them are located towards
The help page for the subplot function in the TeachingDemos package shows an
example of adding an image to an existing plot.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
801.408.8111
> -Original Message-
> From: [EMAIL PROTECTED] [
Hi Ann,
>> Specifically it is for a MNL ...
>> How are you definng "response" and "predictor"?
You really are still leaving me, and I imagine other people on the list who
read your email, largely in the dark. (But since you don't seem to be
prepared to address anyone, or even the list, I don't s
Thank you very much for the suggestions, however I do not think they will
help. Please let me know if I misunderstand.
1) Just using "colSums" without the "(!is.na)" portion provides:
Error in FUN(X[[1L]], ...) :
'x' must be an array of at least two dimensions
2) In that context, sum appears
Win xp sp2, R v2.7.1
Hi. If I have two numeric columns in a data frame, I can use the paste
command to combine them into a new column separated by a comma.
c3=paste(c1,c2,sep=',')
gives: 1 1 -> "1,1"
Is there any way I can use a new line (\n) as a separator?
i.e.
1 1 -> 1
1
I tried paste
On Sun, Nov 16, 2008 at 7:27 PM, Stavros Macrakis <[EMAIL PROTECTED]> wrote:
> Thank you for the pointer to your paper "Getting Started with Lattice
> Graphics" -- it looks very useful.
>
> I think we're talking past each other on the question of the semantics
> of formula operators, and it's proba
My question is probably pretty basic, but since I'm really new to R, here it
goes
I have two separate data frames that include class names and various other
information on classes. I'm trying to create a match based on class names
and if a match exists to create a third data frame with the c
[Using R 2.8.0 / Win XP / ]
I just added a CITATION file to the heplots package--- appended below.
From the document ion for ?CITATION, there can be *one or more* calls to
citEntry() within the CITATION file, and each should produce an object
of class "citation".
With just a single citEntry(), ci
I'm not sure I understand what you're looking for in the result. What
exactly do you mean by a "match"? What do you want in the third table
besides the class names?
Do you just want a list (not a data frame) of those class names which
table A and table B have in common? Then how about
intersect
Dear Sephanie,
Try the following:
# Data set including some NA values
set.seed(123)
X=cbind(c(3, 2, 2, 1, 2, 5, 5, 4, 1, 1),matrix(rnorm(100),ncol=10))
colnames(X)=c('day',paste('X',1:10,sep=""))
X[2,2]<-X[5,10]<-X[3,8]<-NA
X
# mean excluding NA values
# in parenthesis are the number of observati
Alina
I suspect you want match (or %in%):
?match
HTH
Peter Alspach
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Alina Sheyman
> Sent: Tuesday, 18 November 2008 9:04 a.m.
> To: r-help@r-project.org
> Subject: [R] looking for matches
>
I have a data.frame containing survey data that is already organized
like a table:
> foo
Excellent Very Good Good Fair Poor Very Poor
Question 1 8 7300 0
Question 2 5 5710 0
Question 3 71010
Dear slurpy,
Perhaps:
cat('', 1,'\n', 1,'\n')
See ?cat for more information.
HTH,
Jorge
On Mon, Nov 17, 2008 at 3:03 PM, slurpy <[EMAIL PROTECTED]> wrote:
>
> Win xp sp2, R v2.7.1
> Hi. If I have two numeric columns in a data frame, I can use the paste
> command to combine them into a new c
Dear Dr. Carbon,
Perhaps (not lattice):
barplot(t(foo),beside=TRUE,col=rainbow(6),ylim=c(0,max(foo)),names.arg=paste('Question',1:3),
ylab='Your title here',xlab='Your title here')
legend('topleft',colnames(foo),pch=15,col=rainbow(6),ncol=2,cex=1.1)
HTH,
Jorge
On Mon, Nov 17, 2008 at 3:40 PM
slurpy gmail.com> writes:
>
>
> Win xp sp2, R v2.7.1
> Hi. If I have two numeric columns in a data frame, I can use the paste
> command to combine them into a new column separated by a comma.
> c3=paste(c1,c2,sep=',')
> gives: 1 1 -> "1,1"
> Is there any way I can use a new line (\n) as a sepa
Hello,
Does anyone know of R code that exists to correct for sample selection
with panel data as in:
J.M. Wooldridge (1995), “Selection Corrections for Panel Data Models
Under Conditional Mean Independence Assumptions,” Journal of
Econometrics 68, 115-132.
Many thanks for your consideratio
the following is a trivialized version of some functional code i tried
to use in r:
(funcs = lapply(1:5, function(i) function() i))
# a list of no-parameter functions, each with its own closure environment,
# each supposed to return the corresponding index when applied to no
arguments
sapply(func
Thank you. Here's my version, using melt instead of do.call(make.groups...
library(reshape)
fgl2 = melt(fgl[,-10])
fgl2$type = fgl$type
bwplot(value ~ type | variable, data = fgl2)
Steve
Deepayan Sarkar wrote:
On Mon, Nov 17, 2008 at 11:15 AM, Chuck Cleland <[EMAIL PROTECTED]> wrote:
On 11
On Mon, 17 Nov 2008, Michael Friendly wrote:
[Using R 2.8.0 / Win XP / ]
I just added a CITATION file to the heplots package--- appended below.
From the document ion for ?CITATION, there can be *one or more* calls to
citEntry() within the CITATION file, and each should produce an object
of clas
Thanks a lot. That works. I had no clue that the '\n' was still active and
I could use cat to see the result.
On Mon, Nov 17, 2008 at 3:58 PM, Oliver Bandel <[EMAIL PROTECTED]>wrote:
> slurpy gmail.com> writes:
>
> >
> >
> > Win xp sp2, R v2.7.1
> > Hi. If I have two numeric columns in a data
so just to show how scheme would do that:
(map
(lambda (function) (function))
(map
(lambda (index) (lambda () index))
'(1 2 3 4 5)))
;; *this* is the scheme semantics
vQ
xxx wrote:
> I agree, something is wrong here, especially since the R community
> likes to say that R semant
In lattice
#toy data
library(ggplot2)
library(lattice)
x <- rnorm(100)
y <- rnorm(100)
k <- sample(c("Weak","Strong"),100,replace=T)
j <- sample(c("Tall","Short"),100,replace=T)
w <- data.frame(x,y,j,k)
xyplot(y~x|j+k,scales=list(y=list(relation="free")))
will give you a scale in each subplot w
Hello list member:
I've recently had a problem in that I'm unable to quit an R Session. I
noticed this after the update to 2.8.0, but I believe I also noticed it
on another machine, in the previous version. It occurs on both linux
and Mac platforms. It only occurs when I start R in some parti
> "N" == Nick <[EMAIL PROTECTED]>
> on Sat, 15 Nov 2008 13:39:44 -0800 (PST) writes:
N> Code: # svm is an S3 class, so: setOldClass("svm")
N> # Create the super class model: setClass("model")
N> # For svm from e1071 this works (well gives no error):
N> setIs("svm", "
Wacek,
I think when people say that R semantics are derived from Scheme, all
they mean is that R supports lexical closures. But R has other
features which are very un-Scheme-like, and when they interact with
lexical closures, you get behavior you don't find in other functional
languages.
R passe
On Mon, 17 Nov 2008, Wacek Kusnierczyk wrote:
the following is a trivialized version of some functional code i tried
to use in r:
(funcs = lapply(1:5, function(i) function() i))
# a list of no-parameter functions, each with its own closure environment,
# each supposed to return the correspondin
Stavros Macrakis wrote:
> Wacek,
>
> I think when people say that R semantics are derived from Scheme, all
> they mean is that R supports lexical closures. But R has other
> features which are very un-Scheme-like, and when they interact with
> lexical closures, you get behavior you don't find in o
On Mon, Nov 17, 2008 at 3:42 PM, steve <[EMAIL PROTECTED]> wrote:
> Thank you. Here's my version, using melt instead of do.call(make.groups...
>
> library(reshape)
> fgl2 = melt(fgl[,-10])
> fgl2$type = fgl$type
> bwplot(value ~ type | variable, data = fgl2)
Or even more succintly:
fgl2 <- melt(f
To solve the immediate problem of quitting R, you could try to define
a dummy of the missing function, e.g.
validServerIsRunning <- function(...) FALSE;
and then quit(). This does of course ignore any potential problems of
not "finalizing" the 'fame' package correctly. Alternatively, since
the
On Mon, Nov 17, 2008 at 3:49 PM, Etches Jacob <[EMAIL PROTECTED]> wrote:
> In lattice
>
> #toy data
> library(ggplot2)
> library(lattice)
> x <- rnorm(100)
> y <- rnorm(100)
> k <- sample(c("Weak","Strong"),100,replace=T)
> j <- sample(c("Tall","Short"),100,replace=T)
> w <- data.frame(x,y,j,k)
>
>
## I got it. IV(s) of interaction should be orthogonal to main effect IV(s).
## Type III ANOVA / Interaction alone
x_interaction<-cbind(
(drug==2)&(disease==2)
,(drug==3)&(disease==2)
,(drug==4)&(disease==2)
,(drug==2)&(disease==3)
,(drug==3)&(disease==3)
Fair enough. But I find my interactive data analysis jobs quickly get
big enough (data manipulation, a series of model fits, some customised
output) for the analysis script to turn into something that looks like a
program. Of course, YMMV. I also get annoyed at code that uses = for
assignment outsi
Hi All, thank you so very much for your help, i have now got it
working! I thought that i had replied already but i don't think it got
through so this is a repost of it for anyone who does a search on this
topic...
After adding the directory to the path variable, i should have
restarted my laptop.
1 - 100 of 123 matches
Mail list logo