Hi
Most methods can somehow deal with NA values. You can find many packages
regarding discriminant analysis by simply searching CRAN.
discriminant analysis site:r-project.org
Or you can look to CRAN task view Multivariate statistics
maybe package pls or DiscriMiner does what you want.
You ca
Thanks for this guidance.
In light of your advice I have reduced the titanic to a dingy by reducing the
size of my sample to just 100 instances, making this an optimisation of 100
parameters. I am, however, seeing similar output when I output the function
value at each evaluation. This is the
The CRAN master machines are being moved to the new WU campus in Vienna
today. The services will be partially or completely unavailable,
including www.r-project.org, ftp uploads for packages and the web
submission form. We expect services to be available again this evening
European time.
The
Hi,
What should I do if I want to use "non-default" R-options for the
functions inside a package, but not affect the options of the rest of
the session? Specifically, I'm thinking of using "sum contrasts" instead
of the default, in functions like "lm", etc. when they are called by
other functions
a 'split' function should not take hours to complete. I will not even wait a
couple of minutes for one until I look at other approaches. Can you at least
show what the command looks like and what the structure of the data is. Are
you trying to split a dataframe on multiple criteria? I tend t
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Alok Jadhav
> updated code pasted here. So I guess CC is not an option at
> all? In that case is there any other package that would be useful?
I'm guessing here (based on experience of PHP's sendmail inter
Hi all,
several packages print messages during loading. How do I avoid to see
them when the packages are in the defaultPackages?
Here is an example.
With this in ~/.Rprofile
,[ ~/.Rprofile ]
| old <- getOption("defaultPackages")
| options(defaultPackages = c(old, "filehash"))
| rm(old)
`---
Look at the package chemometrics, it can certainly handle your number of
variables (p > n is what that's called and it requires special considerations).
I don't recall about missing values. The authors of that package also have a
very helpful text. Good Luck. Bryan
On Jul 14, 2013, at 10:50
On 13-07-15 5:27 AM, Helios de Rosario wrote:
Hi,
What should I do if I want to use "non-default" R-options for the
functions inside a package, but not affect the options of the rest of
the session? Specifically, I'm thinking of using "sum contrasts" instead
of the default, in functions like "lm
Dear list,
I have got a weird encoding problem with the xmlEventParse() function
from the 'XML' package.
I tried finding an answer on the web for several hours and a Stack
Exchange question came back without success :(
So here's the problem. I created a small XML test file, which looks like
Hello,
You should Cc the list.
What I meant is that, at an R prompt you should type
install.packages(pkgname)
substituting 'pkgname' by the name of the package you're using.
Rui Barradas
Em 15-07-2013 12:01, waqas shafqat escreveu:
Sorry sir
i could not find the install.packages(pkgname) in
>>> El día 15/07/2013 a las 13:30, Duncan Murdoch
escribió:
> On 13-07-15 5:27 AM, Helios de Rosario wrote:
>> Hi,
>>
>> What should I do if I want to use "non-default" R-options for the
>> functions inside a package, but not affect the options of the rest
of
>> the session? Specifically, I'm thin
> Hi all,
>
> several packages print messages during loading. How do I avoid to
see
> them when the packages are in the defaultPackages?
>
> Here is an example.
>
> With this in ~/.Rprofile
> ,[ ~/.Rprofile ]
> | old <- getOption("defaultPackages")
> | options(defaultPackages = c(old, "file
Hi Helios,
"Helios de Rosario" writes:
>> Hi all,
>>
>> several packages print messages during loading. How do I avoid to
> see
>> them when the packages are in the defaultPackages?
>>
>> Here is an example.
>>
>> With this in ~/.Rprofile
>> ,[ ~/.Rprofile ]
>> | old <- getOption("defaul
With minor corrections, the original problem can be "solved" with nlxb
from nlmrt package.
> coef(modeln)
a b
-0.8470857 409.5190808
ssquares = 145585533
but since
> svd(modeln$jacobian)$d
[1] 5.128345e+04 6.049076e-14
>
I may have made nlmrt too robust.
JN
On 13-07
Hi
I am trying to use ANOVA for mixture experiment.
The data can be seen below
Run Blend X1 X2 X3 Response
1 Pure 0.0 0.0 1.0 43.6
2 Pure 0.0 0.0 1.0 42.0
3 Pure 0.0 0.0 1.0 43.0
4 Binary 0.0 0.5 0.5 30.0
5 Binary 0.0 0.5 0.5 29.4
6 Binary 0.0 0.5 0.5 33.6
7 Pure 0.0 1.0 0.0 17.6
8 Pure
Arun, thanks for clearing that up. I was making assumptions about R that I
shouldnât have when it comes to column(variable) assignments. Hope everyone
has a good Monday. BNC
From: arun kirshna [via R] [mailto:ml-node+s789695n4671473...@n4.nabble.com]
Sent: Friday, July 12, 2013 11:59 PM
To:
On 15/07/2013 8:49 AM, Andreas Leha wrote:
Hi Helios,
"Helios de Rosario" writes:
>> Hi all,
>>
>> several packages print messages during loading. How do I avoid to
> see
>> them when the packages are in the defaultPackages?
>>
>> Here is an example.
>>
>> With this in ~/.Rprofile
>> ,[ ~
Great! Thank you, Nello. It works on my larger matrices.
On Mon, Jul 15, 2013 at 2:40 AM, Blaser Nello wrote:
> For small matrix you could use a for-loop.
>
> for (i in 1:nrow(randomized)){
> randomized[i,randomized[i,]!=0] <- sample(original[i,original[i,]!=0])
> }
> randomized
>
> If you ha
On 07/15/2013 06:25 AM, Duncan Murdoch wrote:
On 15/07/2013 8:49 AM, Andreas Leha wrote:
Hi Helios,
"Helios de Rosario" writes:
>> Hi all,
>>
>> several packages print messages during loading. How do I avoid to
> see
>> them when the packages are in the defaultPackages?
>>
>> Here is an exam
> But the ANOVA results is incorrect when I use
>
> aov.out = aov(Response~-1+x1*x2*x3).
First, this doesn't work at all. Your variables are in upper case and your data
environment is not specified, so you must have done something else. I can get
an answer using something like
anova(lm(Response~
Looking for a function which returns the stem of the filename given a path.
i.e.
> file.stem("/the/path/to/afile.txt")
> afile
regards
--
Witold Eryk Wolski
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://sta
Dear Petr,
Sorry for the delay. I've been out.
Unfortunately, your code doesn't work either even when using fromLast = T.
Thank you for your help and your time.
Santi
>
> From: PIKAL Petr
>To: Santiago Guallar
>Cc: r-help
>Sent: Wednesday, July 10, 2013 8:
Read the Posting Guide... use plain text email.
help.search("file")
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#. ##.#. Live Go...
Hello,
You can use ?basename to write a file.stem function:
basename("/the/path/to/afile.txt")
file.stem <- function(x){
bn <- basename(x)
gsub("\\..*$", "", bn)
}
file.stem("/the/path/to/afile.txt")
Hope this helps,
Rui Barradas
Em 15-07-2013 15:23, Witold E Wolski escrev
Hello everyone on the list
I want to run a FindU function on as many files as in a directory. Due of my
limitations on data management I am calling the function, wait until finish the
process, call for the next file and so on. I have almost 120 CSV files on the
same directory. This the example
Rui, Thank you
On 15 July 2013 16:32, Rui Barradas wrote:
> Hello,
>
> You can use ?basename to write a file.stem function:
>
>
> basename("/the/path/to/afile.txt")
>
> file.stem <- function(x){
> bn <- basename(x)
> gsub("\\..*$", "", bn)
> }
> file.stem("/the/path/to/afile.txt")
This should also work:
> mask <- randomized > 0
> mask
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] FALSE FALSE TRUE FALSE FALSE TRUE
[2,] TRUE TRUE FALSE TRUE TRUE TRUE
[3,] FALSE TRUE TRUE FALSE FALSE TRUE
[4,] TRUE TRUE FALSE TRUE FALSE TRUE
> randomized[mask] <- original[mask]
> rand
Hm, so you probably ask to get something what you actually do not want.
AFAIK what I called âfinalâ is the same as you asked for with your toy
data except of column âactâ which you can easily get rid of.
If what I suggested does not work with your real data you shall prepare better
exa
Dear all,
I would like to do t-test across two columns "name" with different "cat"
with overall mean ("val").
(Removing if there is a single observation)
And finally, make a matrix with t-value and p-value associated with a name
(in rows) and cat (in columns)
dput(x)
structure(list(name = struc
Dear all,
I am having problem on plots in R for some cases.
For example:
plot( 1:10 )
text( 1:10, letters[1:10], cex = 1)
works well but :
plot(1:10)
text(1:10, letters[1:10], cex = 0.9)
returns :
Erreur dans text.default(1:10, letters[1:10], cex = 0.9) :
impossible de charger la polic
Bill,
Thanks so much for the tips. For the project that I'm doing now, I don't
need to worry about bullet-proofing it too much, but this is definitely
educational and will help tons when the next one comes around.
Cheers,
eric
On Sun, Jul 7, 2013 at 10:22 AM, William Dunlap wrote:
> > x.lab <
I was not clear on exactly what you want to do. I assume that you want to
automate the process given 120 file names. Here is part of a script that
will parse the file names and give you the output file you want. If it is
something different, then clarify your request:
> input <- "C:/RHTests/dat
On 15/07/2013 15:59, Alice Julien-Laferrière wrote:
Dear all,
I am having problem on plots in R for some cases.
Your machine has Cairo: Cairo-based X11 devices do not use those fonts.
So it seems you have done something you have not shown us.
The fix is to use the Cairo-based X11 devices.
It is a (very) large marked point process object that is being split (see
original post), and that is some other kind of structure. So I think all
bets are off, in terms of what kind of duration to expect.
The folks on r-sig-geo might be able to help with suggestions for faster
approaches.
In ter
Hi,
Not sure about the format for the 2nd part.
df1<- ##data
library(plyr)
df2<-ddply(df1,.(name,cat),summarize,
cbind(t.test(val,df1$val)$statistic,t.test(val,df1$val)$p.value))
df3<-cbind(df2[,1:2],data.frame(df2[,3]))
colnames(df3)[3:4]<- c("t-val","p.val")
library(reshape2)
df3m<- melt(df3
Dear all,
Is there a possibility to remove a geom from a ggplot? Background suppose I
have a function which returns a ggplot object after some data re-formatting and
aggregation. While this ggplot object is fine in 90% of the cases it turns out
that for some cases I want to suppress one of the
This may be close to what you want:
> t.val <- by(x, x$cat, function(y) if (min(table(y$name)>1)) {
+ t.test(val~name, y)})
> t.out <- do.call(rbind, sapply(t.val, function(y) c(y$statistic,
+ p.value=y$p.value)))
> t.out
t p.value
p178266580 -0.1156475 0.91440
On Mon, Jul 15, 2013 at 1:08 PM, Thaler,Thorn,LAUSANNE,Applied
Mathematics wrote:
> Dear all,
>
> Is there a possibility to remove a geom from a ggplot? Background suppose I
> have a function which returns a ggplot object after some data re-formatting
> and aggregation. While this ggplot object
On 2013-07-15, at 4:00 AM, r-help-requ...@r-project.org wrote:
> On 7/14/2013 3:05 PM, Rolf Turner wrote:
>> On 15/07/13 08:57, Spencer Graves wrote:
>>
>>
>>>
>>> You may know that "mixed effects" is another term for "repeated
>>> measurements".
>>
>>
>>
>> I must of course preface t
On Jul 15, 2013, at 4:11 AM, Andreas Leha wrote:
> Hi all,
>
> several packages print messages during loading. How do I avoid to see
> them when the packages are in the defaultPackages?
I'm pretty sure this has been asked before (although the wrinkle of adding a
package to default packages ma
Hello again,
Let say I have following data-frame:
MyData <- data.frame(Names1 = paste("XXX", 1:150), Values1 = 1:150 + 10,
Values2 = 1:150)
Now I want to plot this data-frame with googleVis. Therefore I run
following codes:
library(googleVis)
plot(gvisBarChart(MyData, xvar="Names1", yvar=c("Val
Hi,
May be I misunderstood ur question.
The output David got could be also obtained by:
#df1 dataset
library(plyr)
df2<-ddply(df1,.(cat),function(x) if(min(table(x$name))>1){x1<-
t.test(val~name,x);cbind(t=x1$statistic,p.value=x1$p.value)})
df2
# cat t p.value
#1 p17826658
The "plantbreeding" package is on R-forge. See this link for install
instructions:
http://rplantbreeding.blogspot.com/
Kevin
On Sun, Jul 14, 2013 at 12:49 AM, waqas shafqat wrote:
> sir i could not find the plant breeding libraray in Rgui3.0.0
>
> [[alternative HTML version deleted]]
Hi,
May be this also works.
basename(file_path_sans_ext("/the/path/to/afile.txt"))
#[1] "afile"
A.K.
- Original Message -
From: Rui Barradas
To: Witold E Wolski
Cc: r-help@r-project.org
Sent: Monday, July 15, 2013 10:32 AM
Subject: Re: [R] file.stem?
Hello,
You can use ?basename to
HI Irucka,
May be this is what you wanted:
pat<-paste(paste0("http://www.",siter[,1],"..rdb";),collapse="|")
pat
[1]
"http://www.02437100..rdb|http://www.02439500..rdb|http://www.02441500..rdb|http://www.02446500..rdb|http://www.02467000..rdb|
---
#or
pat1<-paste("http://www.";, siter[,1], "..rdb", sep="", collapse = "|")
identical(pat,pat1)
#[1] TRUE
A.K.
- Original Message -
From: arun
To: Irucka Embry
Cc: R help
Sent: Monday, July 15, 2013 2:47 PM
Subject: Re: converting numeric to character and using character pattern
Hello,
I use glm within a function testing for the appearence of the coexistence
of (minor allels in a subset of) snps. And then I extract the
Pr(>|z|)-value for the interaction. Principally it works but sometimes the
function stops because this "value for the interaction" is NA. For
instance, t
I think what you want is
?try ##or
?tryCatch
## The second is more flexible but slightly more complicated.
to trap the error and perhaps refit the model without interaction?
Cheers,
Bert
On Mon, Jul 15, 2013 at 10:45 AM, Hermann Norpois wrote:
> Hello,
>
> I use glm within a function testing
Hi David,
David Winsemius writes:
> On Jul 15, 2013, at 4:11 AM, Andreas Leha wrote:
>
>> Hi all,
>>
>> several packages print messages during loading. How do I avoid to see
>> them when the packages are in the defaultPackages?
>
> I'm pretty sure this has been asked before (although the wrink
Hi Martin,
Martin Morgan writes:
> On 07/15/2013 06:25 AM, Duncan Murdoch wrote:
>> On 15/07/2013 8:49 AM, Andreas Leha wrote:
>>> Hi Helios,
>>>
>>> "Helios de Rosario" writes:
>>>
>>> >> Hi all,
>>> >>
>>> >> several packages print messages during loading. How do I avoid to
>>> > see
>>> >>
Dear Jari,
Thanks a lot for your reply.
Let me describe more about my dataset: it consists of count data for samples,
and shows how many copies of a certain bacteria are in a sample. There are 92
samples (rows), and about 3500 kinds of bacteria (columns). Variable 1 (Var1)
is a factor indicati
Hellos,
I am looking for a package that can scrub values from flow cytometry data
file keywords that have either the "$" or "@" keyword prefix. Are their any
options out there? The options from WEHI, TreeStar, and whoever made "File
Sanitizer" aren't working.
thanks for any help on this,
--
- s
Dear R-Users,
I need a very fast and reliable database solution so I try to serialize a
data.frame (to binary data) and to store this data to an SQLite database.
This is what I tried to do:
library(RSQLite)
con <- dbDriver("SQLite")
db <- dbConnect(con, "test")
dbSendQuery(db, 'CREATE TABLE fr
Group - I'm having problems with the 'cluster' package. Installation
appears successful but attempts to load it with either library() or
require() result in the error message
Error in library(cluster) : there is no package called ‘cluster’
All that appears to be installed is cluster.dll in the
On 15.07.2013 23:51, David Stevens wrote:
Group - I'm having problems with the 'cluster' package. Installation
appears successful but attempts to load it with either library() or
require() result in the error message
Error in library(cluster) : there is no package called ‘cluster’
All that ap
dear R experts: I am curious again about R memory allocation strategies.
Consider an intentionally inefficient program:
ranmatme <- function( lx, rx ) {
m <- matrix(NA, nrow=lx, ncol=rx)
for (li in 1:rx) {
cat("\tLag i=", li, "object size=", object.size(m), "\n")
m[,li]
Dear R users,
I'm currently using the Graphics package to display several hundred of
matrix objects, using a layout and the image() function.
It works well except for large matrices (> 1000*1000) or for a large
number of matrices (there is a limitation around 400 if I remember well)
To solve th
I can give you the answer to #1. If you had put a "print(str(m))" you
would have seen that initially the matrix was setup as logical which
requires 4 bytes per element. On the first assignment of a numeric, the
mode of 'm' is changed to numeric which requires 8 bytes per element; that
is the reas
Hi,
I am currently working with the package JGL.
I tried to run JGL function "JGL(Y,penalty="fused",...)" over my data, but
there is the following error which makes it stop:
"Error in while ((iter==0)|| (iter < maxiter && diff_value>tol))
missing value where TRUE/FALSE needed. "
I would appreciat
I could be wrong, but I would guess that doing what you are describing is very
unusual. Most of the time the data frame is mapped to a table in the database
so the rows can be searched. Storing data frames as BLOBs really seems odd.
Note that there is an R-sig-db mailing list for questions of th
I have a data frame like shown below
sample1 sample2 sample3 sample4 sample5 a P P I P P b P A P P A c P P P
P P d P P P P P e M P M A P f P P P P P g P P P A P h P P P P P
I want to keep only those rows which have all "P" across all the columns.
Since the matrix is large (about 20,000
Hi,
If I understand it correctly,
df1<- read.table(text="
sample1 sample2 sample3 sample4 sample5
a P P I P P
b P A P P A
c P P P P P
d P P P P P
e M P M A P
f P P P P P
g P P P A P
h P P P P P
",sep="",header=TRUE,stringsAsFactors=FALSE)
df1[rowSums(df1=="P")==ncol(df1),]
# sample1 samp
thx, jim. makes perfect sense now.
I guess a logical in R has a few million possible values ;-).
(Joke. I realize that 4 bytes is to keep the code convenient and faster.)
regards,
/iaw
Ivo Welch (ivo.we...@gmail.com)
On Mon, Jul 15, 2013 at 4:26 PM, jim holtman wrote:
> I can give you
You mentioned data.frame at one place and matrix at another. Matrix would be
faster.
#Speed comparison
#Speed
set.seed(1454)
dfTest<- as.data.frame(matrix(sample(LETTERS[15:18],5*1e6,replace=TRUE),ncol=5))
system.time(res<-dfTest[rowSums(dfTest=="P")==ncol(dfTest),])
# user system elapsed
#
There was a small typo in my sendmail command. I was passing header=header
instead of headers=header. This argument was accepted because of `...`
argument which is not the last argument as one would except. A minor mistake
caused me a bit of frustration. This issue is resolved now. Thanks everyone.
I am currently working with the package JGL.
I tried to run JGL function "JGL(Y,penalty="fused",...)" over my data, but
there is the following error which makes it stop:
"Error in while ((iter==0)|| (iter < maxiter && diff_value>tol))
missing value where TRUE/FALSE needed. "
Any special idea for
Please do not multiple-post.
What did the maintainer say when you asked him (see the posting guide)?
On 16/07/2013 05:36, Azam Peyvand wrote:
I am currently working with the package JGL.
I tried to run JGL function "JGL(Y,penalty="fused",...)" over my data, but
there is the following error whic
68 matches
Mail list logo
