Douglas Bates wrote:
> The first VAX at our university, a VAX-11/780, was allegedly a "1
> MIPS" machine but it was pretty difficult to get it to do a million of
> any instruction in one second. I eventually succeeded because it had
> a special instruction for decrement and branch on zero so if yo
Hi,
I was trying to read & filter records from a large file using read.csv.sql.
I was successfully able to do that with a ~1 GB file. However I get the
following error with a >2 GB file which has exactly the same structure as
the first file.
*Error in try({ : *
* RS-DBI driver: (RS_sqlite_import
Hi,
About the forest plot.
Some Phrarma companies demand the report and graphics follow very
restrictive layouts.
SAS allows uses to use one template for graphs and tables.
Margins have to the same for all reports.
Fonts, fontsizes, linewidths, boxing body, cell spacing, cell pad
Hi,
Again I slight mistake on my part. One of the largest Pharmaceutical
Contract companies uses VMS. I erroneously added the VAX in front from my
memory. I don't want to mention the company, but if anyone else there is
familiar with contact pharma companies, one of these uses
VMS(extensively-p
On Mon, Apr 5, 2010 at 3:13 PM, Roger DeAngelis(xlr82sas)
wrote:
>
> Hi,
>
> This is not meant to be critical of R, but is intended as
> a possible source for improvements to R.
> SAS needs the competition.
>
>
> I am reasonably knowledgeable about
>
> R
> SAS-(all products including IML)
>
> SAS
On Mon, Apr 5, 2010 at 8:34 PM, Roger DeAngelis(xlr82sas)
wrote:
>
> Hi,
>
> You are absolutely correct about 32/64 and it appears to be a severe
> penalty. But I think 32 to 32(win/unix) does not incur the penalty. There
> are even more issues between mainframe and unix/windows. The 32 to 64 is
On Mon, Apr 5, 2010 at 12:21 PM, Joris Meys wrote:
> Dear all,
>
> I've read numerous posts about the random and nested factors in lme,
> comparison to proc Mixed in SAS, and so on, but I'm still a bit confused by
> the notations. More specifically, say we have a model with a fixed effect F,
> a r
Hi,
You are absolutely correct about 32/64 and it appears to be a severe
penalty. But I think 32 to 32(win/unix) does not incur the penalty. There
are even more issues between mainframe and unix/windows. The 32 to 64 is a
big hit when querying data dictionaries that have a mixture of 32/64 bit S
Erik Wright wrote:
>
> Hi Charlie,
>
> Thanks, I got it working by looking at your myRunIfConcise function.
>
> SEXP changePercent(SEXP pBar)
> {
> int *rPercentComplete;
> SEXP utilsPackage, percentComplete;
> PROTECT(utilsPackage = eval(lang2(install("getNamespace"),
> Scal
William -
An interesting feature of matrix indexing in R is that
if you provide a two column matrix as a subscript, you are
refering to the elements whose indices are in the rows
of the matrix. This is extremely handy for converting
tables to matrices:
m = cbind(c(1,1,1,2,2,2,3,3,3),c(1,2,
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of Roger DeAngelis(xlr82sas)
> Sent: Monday, April 05, 2010 3:35 PM
> To: r-help@r-project.org
> Subject: Re: [R] SAS and R on multiple operating systems
>
>
> Hi,
>
> I am using
Hi,
One final point about persistent simultaneous environments.
In windows I sit in my development directory(PWD) and simultaneously my
Unix session sits in my production directory(PWD). This simplifies
versioning, promotion to production and batch execution.
My Unix session is persisten
Dear R list,
I have discovered a seemingly peculiar feature when using a matrix to
index itself (yes, this is strange code, which I have now modified to
be more reasonable).
#this makes sense
s <- matrix(1:3,nrow=1)
s[s]#all three elements are shown
#but when I try
s <- matrix(1:2,n
On Mon, 5 Apr 2010, Desmond D Campbell wrote:
Dear Emmanuel,
Thank you.
Yes I broadly agree with what you say.
I think ML is a better strategy than complete case, because I think its
estimates will be more robust than complete case.
For unbiased estimates I think
ML requires the data is MAR,
Hi Charlie,
Thanks, I got it working by looking at your myRunIfConcise function.
SEXP changePercent(SEXP pBar)
{
int *rPercentComplete;
SEXP utilsPackage, percentComplete;
PROTECT(utilsPackage = eval(lang2(install("getNamespace"),
ScalarString(mkChar("utils"))), R_GlobalE
Hi,
One other point.
The connection I have with mutiple servers is persistent the windows SAS
executable is is constant contact with all the SAS server executables.
Also I can submit a job where unix code is interspersed with windows code.
I do execute R and perl from SAS using pipes.
Dear R Users:
I am new to R and I am trying to do a cluster analysis on a
single continuous variable using the Agnes [Agglomerative Nesting (Hierarchical
Clustering) ] in the Package âclusterâ. I was able to apply this
clustering method to
my data:
ward1 <- Agnes(balances, diss= FALSE, m
Hello I am trying to use predict.lm, but I am having trouble getting out of
sample predictions. I am getting the same output if I use the following
three commands:
predict(ModeloLineal,predictors[721:768,])
predict(ModeloLineal,predictors[1:768,])
predict(ModeloLineal)
where ModeloLineal is the o
On Apr 5, 2010, at 5:25 PM, David Winsemius wrote:
>
> On Apr 5, 2010, at 4:58 PM, Ayush Raman wrote:
>
>> Hi all,
>>
>> How can I have a permuted matrix where the second row is the permutation
>> over the first row ; third is the permutation of the second row; forth is
>> the permutation of th
Hi:
I'll use some fake data to show you how to get the plots. To get the data
from Excel into
R, there are several ways to do it: converting the Excel file into csv and
using read.csv() in
R is one method and the XLSReadWrite package is another. Here's a link from
the R Wiki:
http://rwiki.sciviews
Hi,
Here is the SAS command macro that reads what I highlight in my editor and
prints 40 observations from the highlighted dataset after hitting the
function key F4
F4 stores the highlighted text in the clipboard then executes the command
macro, the rsubmit executes the code on the unix zeus s
Hi,
I am using SSH, however I do have to set up a SAS Spawner on the remote
host and use SAS remote library services. I also have to have listeners on
both client and host?
I am not a systems guy, so I do not know exactly how SAS makes the remote
libraries available to windows.
It is a li
On Apr 5, 2010, at 4:58 PM, Ayush Raman wrote:
Hi all,
How can I have a permuted matrix where the second row is the
permutation
over the first row ; third is the permutation of the second row;
forth is
the permutation of the third row and so on ?
Wouldn't any of those permutations just
On Mon, 2010-04-05 at 23:12 +0100, Barry Rowlingson wrote:
> > 3. Highlight and hit F3 and the code is run
> > interactively in unix.
>
> Okay, what's going on here? You have a Windows box (presumably in
> front of you) and a Unix box somewhere on the network. And hitting F1
> runs it on the Win
Dear Emmanuel,
Thank you.
Yes I broadly agree with what you say.
I think ML is a better strategy than complete case, because I think its
estimates will be more robust than complete case.
For unbiased estimates I think
ML requires the data is MAR,
complete case requires the data is MCAR
Anywa
On Mon, Apr 5, 2010 at 9:13 PM, Roger DeAngelis(xlr82sas)
wrote:
>
> Hi,
>
> This is not meant to be critical of R, but is intended as
> a possible source for improvements to R.
> SAS needs the competition.
>
>
> I am reasonably knowledgeable about
>
> R
> SAS-(all products including IML)
> SAS h
Dear JoAnn,
Thank you very much for your reply.
If that is the case I am surprised.
I would have though ML could incorporate study cases with some missingness
in them.
Furthermore I believe ML estimates should generally be more robust than
complete case based estimates.
For unbiased estimates I th
Hi,
You are correct SAS no longer supports OS X under SAS-Proper. I use the
term SAS-Proper for base SAS with SAS-Connect. It does appear that some
improper SAS products are supported under MAC OS?
SAS releases JMP® 8 for Mac, Linux
Users of all major desktop operating systems can now explore
If the columns are in order you can just past in the new names:
Dat <- read.table(textConnection("name x y name1 x1 y1
test1 1 3 test2 4 4
test1 2 2 test2 5 5
test1 3 1 test2 6 6"), header=TRUE)
closeAllConnections()
x.vars <- grep("x", names(Dat))
y.vars <- grep("y", names(Dat))
names.va
> SAS and R run on
>
> Windows(all flavors)
> UNIX(all flavors)
> Apple OSs
According to SAS (http://support.sas.com/kb/33/140.html and
http://support.sas.com/kb/22/960.html), SAS will not run on OS X past 10.4. OS
X 10.5 was released in late 2007, so I don't think it's really fair to say that
Erik Wright wrote:
>
> Hi Charlie,
>
> I like your idea of updating an R progress bar from C, but I don't at all
> understand how to call txtProgressBar from C. I have looked at Writing R
> Extensions and it is equally confusing. Any help would be appreciated.
>
> Thanks!,
> Erik
>
Hi Erik
On Mon, 5 Apr 2010, Roger DeAngelis(xlr82sas) wrote:
SAS and R run on
Windows(all flavors)
UNIX(all flavors)
Apple OSs
I would expect that for more obscure Unices it would be difficult to get SAS,
but basically, yes.
Does R run on natively (no emulation)?
We have quite a few users on these
Hi all,
How can I have a permuted matrix where the second row is the permutation
over the first row ; third is the permutation of the second row; forth is
the permutation of the third row and so on ?
Thanks.
--
Regards,
Ayush Raman
[[alternative HTML version deleted]]
Hi,
I'd like to move from excel to R because our dataset are so large. Here's
what my data looks like:
Transaction Rate Run# Transaction TypeLocationLatency in
Seconds
101 Order
A0
Here is a slight simplification to the first line based on the fact
that na.locf works column by column:
mat <- t(na.locf(t(mydata), fromLast = TRUE, na.rm = FALSE))
On Mon, Apr 5, 2010 at 1:46 PM, Peter Ehlers wrote:
> If I understand correctly what you want (according to your loop),
> you cou
Can anyone tell me why the ncv.test output and the gvlma output would be
contradictory on the question of heteroscedasticity? Below, the ncv.test
output reveals a problem with heteroscedasticity, but the gvlma output says
that the assumptions are acceptable. How is this reconciled?
> ncv.test(de
Hi List,
I want to show the heatmap of a correlation matrix using "heatmap.2", however
always get this warning message (see below) and the column dendrogram is not
showing. It's not really a big deal, but curious how to suppress it and still
let R show what I want to show (i.e. a symmetrical
If I understand correctly what you want (according to your loop),
you could use the na.locf function in pkg:zoo.
library(zoo)
mat <- t(apply(mydata, 1, na.locf, fromLast=TRUE, na.rm=FALSE))
dat <- as.data.frame(mat) ## since apply returns a matrix
-Peter Ehlers
On 2010-04-05 10:52, Anna Ste
Dear Desmond,
a somewhat analogous question has been posed recently (about 2 weeks
ago) on the sig-mixed-model list, and I tried (in two posts) to give
some elements of information (and some bibliographic pointers). To
summarize tersely :
- a model of "information missingness" (i. e. *why* are so
Hi folks,
I have imported data from an Excel spreadsheet. Columns in that spreadsheet
are named "name", "x", and "y", and several sets of those columns appear in
the worksheet. For example:
name x y name x y
test1 1 3 test2 4 4
test1 2 2 test2 5 5
test1 3 1 test2 6 6
When I import these
Hi Charlie,
I like your idea of updating an R progress bar from C, but I don't at all
understand how to call txtProgressBar from C. I have looked at Writing R
Extensions and it is equally confusing. Any help would be appreciated.
Thanks!,
Erik
On Apr 5, 2010, at 10:29 AM, Sharpie wrote:
>
Hi,
This is not meant to be critical of R, but is intended as
a possible source for improvements to R.
SAS needs the competition.
I am reasonably knowledgeable about
R
SAS-(all products including IML)
SAS and R run on
Windows(all flavors)
UNIX(all flavors)
Apple OSs
Does R run on natively (
Thx Erik,
I have no idea what went wrong with the other code snippet, but this one
works.. Appreciate it.
qta<- table(cut(age,breaks = seq(0, 100, by = 10),include.lowest =
TRUE),cut(year,breaks=seq(1950,2010,by=5),include.lowest=TRUE))
M
On 5. apr. 2010, at 21.45, Erik Iverson wrote:
> I d
you can use break :
for(i in 1:5) {
#e.g.
a <- sample(1:10, 1)
# important part :
if(a==3) break
}
a
or while :
a <- 0
while(a != 3){
# an operation that change a :
a <- sample(1:10, 1)
}
Etienne
Le 2010-04-05 10:46, Ravi S. Shankar a écrit :
> Hi R,
>
>
>
> I am running
I don't know what your data are like, since you haven't given a
reproducible example. I was imagining something like:
## generate fake data
age <- sample(20:90, 100, replace = TRUE)
year <- sample(1950:2000, 100, replace = TRUE)
##look at big table
table(age, year)
## categorize data
## see in
On 05/04/2010 3:20 PM, Trafim Vanishek wrote:
Dear all,
My question how is it possible to define a matrix A with 10 rows 1 column,
so that its elements are vectors of undefined length.
I need to have a possibility later to add elements like A[1,1] <-
c(A[1,1],3,4,5)
Thanks a lot for the help!
I may be mistaken, but I don't think that's possible or even should be
possible. A matrix is m x n, where m and n are (kind of fixed) integers. You
cannot have a matrix where m(1) to m(n) (the row lengths) vary. If you want
to do this, you have to use a list instead (I believe).
As a poor workarou
Use lists:
?list (but probably assumes you know what a list is already)
Relevant sections of "An Introduction to R." (Considerable time and effort
have been spent writing this to ease the entry of new users into R. Have you
devoted any time or effort to reading it? )
Bert Gunter
Genentech Noncli
I already did try the regression modeling approach. However the epidemiologists
(referee) turns out to be quite fond of comparing the incidence rates to
different standard populations, hence the need for this labourius approach.
And trying the "cutting" approach I ended up with :
> table (age5)
On 05/04/2010 2:27 PM, Iuri Gavronski wrote:
Hi,
I am trying to compare two 3D plots. For that, I am trying to use the
"pch" parameter in the "points3d" function, but it is not working. Is
it implemented? Any suggestion?
Try reading the help page for points3d. There is no "pch" mentioned.
Dear all,
My question how is it possible to define a matrix A with 10 rows 1 column,
so that its elements are vectors of undefined length.
I need to have a possibility later to add elements like A[1,1] <-
c(A[1,1],3,4,5)
Thanks a lot for the help!
[[alternative HTML version deleted]]
_
You have tempted, and being weak, I yield to temptation:
"Any good ideas?"
Yes. Don't do this.
(what you probably really want to do is fit a model with age as a factor,
which can be done statistically e.g. by logistic regression; or graphically
using conditioning plots, e.g. via trellis graphics
?cut to create categories
?table to make the table
moleps wrote:
Dear R´ers.
I´ve got a dataset with age and year of diagnosis. In order to age-standardize the incidence I need to transform the data into a matrix with age-groups (divided in 5 or 10 years) along one axis and year divided into 5
Dear R´ers.
I´ve got a dataset with age and year of diagnosis. In order to age-standardize
the incidence I need to transform the data into a matrix with age-groups
(divided in 5 or 10 years) along one axis and year divided into 5 years along
the other axis. Each cell should contain the number o
Hi all,
Does any one know how to make Pearson correlation matrix heatmap in R? The
heatmap is a square with highly correlated elements clustered together. And the
heatmap matrix is symmetric with respect to the diagonal line. Many thanks for
your help!
Bill
___
This is the easiest way I have found to do something similar to what you want.
#output of dput() easy way to share data with the list
x <- (structure(list(a = c(1L, 12L, 2L, 0L, 0L, 13L), b = c(4L, 17L,
5L, 2L, 1L, 15L), c = c(7L, 6L, 8L, 4L, 4L, 19L), d = c(9L, 2L,
1L, 7L, 6L, 10L), e = c(2L, 3L,
Hi,
I am trying to compare two 3D plots. For that, I am trying to use the
"pch" parameter in the "points3d" function, but it is not working. Is
it implemented? Any suggestion?
Here goes a reproducible code. I wanted the second plot having
different symbols for the points.
x <- rbind(matrix(rnorm
Dear all,
I've read numerous posts about the random and nested factors in lme,
comparison to proc Mixed in SAS, and so on, but I'm still a bit confused by
the notations. More specifically, say we have a model with a fixed effect F,
a random effect R and another one N which is nested in R.
Say the
Dr. Stevens,
Hi, my name is Trey Scott, and I'm a grad student of Brian McCarthy's. He
referred me to you because of your expertise in handling complex R problems.
We were hoping you could help us solve a nagging problem that is prohibiting
me from producing graphicl output.
Here is a simple m
Thanks Henrique. It works fine.
Cheers!
-Abhi
On Mon, Apr 5, 2010 at 1:41 PM, Henrique Dallazuanna wrote:
> Try this:
>
> DF <- transform(DF, V1 = sprintf('number%d', V1))
>
>
> On Mon, Apr 5, 2010 at 2:38 PM, Abhishek Pratap
> wrote:
> > Hi All
> >
> > I am looking for a way to prefix a consta
Try this:
DF <- transform(DF, V1 = sprintf('number%d', V1))
On Mon, Apr 5, 2010 at 2:38 PM, Abhishek Pratap wrote:
> Hi All
>
> I am looking for a way to prefix a constant value to all the rows in column
> in a data frame.
>
> Eg.
>
> V1
> 2
> 3
> 4
> 5
>
> I want to make it like this
>
> V1
>
Hi All
I am looking for a way to prefix a constant value to all the rows in column
in a data frame.
Eg.
V1
2
3
4
5
I want to make it like this
V1
number2
number3
number4
number5
Thanks!
-Abhi
[[alternative HTML version deleted]]
__
R-help@
Since this may be homework, I'll confine myself to a hint (which
may or may not be the problem; I haven't checked):
The formula you use for z is strongly dependent on the value of 'n'.
-Peter Ehlers
On 2010-04-05 6:06, hix li wrote:
Hi guys,
I have two data sets of prices: endprice0, endpric
On Mon, 5 Apr 2010, pinusan wrote:
Dear R users,
I would like to transform the following "for loop" from R-code to C-code
because it takes really long time to have inc.freqy table.
Unfortunately, I do not have experience to write C code.
Plese, give me some example or advise to transfrom the R
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Duncan Murdoch
> Sent: Monday, April 05, 2010 8:24 AM
> To: anna
> Cc: r-help@r-project.org
> Subject: Re: [R] Deleting many list elements in one time
>
> On 05/04/2010 11:17 AM,
The problem is that your data contains ties, which mess up the nice theory and
result in different people using different approximations.
I don't know where your z-statistic formula comes from, but you can find the
one R uses by looking at the source code in stats:::wilcox.test.default.
To s
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Erik Wright
> Sent: Monday, April 05, 2010 7:59 AM
> To: r-help@r-project.org
> Subject: [R] Rprintf not updating
>
> Hello all,
>
> I am using Rprintf in a C for loop (from .Ca
#I wish to find the "next non-NA" value within each row of a data-frame.
#e.g. I have a data frame mydata. Rows 1, 2 & 3 have soem NA values.
mydata <- data.frame(matrix(seq(20*6), 20, 6))
mydata[1,3:5] <- NA
mydata[2,2:3] <- NA
mydata[2,5] <- NA
mydata[3,6] <- NA
mydata[1:3,]
#this loop acco
Hi R,
I have a column with dates. I need to create a vector say from (current
date-90 days: current date)
For example I need to subtract 90 days from say Sys.Date()-92
If Sys.Date()-92 == "Sunday", Sys.Date()-92+1
if Sys.Date()-92 == "Saturday", Sys.Date()-92+2
i.e if subtractin
Dear R users,
I would like to transform the following "for loop" from R-code to C-code
because it takes really long time to have inc.freqy table.
Unfortunately, I do not have experience to write C code.
Plese, give me some example or advise to transfrom the R to C code.
I have attached the code
Erik Wright wrote:
>
> Hello all,
>
> I am using Rprintf in a C for loop (from .Call) to print a progress
> indicator showing the current percent complete. The loop I am doing is an
> time intensive call to another function. I have noticed that Rprintf does
> not print to the R-window until th
On 05/04/2010 11:17 AM, anna wrote:
Hi guys, here is a simple thing I want to do but it doesn't work:
I have a vector of the element indexes that I want to delete called index
so when I write myList[[index]] <- NULL to delete these elements here is
what I get:
Error in myList[[index]] <- NULL :
Hi guys, here is a simple thing I want to do but it doesn't work:
I have a vector of the element indexes that I want to delete called index
so when I write myList[[index]] <- NULL to delete these elements here is
what I get:
Error in myList[[index]] <- NULL :
more elements supplied than there a
wouldn't using
inside the loop:
if ("a"==3){
# do one thing
} else {
# do another thing
}
Do the trick ?
Contact
Details:---
Contact me: tal.gal...@gmail.com | 972-52-7275845
Read me: www.talgalili.com (Hebrew) | www.biosta
Hello Desmond,
The only way to not drop cases with incomplete data would be some sort
of imputation for the missing covariates.
JoAnn
Desmond Campbell wrote:
Dear all,
I want to do a logistic regression.
So far I've only found out how to do that in R, in a dataset of complete cases.
I'd lik
On Apr 5, 2010, at 9:20 AM, Raadt, Timothy W. wrote:
>
> Hello,
>
> I have a question on the support of the R 2.6 software. We are in the
> process of planning for a hardware refresh and our new machines will be
> running Windows 7 and Internet Explorer 8. My question is if the R 2.6
> so
Raadt, Timothy W. wrote:
>
>
> Hello,
>
> I have a question on the support of the R 2.6 software. We are in the
> process of planning for a hardware refresh and our new machines will be
> running Windows 7 and Internet Explorer 8. My question is if the R 2.6
> software would be supported o
Hello all,
I am using Rprintf in a C for loop (from .Call) to print a progress indicator
showing the current percent complete. The loop I am doing is an time intensive
call to another function. I have noticed that Rprintf does not print to the
R-window until the entire loop has been completed.
On Mon, Apr 5, 2010 at 3:20 PM, Raadt, Timothy W. wrote:
>
> Hello,
>
> I have a question on the support of the R 2.6 software. We are in the
> process of planning for a hardware refresh and our new machines will be
> running Windows 7 and Internet Explorer 8. My question is if the R 2.6
>
Hi R,
I am running a for loop in which I am doing a certain calculation. As an
outcome of calculation I get an out put say "a". Now in my for loop "I"
needs to be initiated to "a".
Based the below example if the output "a"=3 then the second iteration
needs to be skipped. Is there a way to
Gabor Grothendieck wrote:
>
> Passing by value does not necessarily mean physical copying. Check out
> this:
>
>> x <- matrix(1:1000^2, 1000, 1000)
>> gc()
> used (Mb) gc trigger (Mb) max used (Mb)
> Ncells 114520 3.1 35 9.4 35 9.4
> Vcells 577124 4.51901092 14.6
Hello,
I have a question on the support of the R 2.6 software. We are in the process
of planning for a hardware refresh and our new machines will be running
Windows 7 and Internet Explorer 8. My question is if the R 2.6 software would
be supported on a system running Windows 7 and Internet
On Apr 5, 2010, at 8:06 AM, hix li wrote:
Hi guys,
I have two data sets of prices: endprice0, endprice1
I use the Wilcox test:
wilcox.test(endprice0, endprice1, paired = TRUE, alternative =
"two.sided", conf.int = T, conf.level = 0.9)
The result is with V = 1819, p-value = 0.8812.
Then
First thing to notice in your second example is that you have "factor" for
what you think are strings. You should try:
u2<-as.difftime(c(as.character(a),as.character(b)), format ="%Y-%m-%d
%H:%M:%S", units="mins")
and see if this gives you what you are expecting. If you are reading this
in with
Just to follow up on Bert's and Frank's excellent comments. I'm
continued to be amazed by people trying to interpret a single tree.
Besides the variability in the tree structure (try bootstrapping and see
how the trees change), it is difficult to make sense of splits more than
a few levels down (h
On Apr 5, 2010, at 5:32 AM, Nuno Miguel Madeira Veiga wrote:
Hi,
I am working on GLM models. However I am having some problems that I
would like that someone give me some guidance
One of the explanatory variables ERECTANGLE is not present in all
the individual rows.
1 when I delete
I'm new to R and have the following problem with difftime:
if I directly assign date/time strings in difftime I get the expected
result:
> a<-"2010-03-23 10:52:00"
> a
[1] "2010-03-23 10:52:00"
> b<-"2010-03-23 11:53:00"
> u2<-as.difftime(c(a,b), format ="%Y-%m-%d %H:%M:%S", units="mins")
> u2
T
hello,
i need to calculate ci's for each of 4 groups within a dataset, to be able
to infere about differences in the variable "similarity". the problem is
that data within groups is dependent, as assigned by the blocking-factor
"site". my guess was to use a block bootstrap but samples within in t
Passing by value does not necessarily mean physical copying. Check out this:
> x <- matrix(1:1000^2, 1000, 1000)
> gc()
used (Mb) gc trigger (Mb) max used (Mb)
Ncells 114520 3.1 35 9.4 35 9.4
Vcells 577124 4.51901092 14.6 1577448 12.1
> f <- function(x) { y <- max(
Hi guys,
I have two data sets of prices: endprice0, endprice1
I use the Wilcox test:
wilcox.test(endprice0, endprice1, paired = TRUE, alternative = "two.sided",
conf.int = T, conf.level = 0.9)
The result is with V = 1819, p-value = 0.8812.
Then I calculated the z-value of the test: z-v
Thanks to David Winsemius, Peter Ehlers, and Paul Murrell who pointed
out my careless error working with ggplot2's geom_rect(). Not to make
excuses, but when you've done something successfully dozens of times and
suddenly it doesn't work, you're more likely to look for careless errors
on your p
Dear all,
I want to do a logistic regression.
So far I've only found out how to do that in R, in a dataset of complete
cases.
I'd like to do logistic regression via max likelihood, using all the
study cases (complete and incomplete). Can you help?
I'm using glm() with family=binomial(logit).
R contains a whole set of functions to work with text strings, of which
following are definitely worth a look :
?substr
?strsplit
?grep
?gsub
...
Please provide us with an exact problem, and some sample R-code, e.g.
I want to see if a certain word occurs in a vector :
a <- c("a","b","x","c
On 05/04/2010 7:35 AM, Roger Gill wrote:
Dear All,
I would like some advice on creating R packages, passing by reference and oo R.
I have created a package that is neither elegant nor extensible and rather
cumbersome (it works). I would like to re write the code to make the package
distributa
Dear all,
I want to do a logistic regression.
So far I've only found out how to do that in R, in a dataset of complete cases.
I'd like to do logistic regression via max likelihood, using all the study
cases (complete and incomplete). Can you help?
I'm using glm() with family=binomial(logit).
If
Hi,
I am working on GLM models. However I am having some problems that I would like
that someone give me some guidance
One of the explanatory variables ERECTANGLE is not present in all the
individual rows.
1 when I delete the rows for which the variable ERECTANGLE is missing I get
the f
Dear All,
I would like some advice on creating R packages, passing by reference and oo R.
I have created a package that is neither elegant nor extensible and rather
cumbersome (it works). I would like to re write the code to make the package
distributable (should it be of interest) and easy to
Hello together!
I was searching through and different R sites, however I did not find any
useful site for my problem.
Is there any co-occurrence program in R which could work with text strings?
Thanks for any help or hint,
best regards,
Georg
__
R-he
Hi Dilys,
On Fri, Apr 2, 2010 at 3:08 AM, Dilys Vela wrote:
> Hi everyone,
>
> I am doing PCA with labdsv package. I was trying to create a biplot graphs
> in order to observe arrows related to my variables. However when I run the
> script for this graph, the console just keep saying:
>
> *Error
Dear all,
I want to do a logistic regression.
So far I've only found out how, in a dataset of complete cases.
I'd like to do logistic regression via max likelihood, using all the study
cases (complete and incomplete). Can you help?
I'm using glm() with family=binomial(logit).
If any covariate in
100 matches
Mail list logo