Dear Rlers,
in the following dataset I would like to insert a new column that refers to
the data of the previous row.
My question is whether the probability of a female (Id) changes if she had
given birth to a pup in the previous year. So my dataframe consists of the
column Id, year (2003-2007 fo
Karl Ove Hufthammer wrote:
In article <4addc1d0.2040...@yahoo.de>, niederlein-rs...@yahoo.de
says...
In every list entry is a data.frame but the columns are only partially
the same. If I have exactly the same columns, I could do the following
command to combine my data:
do.call("rbind", myLis
-- Forwarded message --
From: alexander russell
Date: Tue, Oct 20, 2009 at 4:34 PM
Subject: Re: [R] When modeling with negbin from the aod package...
To: Matthieu Lesnoff
Hello again,
It seems that, though we have a simple estimate of the variance, phi, with
negbin, some models
I think that Richard may be correct. In my experience it's always a
bad idea to use numerical labels for factors, especially when
importing data from Excel. But why not do a str() on your data to
see whether R thinks that time is a factor or not? And if not, why
not convert time to factor before p
Sorry for this third posting - the second method is the same as the
first after all: the coefficients of the first linear model *is* a
linear transformation of the second. Just got confused with the
pasting, tis all.
On Oct 21, 2:51 pm, andrew wrote:
> Oh dear, that doesn't look right at all.
Oh dear, that doesn't look right at all. I shall have a think about
what I did wrong and maybe follow my own advice and consult the doco
myself!
On Oct 21, 2:45 pm, andrew wrote:
> The following is *significantly* easier to do than try and add in
> dummy variables, although the dummy variable a
The following is *significantly* easier to do than try and add in
dummy variables, although the dummy variable approach is going to give
you exactly the same answer as the factor method, but possibly with a
different baseline.
Basically, you might want to search the lm help and possibly consult a
Well I do not know what could be happening.
I tried to append records to a table on DB2 (community Edition 9.5 on
a Linux machine). And it works greatly. The only odd thing (may be
because of my low skills on DB2) is I could not create neither append
to a table on a different schema than dbuser (a
There are two factors of time, but they are evenly replicated across all the
other factors/levels. The experiment is perfectly balanced except for one
lost sample, which is deleted automatically in the aov. I am very certain
the analysis is correct. I think its merely a discrepancy between how a
Here's an example:
x <- runif(100)
sex <- sample(0:1, 100, replace = TRUE)
y <- sex + 2*x + x*sex + rnorm(100)
dfr <- data.frame(y, x, sex)
f1 <- lm(y ~ x, data = dfr, subset = sex == 0)
f2 <- lm(y ~ x, data = dfr, subset = sex == 1)
plot(x, y, pch = sex+1)
abline(f1, col = 2, lty = 2)
abline(f2,
Hello,
I would like to adapt the following code:
> data$datetime3<-format(data$datetime2, tz="EST")
for different time zones. For example, North America's CST and MDT (but
those codes don't work).
I have read ?Sys.timezone but I'm afraid it isn't very helpful. I'm just
looking for a list of the
Now I am worried that you have a wrong analysis.
the aov function is perfectly happy using either factors or
numeric variables. Are there really only two levels of time,
which is what one degree of freedom for time suggests? Or are there
more than two level, but since aov() sees that as a numeri
interesting to know, that might explain it. When I imported I just saved as
CSV and imported with read.csv, I've never done anything to specify the
integers as factors (this is the first time I've used numbers as names).
Here are the precise commands I use(d):
> phen=read.csv(file="phenolics.csv
>I think its a problem with my data, something about how Rexcel
>imported it
We don't have enough information to be sure. My guess is that your
data in Excel is integers which are intended to be levels of a factor.
Excel doesn't distinguish between integers and integers that might be
factor leve
Problem solved:
Lesson learned (I think):
TukeyHSD doesn't like it when you use numbers as names for the factors. ie
factor "time" cannot be "24" and "48" but "twenty-four" and "forty-eight"
work fine.
On Tue, Oct 20, 2009 at 5:36 PM, Clayton Coffman
wrote:
> I've tried that as well, and I get
On Oct 20, 2009, at 4:00 PM, Luciano La Sala wrote:
Dear R-people,
I am analyzing epidemiological data using GLMM using the lmer
package. I usually explore the assumption of linearity of continuous
variables in the logit of the outcome by creating 4 categories of
the variable, performing
Hi,
Maybe use polygon() ?
Good luck,
Megha.
2009/10/20 mnstn :
>
> Hello All,
>
> My question is regarding the attached plot. I would like to have multiple
> transparent green bands running the length (yaxis) of the plot the width of
> which is determined by the green lines at y=0 in the plot. C
I've tried that as well, and I get an error like:
> TukeyHSD(phenolic.aov)
>Error in rep.int(n, length(means)) : unimplemented type 'NULL' in 'rep'
>In addition: Warning messages:
>1: In replications(paste("~", xx), data = mf) : non-factors ignored: time
>2: In replications(paste("~", xx), data = m
ESS does everything you requested.
ess.r-project.org
ess-remote is the specific feature that allows you to run a program
on another computer from within a buffer on the local computer.
Rich
__
R-help@r-project.org mailing list
https://stat.ethz.ch/ma
Hello all,
After reviewing the IDE/Script Editors article at sciviews.org, I
wanted to pose a quick question here to see if anyone can offer an
opinion or commentary about GUI editors that can be installed in a
Windoze environment that allow editing/saving of remote .R files and
running R programs
Hello All,
My question is regarding the attached plot. I would like to have multiple
transparent green bands running the length (yaxis) of the plot the width of
which is determined by the green lines at y=0 in the plot. Can you suggest a
way to do it?
For those who can't or are unwilling to down
On Mon, 19 Oct 2009, Fulton wrote:
>
> I’m running some logistic regressions and I’ve been trying to include weights
> in the equation. However, when I run the model, I get this warning message:
>
> Here’s what it says: Warning message: In eval(expr, envir, enclos) :
> non-integer #successes in
Hi Clayton,
I don't think you need summary().
TukeyHSD(data1.aov)
should work.
-Ista
On Tue, Oct 20, 2009 at 4:17 PM, Clayton Coffman
wrote:
> I can prove I've done this before, but I recently installed Rexcel (and it
> was easiest to reinstall R and some other bits to make it work) and now it
Thanks Andy,
I'm sorry, I didn't clear myself. I was talking about the y-axis, so
your explanation was very helpful.
On Tue, Oct 20, 2009 at 9:39 AM, Liaw, Andy wrote:
> Are you talking about the y-axis or the x-axis? If you're talking about
> the y-axis, that range isn't really very meaningf
I can prove I've done this before, but I recently installed Rexcel (and it
was easiest to reinstall R and some other bits to make it work) and now its
no longer working.
Before I would do an ANOVA and a tukey post-hoc like this:
>data1.aov=aov(result~factor1*factor2, data=data1)
then...
>TukeyH
Hi Axel,
On Tue, Oct 20, 2009 at 7:09 PM, Axel Leroix wrote:
> I estimated a system of equation bu using SUR method. The function summary(xx)
> gives me summary of estimated equation system. However, this function does not
> give my the value of the durbin watson statistic for each one of my
>
Hi,
I understand that the glm.fit calls LINPACK fortran routines instead of
LAPACK because it can handle the 'rank deficiency problem'. If my data
matrix is not rank deficient, would a glm.fit function which runs on
LAPACK be faster? Would this be worthwhile to convert glm.fit to use
LAPACK
To follow up:
I believe I had a version mismatch between one of the packages and R. I made
sure to install the latest of all, at which
point jgr.exe launched just fine. Thanks to Simon and Liviu for their comments.
Just one note: the Windows launcher does not provide any info if it can't reach
Dear R-people,
I am analyzing epidemiological data using GLMM using the lmer package. I
usually explore the assumption of linearity of continuous variables in the
logit of the outcome by creating 4 categories of the variable, performing a
bivariate logistic regression, and then plotting the co
Apologies for a mistake in the announcement. To clarify, for the 2010
Award, nomination deadline is Dec 15, 2009, not Dec 15, 2010 as the
announcement had implied.
From: f...@live.com
To: r-help@r-project.org
Subject: New Award Announcement ASA Stat Comp/Graph Sections: The Statistical
Computing
I am trying to make a scatterplot with containing three columns. I
know how to plot the two columns but now the third column consists of M
or F (male or female) and I don't know how to separate the data so I
can make two separate regression lines on the same graph.
meta #name of file
plot(me
The Statistical Computing and Graphics Award
The ASA Sections of Statistical Computing and Statistical Graphics
have established the Statistical Computing and Graphics Award to
recognize an individual or team for innovation in computing, software,
or graphics that has had a great impact on sta
The query generated by R looks ok to me. I went ahead and pasted it into
my DB2 Client command editor. And I got the message confirmation message.
INSERT INTO STORAGE.TEST_APPEND2 ( "MACRONAME", "MACROUSER", "MACRO_RT" )
VALUES ( 's_ej_mach_config_vz', 'jones2', 5 )
DB2I The SQL command com
Ashta wrote:
Hi Bill and all,
On Tue, Oct 20, 2009 at 12:09 PM, William Dunlap wrote:
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Peter Ehlers
Sent: Tuesday, October 20, 2009 8:48 AM
To: Ashta
Cc: R help
Subject: Re: [R] Counting
How about
unch
On Tue, Oct 20, 2009 at 5:09 AM, r-quantide wrote:
[snip]
> . Is there any methods/functions to obtain the variance components
> for the station factor too? (please, give me some references, at least.)
Pinheiro and Bates 2000 is (practically) a prerequisite for
intelligent use of the nlme
Now it becomes strange.
One thing I note from this generated query:
INSERT INTO STORAGE.TEST_APPEND2 ( "MACRONAME", "MACROUSER",
"MACRO_RT" ) VALUES ( 's_ej_mach_config_vz', 'jones2', 5 )
The names of the variables are in double quotes; That is a problem.
Can you try to run this query on another
Thanks Daniel for your response. Actually I haven't tried that yet,
cause I'm not so fluent with SAS. I got this file from a colleague. I
was thinking that R can read these files, but maybe the problem is
with SAS.
-Johannes
C can you read that dataset just using your Windows SAS v8?
Dan
Dani
On Tue, Oct 20, 2009 at 5:09 AM, r-quantide wrote:
> Dear all,
>
> I'm working with a model to estimate components of variance by using the
> lme() function.
>
>
>
> The model is as the following:
>
> model=lme(fixed=X~1+as.factor(station),data=myData,na.action=na.exclude,rand
> om=~N+1|spliceOrHo
Goran,
You had stated that:
"I have only coded the log likelihood function, and I use 'optim' (in R code)
and 'BFGS', without
gradient, and I ask for a numerically differentiated Hessian, which I use for
calculating standard errors and p-values. Tests show that this works
surprisingly well, bu
Dear Arne Henningsen,
I send you this message because I have question with regard to systemfit
package. I hope you answer to my request.
I estimated a system of equation bu using SUR method. The function summary(xx)
gives me summary of estimated equation system. However, this function
Hi Felipe,
Thanks for your message. Your approach works for an Access database. I am
trying to insert records into a table in a DB2 database on a remote AIX
server.
My userid has the authority to insert records, and I am able to do that
using another application (DB2 Client command line).
Sincerely
The problem probably lies in the source-ing part: look at
getwd()
setwd()
HTH,
Giovanni
> Date: Tue, 20 Oct 2009 13:00:02 +0200
> From: Johan Lassen
> Sender: r-help-boun...@r-project.org
> Precedence: list
> DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma;
> DomainK
Type
?mean
?sd
in the R prompt. For your specific goal, type
mean(r[1:25])
in the prompt. Proceed analogously for sd
Just download one of the very many beginner's manuals or introductory
course materials that you get for free on the internet, and get Ricci's
R-refcard.
Daniel
-
On Tue, Oct 20, 2009 at 9:14 AM, Ted Harding
wrote:
> On 20-Oct-09 13:34:49, Peng Yu wrote:
>> fisher.test() gives a very small p-value, which is underflow on my
>> machine. However, the log of it should not be underflow. I'm wondering
>> if there is a way get log() of a small p-value. An approxim
foo is the data frame you provided
mean(foo[1:25,2])
var(foo[1:25,2])
indexing key:
foo[rows,cols]
hth
Stephen Sefick
On Tue, Oct 20, 2009 at 11:04 AM, Sibusiso Moyo wrote:
> This is day one on R for me, I am trying to figure out how to do simple
> computations. For example I have a data set
Hello Caveman!
I am able to append records using other applications. My credentials have
the authority to append. I am able to query to select records, and even
delete records from the remote DB2 table using RODBC. I'm just stuck when
it comes to APPEND.
The DB2 database is on an AIX server, and
Hi Bill and all,
On Tue, Oct 20, 2009 at 12:09 PM, William Dunlap wrote:
>> From: r-help-boun...@r-project.org
>> [mailto:r-help-boun...@r-project.org] On Behalf Of Peter Ehlers
>> Sent: Tuesday, October 20, 2009 8:48 AM
>> To: Ashta
>> Cc: R help
>> Subject: Re: [R] Counting
>>
>> How about
>>
This is day one on R for me, I am trying to figure out how to do simple
computations. For example I have a data set with 200 observations. I am
trying to compute the mean and variance in r for 1:25 (first 25
observations); 1:50 (first 50 obs) ; 1:100th observation etc. Here is the
dataset:
Id va
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of johannes rara
> Sent: Tuesday, October 20, 2009 12:26 AM
> To: r-help@r-project.org
> Subject: [R] Problems importing Unix SAS .ssd04 file to R (Win)
>
> Hello,
>
> I'm trying to
I've shared a document with you:
plotw5.png
http://docs.google.com/Doc?docid=0AfJ5_yv8GrERZGY3NmhuMnNfMWZtcGhjc2d3&hl=en&invite=CPzN4ZsF
It's not an attachment -- it's stored online at Google Docs. To open this
document, just click the link above.
Nice solution, Bill.
-Peter Ehlers
William Dunlap wrote:
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Peter Ehlers
Sent: Tuesday, October 20, 2009 8:48 AM
To: Ashta
Cc: R help
Subject: Re: [R] Counting
How about
unch <- aggregate(x2==x1, by = list(x
> -Original Message-
> From: Ian Willems [mailto:ian.will...@uz.kuleuven.ac.be]
> Sent: Tuesday, October 20, 2009 6:46 AM
> To: William Dunlap; r-help@r-project.org
> Subject: RE: [R] how to get rid of 2 for-loops and optimize runtime
>
> Hi William,
>
> Your programs works perfect and v
Hi,
I seem to be able to run it without any problems:
> kendall.global(mite.small)
$Concordance_analysis
Group.1
W 0.44160305
F 2.37252221
Prob.F 0.04403791
Chi2 15.89770992
Prob.perm 0.0450
attr(,"class")
[1] "kendall.global"
It seems to be .Rando
On Tue, 20 Oct 2009, Etienne Toffin wrote:
Hi,
I'm analyzing experimental results where two different events ("T1" and "T2")
can occur or not during an experiment. I made my experiments with one factor
("Substrate") with two levels ("Sand" and "Clay").
I would like to know wether or not "Subs
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Peter Ehlers
> Sent: Tuesday, October 20, 2009 8:48 AM
> To: Ashta
> Cc: R help
> Subject: Re: [R] Counting
>
> How about
>
> unch <- aggregate(x2==x1, by = list(x1=x1), FUN = sum)
> chgd <- aggregate(x
Keith,
If you are working within a single time zone, including time zone
information with each record does not seem necessary and you probably are
not recording times to the sub-second. The best solution may thus be to use
a simpler date/time class that does not include time zone information. Y
Hi every body:
I need some help with kendall.global. The example in the manual seems not
working well, and cannot used with my data, always the same error.
data(mite)
> mite.hel <- decostand(mite, "hel")
>
> # Reproduce the results shown in Table 2 of Legendre (2005), a single group
> mite.small
How about
unch <- aggregate(x2==x1, by = list(x1=x1), FUN = sum)
chgd <- aggregate(x2!=x1, by = list(x1=x1), FUN = sum)
-Peter Ehlers
Ashta wrote:
Hi All,
Assume that I have the following data set with two variables and I
want count the number of observation with identical values and num
This is not an explanation but it gives you a solution,
Instead of using lda with a formula do it by giving the variables and the
classification factor as arguments, base on your example and data:
outOfSample <- myDat[11:16,]
train <- myDat[1:10,]
outOfSample <- outOfSample[,3:5]
train2 <
Maybe you're getting strange results because you're not supplying a data object
to lda() when you build your fit.
When I do it the "standard" way, predict.lda() uses the new data and produces a
result of length 6 as expected:
myDat <- read.csv("clipboard", sep="\t")
fit <- lda(c1 ~ v1 + v2 +
Dear R-users,
I tried to read recorded data into R, but confronted the problem about
time zone. The data was recorded in GMT+1 (without summer time
changing), and at first I tried to do in the way:
data <- read.zoo("data.txt", header=FALSE, sep=",", format="%H:%M:%S
%d.%m.%Y", strip.white=TRUE, t
Dear all,
I try to use RGtk2:::gdkColorToString, but it throws an error:
Error in .RGtkCall("S_gdk_color_to_string", object, PACKAGE = "RGtk2") :
gdk_color_to_string exists only in Gdk >= 2.12.0
I know what it means, but don't know to solve this problem because I
don't know where I can download
On Sun, Oct 18, 2009 at 10:29 AM, John Kane wrote:
> Thanks Stefan, the annotate approach works beautifully. I had not got that
> far in Hadley's book apparently :(
>
> I'm not convinced though that the explaination
>
>> you shouldn't use aes in this case since nampost,
>> temprange, ... are not
Dear list,
I've a package, myPkg, which was developed on R v 2.6.1.
I have installed R v2.9.2 from source, with a shared library and tcktk
support. The
package, myPkg, depends on RSPerl, and imports it in its NAMESPACE.
I've downloaded and installed RSPerl using command:
>R CMD INSTALL -c --con
When I import a simple dataset, run LDA, and then try to use the model to
forecast out of sample data, I get a forecast for the training set not the
out of sample set. Others have posted this question, but I do not see the
answers to their posts.
Here is some sample data:
DateNames v1
Hi all!!!
How can I create a legend to my plot that automatically reads the values
from two vectors?
betav<-c(0.78,0.94,0.88,0.41,0.59,4.68)
etav<-c(235.6,59.5,31.2,8.7,3.2,1174)
I want my legend to show 6 lines
beta(in greeks)=0.78,eta(in greeks)=235.6
beta(in greeks)=0.94,eta(in greeks)=6.59
.
Yes that worked perfectly. Many thanks.
From: Peter Ehlers [ehl...@ucalgary.ca]
Sent: 19 October 2009 19:17
To: Weber, Sam
Cc: r-help@R-project.org
Subject: Re: [R] How to get slope estimates from a four parameter logistic with
SSfpl?
Weber, Sam wrote:
>
nobody?
emkayenne wrote:
>
> Hello!
> I am currently using R to deal empirically with copulas. I am using
> financial return data to construct copula models that seem appropriate for
> my data sets. Is there anyone who has done something similar and who is
> interested in talking (or writin
Hi Ashta,
Take a look at ?rle, e.g.
> rle(x1)
Run Length Encoding
lengths: int [1:4] 2 3 1 1
values : num [1:4] 1 0 1 0
> rle(x1)$lengths
[1] 2 3 1 1
> rle(x1)$values
[1] 1 0 1 0
HTH,
Jorge
On Tue, Oct 20, 2009 at 10:10 AM, Ashta <> wrote:
> Hi All,
>
> Assume that I have the following dat
Hi,
I'm analyzing experimental results where two different events ("T1"
and "T2") can occur or not during an experiment. I made my experiments
with one factor ("Substrate") with two levels ("Sand" and "Clay").
I would like to know wether or not "Substrate" affects the occurrence
probability
Nice! I was going to recommend
merge(merge(myList[[1]], myList[[2]], all=TRUE, sort=FALSE),
myList[[3]], all=TRUE, sort=FALSE)
but rbind.fill is better.
-Ista
On Tue, Oct 20, 2009 at 10:05 AM, Karl Ove Hufthammer wrote:
> In article <4addc1d0.2040...@yahoo.de>, niederlein-rs...@yahoo.de
> says
On 20-Oct-09 13:34:49, Peng Yu wrote:
> fisher.test() gives a very small p-value, which is underflow on my
> machine. However, the log of it should not be underflow. I'm wondering
> if there is a way get log() of a small p-value. An approximation is
> acceptable in this case. Thank you!
>
>> fishe
Hi All,
Assume that I have the following data set with two variables and I
want count the number of observation with identical values and number
of time each factor changed from x1 to x2.
x1 x2
1 1
1 0
0 1
0 1
0 0
1 1
0 1
The output should be
x1 changed
In article <4addc1d0.2040...@yahoo.de>, niederlein-rs...@yahoo.de
says...
> In every list entry is a data.frame but the columns are only partially
> the same. If I have exactly the same columns, I could do the following
> command to combine my data:
>
> do.call("rbind", myList)
>
> but of cour
Sorry for being late in responding to this thread, but I was made
aware of it only two weeks ago. In my package 'eha' there is a
function 'aftreg', which performs what is asked for, given that the
time-varying covariates are step functions of time and that the
observation period for each individual
Hello there,
with the following dummy code I'd like to give you an idea how my data
looks like:
df1 <- data.frame(a = rnorm(10), b = rnorm(10))
df2 <- data.frame(a = rnorm(10), c = rnorm(10))
df3 <- data.frame(a = rnorm(10), b = rnorm(10), c = rnorm(10))
myList <- list(df1, df2, df3) # (myLis
Hi William,
Your programs works perfect and very fast for the table I'm using right now
(only one match per row)
If I want to reuse this code other tables, it can match with more than one row.
Is it possible to adapt your code easily, if I have to sum the values of last
month from different rows
Are you talking about the y-axis or the x-axis? If you're talking about
the y-axis, that range isn't really very meaningful. The partial
dependence function basically gives you the "average" trend of that
variable (integrating out all others in the model). It's the shape of
that trend that is "i
fisher.test() gives a very small p-value, which is underflow on my
machine. However, the log of it should not be underflow. I'm wondering
if there is a way get log() of a small p-value. An approximation is
acceptable in this case. Thank you!
> fisher.test(rbind(c(1,10),c(10,1000)))
Hi Joris,
The amount of a month ago is normally one value from another row.
But I used 'sum<-sum + dataset[i,22]' because I would like to reuse the code
also for other tables. In some tables it is possible that the value of last
month is the sum of values from different rows.
Thank u for your t
While developing updates to optimization tools for R (I was behind 3 of
the 5 codes in optim(), but Brian Ripley did the implementation), I've
been seeing this kind of error when the objective function cannot be
computed so returns NA. Examples: attempts to divide by 0 or sqrt(-ve)
or log(0). A
Is the following helpful?
pdd<-deriv(~a+(b-a)/(1+exp((c-t)/d)),"d")
> pdd
expression({
.expr1 <- b - a
.expr2 <- c - t
.expr4 <- exp(.expr2/d)
.expr5 <- 1 + .expr4
.value <- a + .expr1/.expr5
.grad <- array(0, c(length(.value), 1L), list(NULL, c("d")))
.grad[, "d"] <-
Hi Johan,
Although it's not clear to my why you're putting this in a function in
the first place...
The way you've written your function it will only work if your working
directory is the same as the directory where file1.r and file2.r are
stored, and I suspect this is the problem. You can 1) set
Dear all,
I'm working with a model to estimate components of variance by using the
lme() function.
The model is as the following:
model=lme(fixed=X~1+as.factor(station),data=myData,na.action=na.exclude,rand
om=~N+1|spliceOrHoming)
Where X is the response measured variable, station is trea
Dear R community,
You may have the solution to how to construct a function using the function
"source()" to build the function; i.e.
myfunction <- function(...){
source('file1.r')
source('file2.r')
}
After compiling and installing the "myfunction" in R, then calling the
"myfunction" gives an err
Paul Hiemstra wrote:
Remko Duursma wrote:
Hi all,
I am having some problems calling a python script from R that resides
in a folder that is in the path (WindowsXP):
Hi Remko,
Some suggestions:
1. Try to see if the path that R has from a call to system is correct
(i.e. the same as
Remko Duursma wrote:
Hi all,
I am having some problems calling a python script from R that resides
in a folder that is in the path (WindowsXP):
Hi Remko,
Some suggestions:
1. Try to see if the path that R has from a call to system is correct
(i.e. the same as from cmd):
system("path")
Hi all,
I am having some problems calling a python script from R that resides
in a folder that is in the path (WindowsXP):
> system("quickPadTool.py")
Warning message:
In system("quickPadTool.py") : quickPadTool.py not found
# I also tried 'shell' (and shell.exec as well).
> shell("quickPadTool.
Hi Michael,
>> How do you control what is the (intercept) in the model returned by the
>> lme function and is there a way to still be able to refer to all groups
>> and
>> timepoints in there without referring to intercept?
Here is some general help. The intercept is controlled by the contrast
If 'X' is a data frame containing columns 'x', 'y' and 'value, try
m <- with(X, tapply(value, list(y,x), all))
z <- im(m, xcol=sort(unique(X$x)), yrow=sort(unique(X$y)))
w <- as.owin(z)
This will only work if the x, y values form a rectangular grid in some
order.
Adri
On 10/20/2009 06:10 PM, Mao Jianfeng wrote:
Dear R-help listers,
I am now asking for helps on how to draw stacked ellipses to
illustrate the shared and specific of multiple objects using R.
My problem comes from my population genetics study. Now, I genotyped
three species, and I get known about
Hi,all
with my data,there are more than 1000 quantitative results of antibody
concentrations, there may be 2 components(positive and negative), or 3
components (may be strong positive, positive, and negative), or 4-6
components. Could you tell me how to determine the number of components of
t
Hi everybody
I am trying to install RPsSQL and get the following error message:
When I do ./configure form the untarred source directory I get
loading cache ./config.cache
checking for crypt in -lcrypt... no
No crypt function found
When I use the Package installer in R I get
install.packages("/
Hello,
I'm trying to import a SAS file made using SAS on Unix. Currently I'm
using SAS on Windows and I'm trying to import that .ssd04 file to R.
The file name of the file is testfile.ssd04 and it is located in
'M:\sasuser'. I'm using Windows XP and R 2.91. Basically what I'm
doing is
###
I use Lyx (www.lyx.org) with the Sweave noweb report or article class for
the same purpose (if I understand you correctly). LyX is a LateX front-end
where you can embed long and short R code. You can write the complete R
program inside the document. For latex tables the most elegant solution is
som
Dear R-help listers,
I am now asking for helps on how to draw stacked ellipses to
illustrate the shared and specific of multiple objects using R.
My problem comes from my population genetics study. Now, I genotyped
three species, and I get known about the amount of shared and specific
haplotypes
96 matches
Mail list logo
