[R] New methods for generic functions show and print : some visible with ls(), some not

2010-02-26 Thread Joris Meys
en by ls(). On the other hand, the method "print" does! So when I use rm(list=ls()) I still have the "initialize" and "show" method, but the "print" method is gone. Am I forgetting something somewhere? It's rather inconvenient to have to run the definition fil

Re: [R] New methods for generic functions show and print : some visible with ls(), some not

2010-02-27 Thread Joris Meys
Thank you both for your answers. On Fri, Feb 26, 2010 at 7:58 PM, Duncan Murdoch wrote: > > You aren't seeing the print method, you are seeing a newly created print > generic function. As Uwe mentioned, print() is not an S4 generic, so when > you create your print method, a new S4 generic also g

[R] weighted.median function from package R.basic

2010-03-30 Thread Joris Meys
ction in the package limma too, but I didn't use that before. Anybody who knows what happened to R.basic? Cheers Joris -- Joris Meys Statistical Consultant Ghent University Faculty of Bioscience Engineering Department of Applied mathematics, biometrics and process control Coupure Links

Re: [R] list index rules evaluation behavior

2010-03-30 Thread Joris Meys
p1745398p1745398.html > Sent from the R help mailing list archive at Nabble.com. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/pos

Re: [R] Create a new variable

2010-03-30 Thread Joris Meys
mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical Consultant Ghent University Faculty of Biosci

[R] compare multiple values with vector and return vector

2010-04-02 Thread Joris Meys
pretty huge datasets. Anybody an idea? thank you in advance. Cheers Joris -- Joris Meys Statistical Consultant Ghent University Faculty of Bioscience Engineering Department of Applied mathematics, biometrics and process control Coupure Links 653 B-9000 Gen

Re: [R] compare multiple values with vector and return vector

2010-04-02 Thread Joris Meys
010 at 12:55 AM, wrote: > Hi Joris > > Try Is.element function: is.element (x,y) > > Regards > mohamed > > Joris Meys a écrit : > > Dear all, >> >> I have a vector, and for each element I want to check whether it is equal >> to >> any eleme

Re: [R] Co-occurrence

2010-04-05 Thread Joris Meys
ps://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical Consultant Ghent University Faculty of Bioscience Engineering Depa

[R] use of random and nested factors in lme

2010-04-05 Thread Joris Meys
clear to me on this subject. What confuses me most, is the use of the intercept in the random factor. Does this mean the intercept is seen as random, has a random component or is it just notation? In different mails from this list I found different explanations. Thank you in advance. Cheers Joris

Re: [R] Error in pf(q, df1, df2, lower.tail, log.p) : Non-numeric argument to mathematical function

2010-04-26 Thread Joris Meys
m. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible

Re: [R] Tinn-R

2010-04-26 Thread Joris Meys
rovide commented, minimal, self-contained, reproducible code. > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > htt

[R] Overview of error and warning messages for teaching + examples

2009-09-30 Thread joris meys
n R. Off course this mailing list is an extensive source of all kind or errors and examples, but I was wondering if there is a somewhat more structured overview. All help appreciated. Kind regards Joris Meys __ R-help@r-project.org mailing list https

Re: [R] Overview of error and warning messages for teaching + examples

2009-09-30 Thread joris meys
ne number: 972-52-7275845 > FaceBook: Tal Galili > My Blogs: > http://www.talgalili.com (Web and general, Hebrew) > http://www.biostatistics.co.il (Statistics, Hebrew) > http://www.r-statistics.com/ (Statistics,R, English) > > > > > On Wed, Sep 30, 2009 at 2:50 PM, jor

Re: [R] Overview of error and warning messages for teaching + examples

2009-10-01 Thread joris meys
Hi Ben, other members of the list also pointed me to this possibility. I'll certainly add the manual. Kind regards Joris On Thu, Oct 1, 2009 at 4:14 AM, Ben Bolker wrote: > Tal Galili gmail.com> writes: > >> >> Hi Joris, >> >> Good luck with your work. >> >> I know how to set up a wiki. But I w

Re: [R] Minimum and blanks

2009-10-01 Thread joris meys
The minimum of a character vector is returned according to the unicode or ascii code values, but I guess that's not what Chris is looking for. Chris, please always try to provide a minimal working example of code to see what is wrong. Something goes wrong with the input in your example. As Peter

Re: [R] Missing functions

2009-10-01 Thread joris meys
r: package 'psych' does not have a name space > > is there a way to clean up and reinstall the packages? > > thanks, > Kenny > > On Thu, Oct 1, 2009 at 4:02 PM, joris meys wrote: >> >> Hi, >> >> It installed perfectly for me, both packages

Re: [R] GAM question

2009-10-01 Thread joris meys
Hi Daniel, I presume you mean coefficients by values. Try ?gamObject and read the info in the help file. A GAM object (returned by the function gam() ) contains an element coefficients that is easily accessible. That one contains all estimated coefficients, including those for the smooth functions

Re: [R] confint fails in quasibinomial glm: dims do not match

2009-10-02 Thread joris meys
Confint doesn't work if you have a multi-dimensional dependent variable. Kind regards Joris On Fri, Oct 2, 2009 at 4:29 AM, smith_cc wrote: > > I am unable to calculate confidence intervals for the slope estimate in a > quasibinomial glm using confint(). Below is the output and the package info

[R] Extracting year from a date object

2009-10-06 Thread joris meys
Hi all, this one left me a bit puzzled, as I don't seem to find a function to perform this easily. I must have overlooked the obvious, so sorry in advance. I have a list of dates in numerical format (i.e. 34576), defined as the number of days that passed since january 1st 1900. So I apply the fun

Re: [R] Unable to load 'doBy' package

2009-10-06 Thread joris meys
Hi Lauren, from the error message it looks like you have a problem with the package "Hmisc". doBy is dependent on that one. Can you check whether Hmisc is installed in your R version? For me, doBy loads without problems. On Tue, Oct 6, 2009 at 12:18 PM, Lauren Szathmary wrote: > Hi all, > > I am

Re: [R] Plot

2009-10-06 Thread joris meys
Actually, it's a bit more complex than just plot(). First you have to make Days an ordered factor, otherwise you get indeed that error. days <- c("Mon", "Tue", "Wed", "Thu", "Fri", "Sat","Sun") Days <- factor(days,levels=days,ordered=T) Then you want to make a plot with 2 axes, you should check :

Re: [R] foreach loop - rejection method

2009-10-08 Thread joris meys
First, you should define your function as : test <- function(cand2,phi,lambda, whatever-arguments-you-want-to-use-further){... insert code here ...} All variables you use inside a function only exist within that function. Your parameters/arguments is the interface between the function and the outs

Re: [R] I can not install DAAG package . help

2009-10-08 Thread joris meys
Quite a strange error. The function installed.packages should be in the utils package. Try ?installed.packages and see if you get a help file. If it doesn't find that function, there is a problem with your R installation. It might be something went wrong with downloading the package itself. DAAG i

Re: [R] Need a vectorized way to avoid two nested FOR loops

2009-10-08 Thread joris meys
Neat piece of code, Jim, but it still uses a nested loop. If you order the matrix first, you only need one passage through the whole matrix to find the information you need. Off course I don't take into account the ordering. If the ordering algorithm doesn't work in linear time, then it doesn't re

Re: [R] Display more than one plot

2009-10-10 Thread joris meys
Hi Emkay, If you want to look at different plots together, you can also plot them side by side in the same plot window. You can specify this using for example: par(mfcol=c(2,2)) ( see ?par and check mfrow and mfcol) or layout(matrix(1:4,2,2)) (see ?layout and ?matrix) eg : x <- c

Re: [R] Display more than one plot

2009-10-10 Thread joris meys
which graph where. It's pretty straight forward. Kind regards Joris On Sun, Oct 11, 2009 at 12:50 AM, joris meys wrote: > Hi Emkay, > > If you want to look at different plots together, you can also plot > them side by side in the same plot window. > > You can sp

Re: [R] Error: cannot allocate vector of size 1.2 Gb

2009-10-12 Thread joris meys
Dear Roman, could you give us the trace given by traceback() ? I suspect the error is resulting from the permutations and/or jackknife procedure in the underlying functions specaccum and specpool. You can take a look at the package R.huge, but that one is deprecated already. There are other packa

Re: [R] Error: cannot allocate vector of size 1.2 Gb

2009-10-12 Thread joris meys
Thank you for the tips as well (R.huge looks promising)! > >> traceback() > 7: vector("integer", length) > 6: integer(nbins) > 5: tabulate(bin, pd) > 4: as.vector(data) > 3: array(tabulate(bin, pd), dims, dimnames = dn) > 2: table(y[, factor]) > 1: accumcomp(

Re: [R] is that possible to graph 4 dimention plot

2009-10-12 Thread joris meys
I'm basically put off by the question itself. Plotting a 4-dimensional graph is rather "complicated" if the world has only 3 dimensions. A 4-dimensional representation is typically a movie (with time as the 4th dimension). You could try to project a heatmap on a 3D surface graph, but I doubt this w

Re: [R] Error: cannot allocate vector of size 1.2 Gb

2009-10-12 Thread joris meys
gards Joris -- Forwarded message ------ From: romunov Date: Mon, Oct 12, 2009 at 3:14 PM Subject: Re: [R] Error: cannot allocate vector of size 1.2 Gb To: joris meys I hope this workspace file attached is what you were looking for. Cheers, Roman On Mon, Oct 12, 2009 at 2:07 PM,

Re: [R] R commander.

2009-10-12 Thread joris meys
As the error says, you have different row numbers in your variables. The variable $Chromosome has no values. try : ann <- data.frame( ann [-3] ) Cheers Joris On Mon, Oct 12, 2009 at 8:29 AM, Ilyas . wrote: > i have two RData files,,i want to print them to check the format of the > tables in the

Re: [R] SPSS long variable names

2009-10-13 Thread joris meys
Sorry to be so blunt, but I cannot believe PSPP can't save a dataset as a .csv file for example. That should be the prefered format to transport a dataset to any other statistical package, including R. csv files are universal. Cheers Joris On Tue, Oct 13, 2009 at 3:30 PM, Orvalho Augusto wrote:

Re: [R] SPSS long variable names

2009-10-13 Thread joris meys
, 2009 at 4:35 PM, joris meys wrote: > Sorry to be so blunt, but I cannot believe PSPP can't save a dataset > as a .csv file for example. That should be the prefered format to > transport a dataset to any other statistical package, including R. csv > files are universal. > > Ch

Re: [R] Request Infomation

2009-10-13 Thread joris meys
The easy way out is to send the file containing only the function to the people willing to use it. They can load the function from the file "foo.R" using : source ("path/foo.R") The hard way out is to write a package. For that, see : http://cran.r-project.org/doc/manuals/R-exts.html Cheers Joris

Re: [R] Function to find prime numbers

2009-10-13 Thread joris meys
On Tue, Oct 13, 2009 at 5:30 PM, Barry Rowlingson wrote: > On Tue, Oct 13, 2009 at 2:41 PM, Thomas Lumley > wrote: >> On Tue, 13 Oct 2009, AJ83 wrote: >> >>> >>> I need to create a function to find all the prime numbers in an array. Can >>> anyone point me in the right direction? > >  This almos

Re: [R] histogram

2009-10-13 Thread joris meys
see ?hist On Tue, Oct 13, 2009 at 1:06 PM, Dmitry Gospodaryov wrote: > Dear R developers,How I can build a histogram from matrix: > > 0 0.5 1 > > 0.25 34 43 65 > 1 23 35 54 > 4 22 29 42 > 10 21 22 29 > 20 15 17 20 > > (first string is represented names of columns, > first column is represented na

Re: [R] histogram

2009-10-13 Thread joris meys
On second thought, You'll have to reorganize your data first, and then check the function histogram() in the package lattice. cheers On Wed, Oct 14, 2009 at 12:43 AM, joris meys wrote: > see ?hist > > On Tue, Oct 13, 2009 at 1:06 PM, Dmitry Gospodaryov > wrote: >> Dear

Re: [R] Selecting initial numerals

2009-10-14 Thread joris meys
On Tue, Oct 13, 2009 at 6:48 PM, PDXRugger wrote: > > II just want to create a new object with the first two numerals of the data. > Not sure why this isnt working, consider the following: > > EmpEst$naics=c(238321, 624410, 484121 ,238911, 81, 531110, 621399, > 541613, > 524210 ,236115 ,811121

Re: [R] Selecting initial numerals

2009-10-14 Thread joris meys
Josh, One way would be to convert the numeric vector to a character and use the function substr(). Following code returns a numeric vector with the 2 first digits of every element. naics=c(238321, 624410, 484121 ,238911, 81, 531110, 621399,541613,524210 ,236115 ,811121 ,236115 ,236115 ,621610

Re: [R] SPSS long variable names

2009-10-14 Thread joris meys
On Wed, Oct 14, 2009 at 2:45 PM, Robert Baer wrote: >> The problem is the limit of 8 characters long on variable >> names. > > And again, my answer is that one approach would be to map SHORT names to > long variable LABELS.  This was a common use of labels before variable names > supporte

Re: [R] Time Dependent Cox Model

2009-10-14 Thread joris meys
Well, it might be wise to elaborate a bit more about the variables and what exactly you want e.g. death-time to be. I'd interprete it as time of death, but the fact that it is 0/1, means it is a logical (?) binary variable of some sort. Please ask your question in such a way that somebody who doe

Re: [R] SPSS long variable names

2009-10-14 Thread joris meys
but is not beautfull and uses > more resources. > > The great thing would be to support long variables names on read.spss. > > Thanks guys for everything > Caveman > > On Wed, Oct 14, 2009 at 4:52 PM, joris meys wrote: >> Hi Orvalho, >> >> question : wher

Re: [R] Time Dependent Cox Model

2009-10-14 Thread joris meys
On Wed, Oct 14, 2009 at 9:19 PM, quaildoc wrote: > > Some suggested that go into more detail on what I wanted to accomplish and > the rest of my code.  I want to accomplish exactly what Fox did in this > article( http://www.nabble.com/file/p25897307/appendix-cox-regression.pdf > appendix-cox-regre

Re: [R] help: combine multiple 'pairs plots'

2009-10-15 Thread joris meys
Don't know if it will work, but did you try ?layout ? Cheers Joris On Thu, Oct 15, 2009 at 2:39 PM, Jesse Poland wrote: > I would like to combine multiple pairs plots (each one being a matrix of > correlation plots) into a single graphic. I have tried the par() function, > which works well to co

Re: [R] help: combine multiple 'pairs plots'

2009-10-15 Thread joris meys
Don't know if it will work, but did you try ?layout ? Otherwise you can take a look at ?split.screen too. Cheers Joris On Thu, Oct 15, 2009 at 2:39 PM, Jesse Poland wrote: > I would like to combine multiple pairs plots (each one being a matrix of > correlation plots) into a single graphic. I hav

Re: [R] Solving an equation in R (and using a loop)

2009-10-15 Thread joris meys
I tried : > test <- 0:1 > test2 <- f(0:1) > plot(test,test2) > plot(test,log(test2)) and > test <- seq(0,1,by=0.1) > test2 <- sapply(test,f) > plot(test,log(test2)) > plot(test,test2) with the values you gave. And according to this result, you did nothing wrong. Zero is the onl

Re: [R] loop and plot

2009-10-19 Thread joris meys
Hi Rene, the problem is probably due to the fact that R will send all plots to the same graphical output window. Each next plot just replaces the previous one. if it's only a few plots, you can divide the graphical window with the commands par(mfrow=...) (see ?par) or layout(matrix(...)) (see ?la

Re: [R] loop and plot

2009-10-19 Thread joris meys
,3,1,2,3,1),colour=c("a","b","c","c","a","b","b") > )) > kk = function(f) >             { >              ls=as.character(f) >              for (i in length(f)) >              { >              pie(dataset.ta

Re: [R] Convert data frame entries to numbers

2009-10-19 Thread joris meys
2 problems : test seems to be a data frame or list with one variable. So you have to specify : as.numeric(test$classcol) But this will make the internal factor levels the real values, not the numbers you specified. What you need, is as.numeric(as.character(test$classcol)) Cheers Joris On Mon, O

Re: [R] how to get rid of 2 for-loops and optimize runtime

2009-10-19 Thread joris meys
Hi Ian, first of all, take a look at the functions sapply, mapply, lapply, tapply, ... : they are the more efficient way of implementing loops. Second, could you elaborate a bit further on the data set : the amount of the month ago, is that one value from another row, or the sum of all values in

[R] Labels in horizontal dendrogram not placed correctly?

2009-11-16 Thread joris meys
Hi all, I tried plotting a horizontal dendrogram, but it seems as if the labels are not taken into account in the function plot.dendrogram(). A minimal example : Test <- data.frame( x1x = c(1:10), x2x = c(2:11), x3x = c(11:2) ) TestDist <- daisy(data.frame(t(Test))) TestA

[R] Fwd: dendrogram

2009-11-16 Thread joris meys
I was missing something. Thx Dennis. -- Forwarded message -- From: Dennis Murphy Date: Tue, Nov 17, 2009 at 12:34 AM Subject: Re: dendrogram To: jorism...@gmail.com Hi, There are a couple of things you could do to pull the text back into the dendrogram plot, both of which I fo

[R] categorisation of continuous variables in R

2009-11-23 Thread joris meys
Dear all, I'm looking for a function comparable to switch, to categorize a continuous variable in a few levels. Off course that can be done with a series of ifelse statements, but that looks rather clumsy. I looked at switch, but couldn't figure out how to use it for this. I guess that's not possi

Re: [R] categorisation of continuous variables in R

2009-11-23 Thread joris meys
Never mind, found the function : cut(test,breaks=c(0,10,50,90,100),labels=c("lowest","low","high","highest"),include.lowest=T,right=F) Cheers Joris On Mon, Nov 23, 2009 at 2:14 PM, joris meys wrote: > Dear all, > > I'm looking for a funct

Re: [R] Help: Beanplots calculating wrong average

2009-11-24 Thread joris meys
Hi, The log of the mean is not the same as the mean of the logs, that's a no-brainer. Guess you use the beanplot from the package with the same name. > beanplot(As1988,log="") gives the correct plot. Next time, could you provide a minimal code example we can run ourselves? If we don't know what

Re: [R] Re adline()

2009-11-24 Thread joris meys
I can't reproduce the error you have. With me, it runs all fine. In fact, readline apparently flushes automatically, as I didn't have to use the flush.console() at all. test.r --- cat("1- 24horas\n") cat("2- 12horas\n") cat("3- 8horas\n") selection<-readline(prompt="\nSelecciona numero de h

Re: [R] Labels in horizontal dendrogram not placed correctly?

2009-11-24 Thread joris meys
; On Mon, Nov 16, 2009 at 07:28, joris meys wrote: >> Hi all, >> >> I tried plotting a horizontal dendrogram, but it seems as if the >> labels are not taken into account in the function plot.dendrogram(). >> >> A minimal example : >> Test <- data.f

Re: [R] Re adline()

2009-11-25 Thread joris meys
I run your script from the windows console (evidently), and as said, I can't reproduce your error. Can you give me the exact command you use to run it from the windows console? Cheers Joris On Tue, Nov 24, 2009 at 7:33 PM, yonosoyelmejor wrote: > > Sure,but my problem is that the script is run f

Re: [R] Arrhenius Plot 2 with lattice

2009-11-26 Thread joris meys
Please provide minimal and self-contained code. Without being able to read in the dataset, it's impossible to check the code. Cheers Joris On Thu, Nov 26, 2009 at 9:13 AM, Markus Häge wrote: > Hello, > > sorry for incomplete code... > > with this I read the file and calculate my stuff. I have a

Re: [R] R help with princomp and pam clustering

2009-11-26 Thread joris meys
On Thu, Nov 26, 2009 at 1:04 AM, Tyler82 wrote: > > Hi all! > I am working with R package cluster and I have a little problem: > let's say I have two datasets...first one ("A") is divided into 4 clusters > by means of Pam algorythm. > Let's say I want to project the second database ("B") onto the

Re: [R] R help with princomp and pam clustering

2009-12-03 Thread Joris Meys
Hi Tyler, sorry, I missed your response. Don't know if it's solved already, but some remarks. the predict uses a princomp object, so your second command should work. I guess the problem is either the data format of B (should contain exactly the same amount of columns, with the same names), or the

[R] S4 Classes, nested objects and references

2009-12-03 Thread Joris Meys
Hi all, I'm currently programming my first complete package in S4. (thanks to Christophe Genolini for the nice introduction he wrote). I have an object "Data" with a number of slots. One of those slots is "meteo". Now "Meteo" is on itself a class with again a number of slots (like rainfall, temper

Re: [R] Formatting of numbers on y axis

2009-12-03 Thread Joris Meys
The problem is easily solved by plotting salaries$salary/100 Cheers Joris On Thu, Dec 3, 2009 at 5:16 PM, Wells Oliver wrote: > Hello all. I have the following: > > plot(salaries$yearID, salaries$salary, type='n', xaxt='n', xlab='', > yaxt='n', ylab='') > axis(1, at=unique(salaries$yearID),

[R] modelling probabilities instead of binary data with logistic regression

2009-03-24 Thread joris meys
Dear all, I have a dataset where I reduced the dimensionality, and now I have a response variable with probabilities/proportions between 0 and 1. I wanted to do a logistic regression on those, but the function glm refuses to do that with non-integer values in the response. I also tried lrm, but th

Re: [R] modelling probabilities instead of binary data with logistic regression

2009-03-24 Thread joris meys
t; More seriously, you could transform the proportions to logits > logit <- log(p/(1-p)) > and fit a linear regression. > > Kjetil > > On Tue, Mar 24, 2009 at 3:30 PM, joris meys wrote: > >> Dear all, >> >> I have a dataset where I reduced the dimensionality, a

Re: [R] modelling probabilities instead of binary data with logisticregression

2009-03-25 Thread joris meys
e plural of anecdote is not data. > ~ Roger Brinner > > The combination of some data and an aching desire for an answer does not > ensure that a reasonable answer can be extracted from a given body of > data. > ~ John Tukey > > -Oorspronkelijk bericht- > Van: r-help-bou

Re: [R] Hopefully a simple question

2009-03-25 Thread joris meys
Dear Jacy, As Jorge said. To clarify : what you ask, is actually to compare the first element of genotype with the first element of FREQ, the second of genotype with the second of FREQ and so on. If the lengths differ, the shorter vector is recycled. This toy example illustrates that principle :

[R] metafor: using mixed models

2010-01-04 Thread Joris Meys
Dear all, I'm currently applying a mixed model approach to meta analysis using the package metafor. I use the "model.matrix()" function to create dummy variables. The option btt gives me the combined test for the dummies. Problem is, I don't know which indices I have to use, and can't really figur

Re: [R] Tinn-R user guide (latex sources) available on GitHub

2013-11-27 Thread Joris Meys
ing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consultant Ghent University Faculty of Bioscience Eng

Re: [R] Should there be an R-beginners list?

2013-11-27 Thread Joris Meys
you completely. > > cheers, > > Rolf Turner > > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/ > posting-guide.html

Re: [R] Get a list of installed commands

2010-06-13 Thread Joris Meys
> R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consul

Re: [R] How to output text to sink from data frame line by line without column names

2010-06-13 Thread Joris Meys
at I would like in test.txt  is: >> >> "I will be including text of various sorts in this file so cannot use >> print table or similar command" >> col1 >> col2 >> col3 >>  a  234 >>  B   44 >>  C  566 >> >> Many th

Re: [R] Clustering algorithms don't find obvious clusters

2010-06-13 Thread Joris Meys
; >>> Henrik >>> >>>        [[alternative HTML version deleted]] >>> >>> __ >>> R-help@r-project.org mailing list >>> https://stat.ethz.ch/mailman/listinfo/r-help >>> PLEASE do read the posting guide

Re: [R] Large Data

2010-06-14 Thread Joris Meys
And this one is only from last week. Please, read the posting guides carefully. Cheers Joris -- Forwarded message -- From: Joris Meys Date: Sat, Jun 5, 2010 at 11:04 PM Subject: Re: [R] What is the largest in memory data object you've worked with in R? To: Nathan Stephen

Re: [R] Large Data

2010-06-14 Thread Joris Meys
e at Nabble.com. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, repro

Re: [R] logistic regression with 50 varaibales

2010-06-14 Thread Joris Meys
> > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducib

Re: [R] Html help

2010-06-14 Thread Joris Meys
> R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consul

Re: [R] how to ignore rows missing arguments of a function when creating a function?

2010-06-14 Thread Joris Meys
; 4950  0    NA 1 1    9    90 > 4951 NA    NA     0 1   15   151 > > It works perfectly if I do > cl2(testdata, lm(income ~ transport + reg_f ), female) > > but not for cl2(testdata, lm(dist ~ transport + reg_f ), female) > Or any other case ab

Re: [R] lme, spline

2010-06-15 Thread Joris Meys
ASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consultant Ghent University Faculty of Bioscience Engineering Department of Applied mathematics, biome

Re: [R] Prime Numbers Pkgs

2010-06-15 Thread Joris Meys
_ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Jo

Re: [R] R software

2010-06-16 Thread Joris Meys
gt; > > -- > Gulbanu > >        [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-gui

Re: [R] generating samples from multivariate distributions

2010-06-16 Thread Joris Meys
; R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consultant Ghent Uni

Re: [R] Question

2010-06-16 Thread Joris Meys
; > > > >        [[alternative HTML version deleted]] > > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and

Re: [R] error message stops "for" command

2010-06-16 Thread Joris Meys
ing-guide.html > and provide commented, minimal, self-contained, reproducible code. > > -- Joris Meys Statistical consultant Ghent University Faculty of Bioscience Engineering Department of Applied mathematics, biometrics and process control tel : +32 9 264 59 87 joris

Re: [R] Multiple plots in a single page and stripplot()

2010-06-17 Thread Joris Meys
mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consultant Ghent University Faculty of Bioscience

Re: [R] Subtraction of group means using AGGREGATE and MERGE

2010-06-17 Thread Joris Meys
> Many thanks, > > Ben Cocker > MSc Statistics at UCL, London, UK > > ______ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/post

Re: [R] big big problem

2010-06-17 Thread Joris Meys
   EC02   9 10 > > M2        2007    EC02    6 10 > > > > I put 1 instead of 0 in the code > > series=ddply(X,.(variable,CLUSTER),transform,series=rev(value2-c(0,cumsum(rev(value1[-1]) > > but

Re: [R] Plotting different symbols in R for different values in a vector

2010-06-17 Thread Joris Meys
__ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Jori

Re: [R] simulating data from a multivariate dist

2010-06-17 Thread Joris Meys
amp; Regards, > Suman Dhara > >        [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/po

Re: [R] Problems using allEffects() (package effect)

2010-06-17 Thread Joris Meys
_ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consultant Ghen

Re: [R] R licensing query

2010-06-17 Thread Joris Meys
Worst case scenario you can install R as a user as well, you don't need administration rights. Regarding data analysis in Excel : http://www.pages.drexel.edu/~bdm25/excel2007.pdf : peer reviewed Other sources : http://www.coventry.ac.uk/ec/~nhunt/pottel.pdf http://www.forecastingprinciples.com/f

Re: [R] RKWARD

2010-06-17 Thread Joris Meys
l > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consultant Ghent University Faculty of Bioscience Engineering Department of Applied mathematics, biometrics and process control tel : +32 9 264 59 87 joris.m...@ugent.be ---

Re: [R] Problems using allEffects() (package effect)

2010-06-18 Thread Joris Meys
ects(mod=lm1, xlevels=list(temps=c(1:3), Traitement=1:3)) > > And here is the error I get: > Error in apply(mod.matrix[, components], 1, prod) : index out of range > > In the example given in the help page there are factor variables and > numeric variables so I don't think the proble

Re: [R] Problems using allEffects() (package effect)

2010-06-18 Thread Joris Meys
; > But the data is not at all balanced : >> table(data$Série, data$Traitemen) > >     1  2  3 >  7  0 28 31 >  8 24 17  0 > > In batch 7 there is no mouse treated with treatment 1, and in batch 8 > there is no mouse treated with treatment 3. Maybe the error

Re: [R] double integral

2010-06-18 Thread Joris Meys
ead the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consultant Ghent University Faculty of Bioscience Engineering Department of Applied mathemati

Re: [R] Extract estimates from each dataset: MI package

2010-06-18 Thread Joris Meys
stinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Joris Meys Statistical consultant Ghent University Faculty of Bioscience Engineering Department of Applied mathemat

Re: [R] Reading data from xls..........please help

2010-06-18 Thread Joris Meys
se it hasn't found perl. >> Do you have perl on your system? Do you need to specify the path to >> perl, as in the examples for Windows in help(read.xls)? >> >> Barry >> > >        [[alternative HTML version deleted]] > > ___

Re: [R] Extract estimates from each dataset: MI package

2010-06-18 Thread Joris Meys
ot work and i keep getting an error "object of type 'closure' is not > subsettable" > > Hope this makes my question clear. > > Thanks > > > -- > View this message in context: > http://r.789695.n4.nabble.com/Extract-estimates-from-each-dataset-MI-package-tp2259864p2260019.html >

Re: [R] plotting radial dendrograms

2010-06-18 Thread Joris Meys
? > > Regards, > > Carson > > -- > Carson J. Q. Farmer > ISSP Doctoral Fellow > National Centre for Geocomputation > National University of Ireland, Maynooth, > http://www.carsonfarmer.com/ > > __ > R-help@r-projec

Re: [R] mgcv, testing gamm vs lme, which degrees of freedom?

2010-06-18 Thread Joris Meys
t; > ## much higher edf according to summary: > summary(f2$gam) > > --- > > ______ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.

  1   2   3   4   5   >