2011/4/8 Juan Carlos Borrás :
> #Use the indexes of S in a sapply function.
>
> N <- 10
> S <- sample(c(0,1), size=N, replace=TRUE)
> v1 <- sapply(c(1:N-1), function(i) S[i]&&S[i+1])
You can achieve the same v1 using
v1.2 <- S[2:N-1] & S[2:N]
.. or if you insist on having NA as the first elemen
Dear all,
i have a dataset of about 400 records , with a variable that has two levels 40
bad and 360 good among other variables,how do i come up with10 random samples
that have the composition of as the main sample but maintaining the 40 bad 360
good with replacement, i recently discovered th
Thanks,
The first literal solution does not do what I want: the 'meanplot' was
not a background but a smaller object in the middle of the barcharts.
The second solution is exactly what I want. I just changed the plotting
order and the global ylim in order to have a better plot.
Thanks for you
Hi,
Sorry for my ignorance, but how can I see if a function is 'compiled code' or
plain R? E.g. the daisy
function from the cluster package.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Dear R helpers,
Thanks a lot for your kind help. The xtab function suggested by Mr. Henrique
Dallazuanna produces following output for me
> xtabs(values ~ name + period + cy_date, rat_dat)
, , cy_date = 31-May-11
period
name 29-Mar-11 30-Mar-11 31-Mar-11
A 10.37 10.
Thanks a lot for the answer.
So in fact, without knowing, I was already doing the correct thing ;-)
Regards,
Ivy
-Oorspronkelijk bericht-
Van: Søren Højsgaard [mailto:soren.hojsga...@agrsci.dk]
Verzonden: donderdag 7 april 2011 17:05
Aan: JANSEN, Ivy; r-help@r-project.org
Onderwerp: S
On Apr 8, 2011, at 06:08 , Russ Abbott wrote:
> Haskell is the prototypical lazy evaluation language. One can compute a
> Fibonacci sequence by the Haaskell equivalent of the following R code.
>
>> fibs <- c(0, 1, rep(0, 8))
>> fibs[3:10] <- fibs + fibs[-1]
>
> This works as follows.
>
> fibs
Kenn,
I find your solution more elegant.
2011/4/8 Kenn Konstabel :
> 2011/4/8 Juan Carlos Borrás :
>> #Use the indexes of S in a sapply function.
>>
>> N <- 10
>> S <- sample(c(0,1), size=N, replace=TRUE)
>> v1 <- sapply(c(1:N-1), function(i) S[i]&&S[i+1])
>
> You can achieve the same v1 using
>
>
Hi,
this response uses the previous responses with an example:
#Assume you have 100 observations
n=100
#Simulate a time series of prices
error=rnorm(n,0,3)
raw.price=rpois(n,100)
lag.price=c(rpois(1,100),raw.price[1:99])
price=lag.price+error
#Say you want the moving average based on this
#and
If you want perfect equality, split the data in good and bad and sample from
the two samples individually.
On average, however, random sampling from the entire data will reproduce the
proportion of good and bad in the data.
hth,
Daniel
--
View this message in context:
http://r.789695.n4.nabbl
Hi:
Here's one approach using the reshape package:
library(reshape)
> rat_melt <- cast(rat_dat, period ~ name)
> rat_melt
period A B C D
1 29-Mar-11 10.37 12.61 14.08 11.17
2 30-Mar-11 10.46 12.65 14.12 11.15
3 31-Mar-11 10.42 12.62 14.10 11.11
# To get the dates in reverse
Thank you for your response, but these changes doesn't seem to change
anything, outcomes of effect command is still the same - error.
Tomas
On Fri, Apr 1, 2011 at 5:03 AM, John Fox wrote:
> Dear Tomas,
>
> Write the model as
>
> mreg01 = lm(enep1 ~ enpres * proximity1), data=a90)
>
> That is,
Hi
r-help-boun...@r-project.org napsal dne 08.04.2011 09:31:44:
> Dear all,
> i have a dataset of about 400 records , with a variable that has two
levels
> 40 bad and 360 good among other variables,how do i come up with10
random
> samples that have the composition of as the main sample but
Hi,
I am not perfectly sure what you want to do, but here is what I would do
to maintain good/bad ratio in the sample (as Daniel posted, split the
data and sample from the groups):
df <- data.frame(V1 = 1:400, V2 = c(rep("good",360), rep("bad",40)))
isGood <- which(df$V2=="good")
isBad <- whi
# The code demonstrating the final version I am going to use is as follows
rm(list=ls()) # Beware of this one so it doesn't spoil your workspace
N <- 100
M <- 2
x <- matrix(data=rnorm(N*M, 0, 3)-10, ncol=M, nrow=N)
y <- matrix(c(1,-2,-2,1), ncol=M, nrow=M)
z <- data.frame(x %*% y)
colnames(z) <-
Thank you for your reply.
> forest(log(rr), ci.lb=lci, ci.ub=uci, xlab="Log Relative Risk")
Error in sei^2 : 'sei' is missing
Regards,
Cheba
2011/4/7 Viechtbauer Wolfgang (STAT) <
wolfgang.viechtba...@maastrichtuniversity.nl>
> Actually, Scott suggested ggplot2, which is available on CRAN:
>
>
What version of the metafor package are you using? The option of directly
providing the CI bounds was added in version 1.5-0, so make sure you are using
the newest version.
Best,
--
Wolfgang Viechtbauer
Department of Psychiatry and Neuropsychology
School for Mental Health and Neuroscience
Maa
I have just installed the package as
> install.packages("metafor")
--- Please select a CRAN mirror for use in this session ---
.
downloaded 403 Kb
package 'metafor' successfully unpacked and MD5 sums checked
> library("metafor")
Loading required package: Formula
Loading required package: nlm
Hi,
I have a very simple doubt.
Look:
> teste <- c("A","B","C")
> teste2 <- paste(teste[1],teste[2],teste[3],sep="+")
> teste2
[1] "A+B+C"
>
How to make it automatic, like I try to use paste(teste,sep="+") but the
paste dont get the teste elements separately to join again in a unique
elemen
paste(teste, collapse="+")
ciao!
mario
On 08-Apr-11 12:44, Ronaldo Reis Junior wrote:
teste <- c("A","B","C")
--
Ing. Mario Valle
Data Analysis and Visualization Group| http://www.cscs.ch/~mvalle
Swiss National Supercomputing Centre (CSCS) | Tel: +41 (91) 610
On 04/08/2011 08:44 PM, Ronaldo Reis Junior wrote:
Hi,
I have a very simple doubt.
Look:
> teste <- c("A","B","C")
> teste2 <- paste(teste[1],teste[2],teste[3],sep="+")
> teste2
[1] "A+B+C"
>
How to make it automatic, like I try to use paste(teste,sep="+") but the
paste dont get the test
On 04/08/2011 09:44 AM, Stefan Petersson wrote:
> Hi,
>
> Sorry for my ignorance, but how can I see if a function is 'compiled code' or
> plain R? E.g. the daisy
> function from the cluster package.
>
> __
> R-help@r-project.org mailing list
> https://s
Thanks Daniel,
Its a long way but it will work.
--- On Fri, 4/8/11, Daniel Malter wrote:
From: Daniel Malter
Subject: Re: [R] random sampling with levels and with replacement
To: r-help@r-project.org
Date: Friday, April 8, 2011, 10:08 AM
If you want perfect equality, split the data in
Dear all,
I've tried an lm model. I want to check indepence of the residuals, so I was
trying the following:
library(ape)
plot.correlogram(residuals(model1))
I obtain this error:
Error en x$p.values : $ operator is invalid for atomic vectors
Any idea? Thanks in advance, u...@host.com
--
Hi,
try
paste(teste,collapse="+")
cheers.
Am 08.04.2011 12:44, schrieb Ronaldo Reis Junior:
> Hi,
>
> I have a very simple doubt.
>
> Look:
>
>> teste <- c("A","B","C")
>
>> teste2 <- paste(teste[1],teste[2],teste[3],sep="+")
>> teste2
> [1] "A+B+C"
>>
>
> How to make it automatic, like I tr
I have the following questions about the variance of the random effects in the
survreg() function in the survival package:
1) How can I extract the variance of the random effects after fitting a model?
For example:
set.seed(1007)
x <- runif(100)
m <- rnorm(10, mean = 1, sd =2)
mu <- rep(m, rep(
> Exist any function to associated with paste make this result more
> automatically.
you can try with gsub.
Alfredo
2011/4/8 Ronaldo Reis Junior :
> Hi,
>
> I have a very simple doubt.
>
> Look:
>
>> teste <- c("A","B","C")
>
>> teste2 <- paste(teste[1],teste[2],teste[3],sep="+")
>> teste2
>
I installed R 2.12 and it works now, thank you Wolfgang.
Servus
Cheba
> library(metafor)
Error: This is R 2.11.0, package 'metafor' needs >= 2.12.0
2011/4/8 cheba meier
> I have just installed the package as
>
> > install.packages("metafor")
> --- Please select a CRAN mirror for use in this s
Dear useRs:
I try to plot an rpart object but cannot get a nice tree structure plot. I
am using plot.rpart and text.rpart (please see below) but the branches that
connect the nodes overlap the text in the ellipses and rectangles. Is there
a way to get a clean nice tree plot (as in the Rpart Mayo
On Apr 7, 2011, at 9:46 PM, mavkoup wrote:
I have a for loop with counter i and I'm producing plots where I
ultimately
want to label the axes:
"Entry a_i from P^q", where a, P and q are just letters, and i is
the value
of the for loop counter.
I've tried various combinations of paste, ex
Hi:
I didn't see anything on first blush from the mod1 or summary(mod1) objects,
but it's not too hard to compute:
> names(mod1)
[1] "coefficients" "icoef" "var"
[4] "var2" "loglik""iter"
[7] "linear.predictors" "frail" "fvar"
[10] "df"
Try this:
reshape(rat_dat, direction = 'wide', idvar = 'period', timevar = 'name')
On Fri, Apr 8, 2011 at 4:44 AM, Sandeepa Ramakrishnan <
sandeepa_ramakrish...@yahoo.com> wrote:
>Dear R helpers,
>
> Thanks a lot for your kind help. The xtab function suggested by Mr.
> Henrique Dallazua
Dear ExpeRts,
I am trying to plot a bunch of growth curves and would like to get
some more control over groups and line colors than I seem to have.
Example:
# make some data
dat <- Orange
dat$group <- ifelse(dat$Tree%in%c('1','4','5'), 'A', 'B')
# plot
xyplot(circumference~age, dat, gr
On Apr 8, 2011, at 4:24 AM, Tomii wrote:
Thank you for your response, but these changes doesn't seem to change
anything, outcomes of effect command is still the same - error.
Tomas
On Fri, Apr 1, 2011 at 5:03 AM, John Fox wrote:
Dear Tomas,
Write the model as
mreg01 = lm(enep1 ~ enpres *
Use 'diff' to determine where the changes are:
> S <- sample(0:1,30,TRUE)
> S
[1] 0 0 1 1 0 1 1 1 1 0 0 0 1 0 1 0 1 1 0 1 1 0 1 0 0 0 0 0 1 0
> which(diff(S) == -1)
[1] 4 9 13 15 18 21 23 29
>
then use the indices for the other processing.
On Thu, Apr 7, 2011 at 10:30 PM, Worik R wrote:
> Fr
On Apr 8, 2011, at 8:06 AM, Philipp Pagel wrote:
Dear ExpeRts,
I am trying to plot a bunch of growth curves and would like to get
some more control over groups and line colors than I seem to have.
Example:
# make some data
dat <- Orange
dat$group <- ifelse(dat$Tree%in%c('1','4','5')
Philipp,
I would do the following with ggplot2:
# Set up data
require(ggplot2)
dat <- Orange
dat$group <- ifelse(dat$Tree%in%c('1','4','5'), 'A', 'B')
# Specify the ggplot group aesthetic as Tree
g1 <- ggplot(data = dat, aes(x = age, y = circumference, group=Tree))
# Specify the geom_point and
Hi Tudor,
What type of tree did you create? (classification or regression?)
Could you give an explanation as to what is not looking nice?
Contact
Details:---
Contact me: tal.gal...@gmail.com | 972-52-7275845
Read me: www.talgal
Hadley -
Thanks for your answers. I didnt think of it that way and what you say makes
complete sense. Truthfully, however, I dont care to maintain an equal area for
each increment. When discussing frequency distributions on the rose, having a
diagram with equally spaced intervals would be mor
Dear Tomas,
On Fri, 8 Apr 2011 10:24:45 +0200
Tomii wrote:
> Thank you for your response, but these changes doesn't seem to change
> anything, outcomes of effect command is still the same - error.
It's not really possible to help you with so little information; if you send a
reproducible examp
Hi Chris,
Sorry for the late response. The following selects all students that have
never been suspended:
> with(susim, tapply(sus, id_r, function(x) all(x == 0)))
999881 999886 999890 999892 999896 999897
FALSE TRUE TRUE TRUE FALSE TRUE
> r <- with(susim, tapply(sus, id_r, function(x
Hello R users,
I have a problem to delete rows in a table which are not duplicated in order
of an id number
a short example:
x <- data.frame(cbind(id=c(1,2,2,2,3,3,4,5,6,6), value=1:10))
x_new <- x[which(duplicated(x$id)),]
> x_new
id value
3 2 3
4 2 4
6 3 6
10 610
Hi:
This seems to 'work':
xyplot(circumference~age, dat, groups=Tree, type='l',
col.line = c('red', 'blue', 'blue', 'red', 'red'))
After a little more fiddling around, this also works, and seems a bit less
kludgy:
dat$group2 <- factor(dat$group, labels = c('red', 'blue'))
xyplot(circu
paul4 wrote:
>
> Hello all,
>
> I am new in R language . I want to forecast a univariate time series by
> arima model. Could you please help me finding a good document that
> describes this problem by a simple example and step by step ?
>
?arima
and follow the links in the "See also" section
As I understand it, you are trying to subset the data frame to include only
rows with a non-unique id.
Try this:
x <- data.frame(cbind(id=c(1,2,2,2,3,3,4,5,6,6), value=1:10))
id.table <- table(x$id)
x_new <- subset(x, id %in% id.table[id.table > 1])
Jeremy
___
Running a binary logit model on the data
df <- data.frame(y=sample(letters[1:3], 100, repl=T),
x=rnorm(100))
reveals some residual deviance:
summary(glm(y ~ ., data=df, family=binomial("logit")))
However, running a multinomial model on that data (multinom, nnet)
reveals a residual deviance:
Sorry, I left out the names() function in the last step.
Try this instead:
x <- data.frame(cbind(id=c(1,2,2,2,3,3,4,5,6,6), value=1:10))
id.table <- table(x$id)
x_new <- subset(x, id %in% names(id.table[id.table > 1]))
Jeremy
__
R-help@r-project.org mai
Sir or Madam:
I am new to R and the use of quantile regeression. In addition, I am a
finance person not a true statistcian. Basic regression form is Y =
(Coefficient * Variable) + Error Term
I have results from a quantile regression where I used the Barro and
Roberts method with bootstrapping f
Dear Sir\ Madam
I am trying to maximise a complicated loglikelihood function via EM algorithm
which consists of two step "E-step and M-step"and in this case I need to
maximize the expected of like lihood function " which I get from E- step" and
take those estimates of parameter to update the E
Hi,
I am trying to add text to the bottom of a lattice bwplot with
multiple panels. I would like to add a label below each boxplot, but
the labels do not come from the data. I've tried the following, code:
f1 <- c(rep(c(rep("a", 3), rep("b", 3), rep("c", 3)), 2))
f2 <- c(rep("A", 9), rep("B", 9
Hi!
I'm very new user of R! I need to learn more about! My problems are to
create a time series with exiting data and to make an study:
1 ) I have column with DD-MM-YR, to transform from string text to date,
and use it as to use as time series, I've made the transformation with
command "as.d
I'm working in Rosario, Argentina, trying to plot some contours.
I found the function "draw.contour" created by James Forester at
R-bloggers. If you take a look at that function you'll see it's very
useful!
But I have a problem with it and James is trying to help me (but also
asked me to seek fo
hello
I am at the very beginneing of using the R program
I just don't understand how one can save a programfile
For exemple, if I type in R 23+456 and want to save this file under a ceratin
name to reload it later, i just don't get the way to do it
I can save it with the save function;I achei
I have been trying ot post through R-help and it asks me to subscribe. I have
subscribed may times but did not get a confirmation. can you check to see what
is wrong with my subscription. I can not post on the R-help mailing list at
this
time.
thanks
[[alternative HTML version
Dear users,
I want to set up R to use one R-script directory and one R-library (directory
for packages) for many end-users. I try to do this by using one network-share.
This works fine as long as we don't use the snowfall package with parallel=TRUE
(sfInit(parallel=TRUE, cpus=4, type="SOCK")
Thanks a lot, Jeremy!
It's working perfectly.
With best regards
--
View this message in context:
http://r.789695.n4.nabble.com/Removing-not-duplicated-rows-tp3436600p3436736.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-p
I'm assuming that English is not your first language. In light of that, I
would suggest that you post your question in a language that we both
speak: R. Please include:
1) E-step and M-step functions.
2) The proposed optim/maxLik call.
Mostly because (at least optim) does not have a parameter c
Hi,
On Fri, Apr 8, 2011 at 10:33 AM, gauri j wrote:
> hi,
> I have problem in following code, error is occurred. I have attached my data
> herewith. and my code is as following,
>> library(lars)
> Loaded lars 0.9-8
> Warning message:
> package 'lars' was built under R version 2.12.2
>> x<- read.t
Thanks
That is it really simple.
Ronaldo
Em 08-04-2011 07:48, Mario Valle escreveu:
> paste(teste, collapse="+")
>
> ciao!
> mario
>
> On 08-Apr-11 12:44, Ronaldo Reis Junior wrote:
>> teste <- c("A","B","C")
>
--
1ª lei - Suas férias começam após a defesa e entrega de sua dis
On Fri, 8 Apr 2011, tudor wrote:
Dear useRs:
I try to plot an rpart object but cannot get a nice tree structure plot. I
am using plot.rpart and text.rpart (please see below) but the branches that
connect the nodes overlap the text in the ellipses and rectangles. Is there
a way to get a clean
On 08/04/2011 9:16 AM, nferr...@fceia.unr.edu.ar wrote:
I'm working in Rosario, Argentina, trying to plot some contours.
I found the function "draw.contour" created by James Forester at
R-bloggers. If you take a look at that function you'll see it's very
useful!
But I have a problem with it and
For simple permutation tests I usually just code it up in regular R without
worrying about any packages. R is powerful enough that it is simple to do a
permutation test in only a few lines ( or sometimes just one long line). And
that way you know exactly what it is doing.
--
Gregory (Greg) L
Brittain Briber yahoo.com> writes:
> Do you two or anyone
> else out there happen to know if it is possible to specify the
> width of the intervals in windrose? Is there
> a parameter that I can pass along that would do this?
> Below is a reproducible script in which you can see
> the uneven
I too am encountering this problem. When I have a large script, if I select
all in the editor and then ctrl-r to run, if it encounters a stop() function
it simply prints an error message and continues to execute the remainder of
the script, as opposed to terminating execution at that line. The qu
Look at the updateusr function in the TeachingDemos package.
--
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-bounces@r-
> project.org] On Behalf Of
On 08/04/2011 9:20 AM, DEBERGH Patrick wrote:
hello
I am at the very beginneing of using the R program
I just don't understand how one can save a programfile
For exemple, if I type in R 23+456 and want to save this file under a ceratin
name to reload it later, i just don't get the way to do
On 08/04/2011 11:47 AM, algorimancer wrote:
I too am encountering this problem. When I have a large script, if I select
all in the editor and then ctrl-r to run, if it encounters a stop() function
it simply prints an error message and continues to execute the remainder of
the script, as opposed
Put an open brace in the first line of your file
and a close brace in the last line.
I encourage people with scripts long enough that this
is a problem to divide up their work into functions
that a shorter script calls. (This is akin to UCSD
Pascal on the Osbourne II that refused to deal with
a f
Is anyone aware of a fast way of doing fisher's exact test for a series of 2
x 2 tables in R? The fisher.test is really slow if n1=1000 and n2 = 1000.
--
Thanks,
Jim.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
ht
Some of the functions that were the first in the TeachingDemos package were
originally written to help me visualize something, so it is not just teachers
demoing, but people demoing to themselves. It has become a bit of a misc
package with several utilities that are useful in themselves, but wh
On 08/04/2011 12:40 PM, William Dunlap wrote:
Put an open brace in the first line of your file
and a close brace in the last line.
This shares an advantage with the method I suggested: if there's a
syntax error somewhere in the script, *nothing* will be executed. It
has the disadvantage of
Dear Peter,
Quantile regression is a nice tool but one that requires some statistical
training in order to use it and interpret the results properly. I suggest
backing up a bit.
Frank
Sheldrick,
Peter (Specialty Casualty UW Support) wrote:
>
> Sir or Madam:
>
> I am new to R and the
Would options(error = recover) be of some help?
--
Jonathan P. Daily
Technician - USGS Leetown Science Center
11649 Leetown Road
Kearneysville WV, 25430
(304) 724-4480
"Is the room still a room when its empty? Does the room,
the thing itself have purpose? Or do
Hi:
After a number of false starts, I finally consulted Deepayan's book and the
example on p. 73, suitably adapted, yielded a solution. Add a variable for
the labels and then...
df$lab <- rep(1:6, each = 3)
bwplot(dv ~ f1 | f2, data = df, ylim = c(0.5, 1),
panel = function(x, y, ..., subscripts)
On 08/04/2011 1:29 PM, Jonathan P Daily wrote:
Would options(error = recover) be of some help?
No, that will probably be very confusing. The problem is that the
Windows GUI uses Ctrl-R as a short form of "cut from the editor, paste
to the console", and it will paste the whole text regardless
Hi R community,
I posted a question on using the R maximum likelihood functions a short
while ago and got an email saying that some of the content was
"unprocessed". Hence, I am reposting the question just to be sure (sorry for
the multiple emails if both reached you).
My question is reagarding t
Dear all,
Do you know how to check redundancy of variables and reduce the dimension in R?
It seems PCA can do this. Is there any other better way to implement in R? Many
thanks!
Bill
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/l
Thank you all for the astoundingly quick responses. I think that the
bounding open/closed braces approach sounds like the easiest solution for
the moment -- though I look forward to seeing this all automated in a future
version of R :)
Incidentally, I have indeed encapsulated much of the code as
1. I am not an expert on this.
2. However, my strong prior would be no, since because it is "exact" it has
to calculate all the possible configurations and there are a lot to
calculate with the values of n1 and n2 you gave.
-- Bert
On Fri, Apr 8, 2011 at 9:43 AM, Jim Silverton wrote:
> Is anyon
Hi,
First: Please make sure you CC R-help when responding to messages on
the list so that more people can help you, and more people can benefit
from help that is provided. You can do so easily by replying with
"reply-all."
Now:
Look at the help for the lars function:
R> ?lars
You'll see that `x
I believe the solutions proposed ignore the recursive nature of the
original problem and hence produce wrong solutions.
P <- c(5, 7, 6.01, 6.01, 7)
m <- rep(6, 5)
S0 <- as.numeric(P>(m*1.005))
Then the original loop from Worik gives
S <- S0
for(i in 2:length(S)){
if(S[i]==0 && S[i-1] == 1){
On 2011-04-08 09:26, Ben Bolker wrote:
Brittain Briber yahoo.com> writes:
Do you two or anyone
else out there happen to know if it is possible to specify the
width of the intervals in windrose? Is there
a parameter that I can pass along that would do this?
Below is a reproducible scri
Try
Sys.setenv("LANGUAGE"="Fr")
Uwe Ligges
On 07.04.2011 10:29, Raji wrote:
Hi R-helpers,
Is it possible to localise the error messages/warnings that comes from
R.My application takes in a locale information.I used the following command
to set the locale in R.But in RGui, i still get the e
Greetings All!
After much trial and error, and an exhaustive search of the archives, I'm
writing to solicit help with a problem.
I want to dynamically modify variable names and values as function arguments
inside a loop.
I have a canned function from a package that takes ellipsis (...) arguments
Hi Achim:
I will give it a try and let you know how it goes.
Thanks.
Tudor
--
View this message in context:
http://r.789695.n4.nabble.com/Rpart-decision-tree-tp3435888p3437212.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@
HI,
Thank you, I will try this, and let you know about the results.
Thanks once again.
Gauri
2011/4/8 Steve Lianoglou
> Hi,
>
> First: Please make sure you CC R-help when responding to messages on
> the list so that more people can help you, and more people can benefit
> from help that is prov
I have used R for years but run into a seemingly simple problem involving
'ifelse'. condensed code like this
a=c(2,NA,NA,NA,2,2,NA,2,NA,2)
b=c(NA,1,1,NA,2,2,2,2,2,2)
#I want to combined a and b into c so that c would be a valid number either a
or b is not missing
c=ifelse(a==1|b==1,1,ifelse
Greetings, all ...
I'm trying to find documentation on the use of the "local()" function, and
am drawing a complete blank; would someone please send me a pointer to where
I can learn what it is intended for, and how to correctly use it?
Many thanks!
Doug
[[alternative HTML version delet
Hi Tal:
Thanks for your quick reply.
I tried to create a classification tree.
The problem is that the branches overlap the text in the ellipses (inner
nodes) and rectangles (terminal nodes) - At times this makes it difficult
for the reader to correctly process the results. One way to overcome
On 08/04/2011 2:11 PM, Doug Elias wrote:
Greetings, all ...
I'm trying to find documentation on the use of the "local()" function, and
am drawing a complete blank; would someone please send me a pointer to where
I can learn what it is intended for, and how to correctly use it?
Try ?local. As
Hi,
The problem is that (a == 1 | b == 1) returns some NAs. NAs are not
treated as strictly TRUE or FALSE. From the documentation (see
?ifelse) "Missing values in 'test' give missing values in the result".
You need to use a construct designed to force a TRUE/FALSE value OR
something that handle
On 04/08/2011 11:11 AM, Doug Elias wrote:
Greetings, all ...
I'm trying to find documentation on the use of the "local()" function, and
am drawing a complete blank; would someone please send me a pointer to where
I can learn what it is intended for, and how to correctly use it?
Hi Doug --
?lo
Hello,
I am interested in using the hscat command in the gstat package to create
h-scatterplots of snow depth measurements over a given area (spatial
autocorrelation). Unfortunately I do not understand how to format my data
for import into R for use with the hscat command. I have checked the gstat
Dear Rxperts!
A simple example where "write.csv" does not seem to accept user specified
arguments.. Why?
write.csv(t(1:10),"./te1.csv",quo=F,col.names=F)
Warning message:
In write.csv(t(1:10), "./te1.csv", quo = F, col.names = F) :
attempt to set 'col.names' ignored
However, write.table does
Please read the documentation where this is clearly explained. See ?write.csv
Best,
Ista
On Fri, Apr 8, 2011 at 4:18 PM, Santosh wrote:
> Dear Rxperts!
>
> A simple example where "write.csv" does not seem to accept user specified
> arguments.. Why?
>
> write.csv(t(1:10),"./te1.csv",quo=F,col.nam
Hi. I am trying to run perl script from R. Here's what happens.
> setwd("C:\\Terri\\perl")
> getwd()
[1] "C:/Terri/perl"
> list.files()
[1] "hello.pl" "lines1_10.txt" "try_my.pl" "try_my.txt"
> system("hello.pl")
Warning message:
In system("hello.pl") : Impossible to run C:\Terri\perl\hel
Dear R users,
I am using package mice and I am getting the error "
Error in if (meth[j] != "") { : argument is of length zero." I have tried using
several different versions of R (even the one that will be coming out this
month) to no avail. I am using RStudio as my interface with R. Also note t
I believe write.csv had been revamped and some options are no longer available.
--- On Fri, 4/8/11, Santosh wrote:
> From: Santosh
> Subject: [R] a bug in "write.csv"?
> To: "r-help"
> Received: Friday, April 8, 2011, 4:18 PM
> Dear Rxperts!
>
> A simple example where "write.csv" does not see
Do you mean a test something such as this?
> fisher.test(matrix(c(502,498,490, 510), nrow = 2))
Fisher's Exact Test for Count Data
data: matrix(c(502, 498, 490, 510), nrow = 2)
p-value = 0.6228
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
As pointed out by Ista, please read the help file.
"write.csv and write.csv2 provide convenience wrappers for writing CSV files.
They set sep, dec and qmethod, and col.names to NA if row.names = TRUE and
TRUE otherwise. "
From: John Kane
To: r-help ; San
1 - 100 of 127 matches
Mail list logo
