Re: [R] Entropy based feature selection in R

2011-08-13 Thread andy1234
Hello everyone, Any thoughts in this one please? The only thing I found was the FSelector package (http://en.wikibooks.org/wiki/Data_Mining_Algorithms_In_R/Dimensionality_Reduction/Feature_Selection#Aviable_Feature_Ranking_Techniques_in_FSelector_Package). Unfortunately though it seems to be

Re: [R] Automating an R function call

2011-08-13 Thread RobertJK
Thanks for your help everyone! I'm happy enough with an asynchronous solution here. Thanks! Robert -- View this message in context: http://r.789695.n4.nabble.com/Automating-an-R-function-call-tp3740070p3740333.html Sent from the R help mailing list archive at Nabble.com. [[alternative

[R] Adjacency Matrix help

2011-08-13 Thread collegegurl69
I have created an adjacency matrix but have not been able to figure something out. I need to put zeros on the diagonal of the adjacency matrix. For instance, location (i,i) to equal 0. Please help. Thanks -- View this message in context:

[R] Finding an average time spent

2011-08-13 Thread erinbspace
Hello R help! I am extremely new to R (as in 3 just days) and I've been using it to do some pretty basic things. I am frustratingly stuck on one point, and am so so so close to figuring it out, but just far enough away to ask for some (perhaps embarrassingly easy) help. I have a dataset,

Re: [R] sapply to bind columns, with repeat?

2011-08-13 Thread Katrina Bennett
Hi Weidong Gu, This works! For my clarity, and so I can repeat this process if need be: The 'mat' generates a matrix using whatever is supplied to x (i.e. coop.dat) using the columns from position 9:length(x) of 6 columns (by row). The 'rem.col' generates a matrix of the first 1:8 columns of 8

[R] post

2011-08-13 Thread bdeepthi
Hello, I was trying to plot multiple graph using par(mfrow=c(3,2)). But this is giving me the following error: Error in axis(side = side, at = at, labels = labels, ...) : X11 font -adobe-helvetica-%s-%s-*-*-%d-*-*-*-*-*-*-*, face 1 at size 8 could not be loaded Could someone decode this

Re: [R] Adjacency Matrix help

2011-08-13 Thread Timothy Bates
diag(adjMatrix) -0 On Aug 13, 2011, at 7:34 AM, collegegurl69 wrote: I have created an adjacency matrix but have not been able to figure something out. I need to put zeros on the diagonal of the adjacency matrix. For instance, location (i,i) to equal 0. Please help. Thanks

Re: [R] Finding an average time spent

2011-08-13 Thread Gabor Grothendieck
On Fri, Aug 12, 2011 at 4:23 PM, erinbspace erin.brasw...@gmail.com wrote: Hello R help! I am extremely new to R (as in 3 just days) and I've been using it to do some pretty basic things. I am frustratingly stuck on one point, and am so so so close to figuring it out, but just far enough away

Re: [R] Adjacency Matrix help

2011-08-13 Thread collegegurl69
Thanks so much for your quick reply. it seems to work. the problem is that it now places actual zeros on the diagonal whereas the rest of the adjacency matrix has dots to represent zeroes. Do you have any ideas on how to change these zeros to dots like in the rest of the adj matrix? Or is it the

[R] optimization problems

2011-08-13 Thread Kathie
Dear R users I am trying to use OPTIMX(OPTIM) for nonlinear optimization. There is no error in my code but the results are so weird (see below). When I ran via OPTIM, the results are that Initial values are that theta0 = 0.6 1.6 0.6 1.6 0.7. (In fact true vales are 0.5,1.0,0.8,1.2, 0.6.)

Re: [R] optimization problems

2011-08-13 Thread Kathie
To be honest, The first derivative of my objective function is very complicated so I ignore this. Could it lead to this sort of problem? Kathie -- View this message in context: http://r.789695.n4.nabble.com/optimization-problems-tp3741005p3741010.html Sent from the R help mailing list archive

Re: [R] Any alternatives to draw.colorkey from lattice package?

2011-08-13 Thread Jim Lemon
On 08/13/2011 04:34 AM, Mikhail Titov wrote: Hello! I’d like to have a continuous color bar on my lattice xyplot with colors lets say from topo.colors such that it has ticks labels at few specific points only. Right now I use do.breaks level.colors with somewhat large number of steps. The

Re: [R] Finding an average time spent

2011-08-13 Thread Jim Lemon
On 08/13/2011 06:23 AM, erinbspace wrote: Hello R help! I am extremely new to R (as in 3 just days) and I've been using it to do some pretty basic things. I am frustratingly stuck on one point, and am so so so close to figuring it out, but just far enough away to ask for some (perhaps

Re: [R] Any alternatives to draw.colorkey from lattice package?

2011-08-13 Thread Felix Andrews
You can just specify the label positions, you don't need to give labels for every color change point: (there is an 'at' for the color changes and a 'labels$at' for the labels) levelplot(rnorm(100) ~ x * y, expand.grid(x = 1:10, y = 1:10), colorkey = list(at = seq(-3,3,length=100),

[R] linear regression

2011-08-13 Thread maggy yan
dear R users, my data looks like this PM10 Ref UZ JZ WT RH FT WR 1 10.973195 4.338874 nein Winter Dienstag ja nein West 26.381684 2.250446 nein SommerSonntag nein ja Süd 3 62.586512 66.304869 ja SommerSonntag nein nein Ost 4

Re: [R] Passing on groups argument to xyplot within a plotting function

2011-08-13 Thread Felix Andrews
The problem is that xyplot tries to evaluate 'groups' in 'data' or in the formula environment. Your local function environment (where the variable named groups is defined) is neither of these. There are a couple of ways to get the evaluation to work out; here is one: pb - list(F1 = 1:8, F2 = 1:8,

[R] Excluding NAs from round correlation

2011-08-13 Thread Julie
Hello, I am quite new to R and I am trying to get a round correlation from a table with dozens of columns. However, all the columns contain several blank places which show to me as NAs. Then, when I type round(cor(data),2), I get no results - everything (except correlation of one column with the

Re: [R] linear regression

2011-08-13 Thread Steven Kennedy
your dataframe needs to be called Nord. If it is not, then replace Nord with the actual name of your dataframe On Sat, Aug 13, 2011 at 10:43 PM, maggy yan kiot...@googlemail.com wrote: dear R users, my data looks like this         PM10       Ref   UZ     JZ         WT   RH   FT   WR 1  

Re: [R] Excluding NAs from round correlation

2011-08-13 Thread Weidong Gu
check ?cor Please note the parameter 'use' Weidong Gu On Sat, Aug 13, 2011 at 9:06 AM, Julie julie.novak...@gmail.com wrote: Hello, I am quite new to R and I am trying to get a round correlation from a table with dozens of columns. However, all the columns contain several blank places

Re: [R] post

2011-08-13 Thread Uwe Ligges
On 13.08.2011 06:52, bdeep...@ibab.ac.in wrote: Hello, I was trying to plot multiple graph using par(mfrow=c(3,2)). But this is giving me the following error: Error in axis(side = side, at = at, labels = labels, ...) : X11 font -adobe-helvetica-%s-%s-*-*-%d-*-*-*-*-*-*-*, face 1 at size 8

Re: [R] optimization problems

2011-08-13 Thread John C Nash
optimx with BFGS uses optim, so you actually incur some overhead unnecessarily. And BFGS really needs good gradients (as does Rvmmin and Rcgmin which are updated BFGS and CG, but all in R and with bounds or box constraints). From the Hessian, your function is (one of the many!) that have pretty

[R] optimization problems

2011-08-13 Thread Ravi Varadhan
Kathie, It is very difficult to help without adequate information. What does your objective function look like? Are you maximizing (in which case you have to make sure that the sign of the objective function is correct) or minimizing? Can you try optimx with the control option

Re: [R] define variables from a matrix

2011-08-13 Thread Dennis Murphy
There may well be more efficient ways to do this, but here's one attempt: foo - function(x, val) if(any(x == val, na.rm = TRUE)) which(x == val) else NA u - apply(A, 1, function(x) foo(x, 20L)) v - apply(A, 1, function(x) foo(x, 100L)) ifelse(u v, v, NA) [1] 3 5 NA NA NA HTH, Dennis On Fri,

Re: [R] linear regression

2011-08-13 Thread Dennis Murphy
Hi: Try something like this, using dat as the name of your data frame: xyplot(PM10 ~ Ref | WR, data = dat, type = c('p', 'r')) The plot looks silly with the data snippet you provided, but should hopefully look more sensible with the complete data. The code creates a four panel plot, one per

Re: [R] define variables from a matrix

2011-08-13 Thread David Winsemius
On Aug 12, 2011, at 7:18 PM, gallon li wrote: I have a following matrix and wish to define a variable based the variable A=matrix(0,5,5) A[1,]=c(30,20,100,120,90) A[2,]=c(40,30,20,50,100) A[3,]=c(50,50,40,30,30) A[4,]=c(30,20,40,50,50) A[5,]=c(30,50,NA,NA,100) A [,1] [,2] [,3] [,4]

Re: [R] Any alternatives to draw.colorkey from lattice package?

2011-08-13 Thread Mikhail Titov
Felix: Thank you! Perhaps I should read documentation more careful as I missed that another `at`. lattice latticeExtra are so marvelous so I hardly want to use anything else. Mikhail On 08/13/2011 07:31 AM, Felix Andrews wrote: You can just specify the label positions, you don't need to give

[R] R's handling of high dimensional data

2011-08-13 Thread andy1234
Hello all, I am looking at doing text classification on very high dimensional data (about 300,000 or more features) and upto 2000 documents. I am quite new to R though, and was just wondering if R and it's libraries would scale to such high dimensions. Any thoughts will be much appreciated.

[R] degrees of freedom does not appear in the summary lmer :(

2011-08-13 Thread xy
Hi , Could someone pls help me about this topic, I dont know how can i extract them from my model!! Thanks, Sophie -- View this message in context: http://r.789695.n4.nabble.com/degrees-of-freedom-does-not-appear-in-the-summary-lmer-tp3741327p3741327.html Sent from the R help mailing list

[R] Plotting and quantiles

2011-08-13 Thread Mark D.
Dear R users, This is most likely very basic question but I am new to R and would really appreciate some tips on those two problems. 1) I need to plot variables from a data frame. Because of some few high numbers my graph is really strange looking. How could I plot a fraction of the samples

[R] Individual p-values for correlation matrices

2011-08-13 Thread January Weiner
Dear all, I am calculating each-against-each correlations for a number of variables, in order to use the correlations as distances. This is easy enough using just cor(), but not if I want to have a p-value for each calculated correlation, and especially if I want to correct them for multiple

[R] Own R function doubt

2011-08-13 Thread garciap
Hi to all the people again, I was writting a simply function in R, and wish to collect the results in a excel file. The work goes as follows, Ciervos-function(K1, K0, A, R,M,Pi,Hembras) {B-(K1-K0)/A T1-(R*Pi*Hembras-M*Pi+B)/(Pi-M*Pi+R*Pi*Hembras) P1-Pi-B R1-P1*Hembras*R M1-P1*M

Re: [R] fit.mult.impute() in Hmisc

2011-08-13 Thread Paul Johnson
On Thu, Mar 31, 2011 at 2:56 PM, Yuelin Li li...@mskcc.org wrote: I tried multiple imputation with aregImpute() and fit.mult.impute() in Hmisc 3.8-3 (June 2010) and R-2.12.1. The warning message below suggests that summary(f) of fit.mult.impute() would only use the last imputed data set.

Re: [R] fit.mult.impute() in Hmisc

2011-08-13 Thread Frank Harrell
For your approach how do you know that either summary or vcov used multiple imputation? You are using a non-rms fitting function so be careful. Compare with using the lrm fitting function. Also repace Design with the rms package. Please omit confidentiality notices from your e-mails. Frank

[R] seeking advice about rounding error and %%

2011-08-13 Thread Paul Johnson
A client came into our consulting center with some data that had been damaged by somebody who opened it in MS Excel. The columns were supposed to be integer valued, 0 through 5, but some of the values were mysteriously damaged. There were scores like 1.18329322 and such in there. Until he tracks

Re: [R] Own R function doubt

2011-08-13 Thread R. Michael Weylandt
It sounds like the data frame produced by Simulate() doesn't set the names you want. You can probably fix this by including colnames(Gestion) c(B,T1,... # etc) immediately after the simulation. Can't confirm this without knowing which of the excel/R interface packages you're using, but I'd be

Re: [R] degrees of freedom does not appear in the summary lmer :(

2011-08-13 Thread Joshua Wiley
Hi Sophie, It is not clear what the degrees of freedom should be in an lmer model, so their not appearing is intentional. There is fairly extensive discussion of this topic in the archives for the R-sig-mixed list. See, for example: http://rwiki.sciviews.org/doku.php?id=guides:lmer-tests

Re: [R] linear regression

2011-08-13 Thread Mikhail Titov
Don't forget to load `lattice` package. `latticeExtra` with `panel.ablineeq` can be also helpful. This was however for plotting. For subset regression by each WR without plotting you'd use something like `lapply` or `sapply`. ans - sapply(unique(data$WR), function(dir) { out -

Re: [R] Plotting and quantiles

2011-08-13 Thread R. Michael Weylandt
I believe you received an informative answer to both these questions from Daniel Maiter one hour and twenty five minutes after sending your question: I repeat it here just in case you didn't get it. -- Q1 is very opaque because you are not even saying what kind of plot

Re: [R] efficient use of lm over a matrix vs. using apply over rows

2011-08-13 Thread darius
Good Bless you Duncan. Your explanation is crisp and to the point. Thank you. -- View this message in context: http://r.789695.n4.nabble.com/efficient-use-of-lm-over-a-matrix-vs-using-apply-over-rows-tp870810p3742043.html Sent from the R help mailing list archive at Nabble.com.

[R] NA in lm model summary, no NA in data table!

2011-08-13 Thread stefy
Dear users of R! I have problems with linear model summary. I do not have any NA values in my data table, but in summary of linear model there are some NA instead of resulst. I don´t know why :( I am interested in ecological factors influencing temperature in ant nest, I have data concerning ant

Re: [R] Excluding NAs from round correlation

2011-08-13 Thread Julie
Thank you, I found this in the help pack: use: an optional character string giving a method for computing covariances in the presence of missing values. This must be (an abbreviation of) one of the strings everything, all.obs, complete.obs, na.or.complete, or pairwise.complete.obs I should

[R] what is Inverse link functions in linear modelling of location

2011-08-13 Thread Laura
Hello I have a problem with the following function (http://www.oga-lab.net/RGM2/func.php?rd_id=ismev:gev.fit): gev.fit(xdat, ydat = NULL, mul = NULL, sigl = NULL, shl = NULL, mulink = identity, siglink = identity, shlink = identity, muinit = NULL, siginit = NULL, shinit = NULL,

Re: [R] Excluding NAs from round correlation

2011-08-13 Thread Julie
The help pack says: use: an optional character string giving a method for computing covariances in the presence of missing values. This must be (an abbreviation of) one of the strings everything, all.obs, complete.obs, na.or.complete, or pairwise.complete.obs If I used everything, the results

Re: [R] degrees of freedom does not appear in the summary lmer :(

2011-08-13 Thread Dennis Murphy
Hi: This is worth reading and bookmarking: http://glmm.wikidot.com/faq HTH, Dennis On Sat, Aug 13, 2011 at 6:31 AM, xy wtemptat...@hotmail.co.uk wrote: Hi , Could someone pls help me about this topic, I dont know how can i extract them from my  model!! Thanks, Sophie -- View this

Re: [R] efficient use of lm over a matrix vs. using apply over rows

2011-08-13 Thread darius
Good Bless you Duncan. Your explanation is crisp and to the point. Thank you. -- View this message in context: http://r.789695.n4.nabble.com/efficient-use-of-lm-over-a-matrix-vs-using-apply-over-rows-tp870810p3742058.html Sent from the R help mailing list archive at Nabble.com.

[R] Casualty Actuarial Society request for proposals for R Workshop

2011-08-13 Thread Kevin Burke
I'm a property-casualty actuary, use R in at my job, and lurk on the list. In conjunction with one of its meetings, the Casualty Actuarial Society (I'm a member) is looking for proposals from people to teach a workshop in R and I thought members of the list might be interested. I've pasted the

Re: [R] Adjacency Matrix help

2011-08-13 Thread Daniel Nordlund
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of collegegurl69 Sent: Saturday, August 13, 2011 1:01 AM To: r-help@r-project.org Subject: Re: [R] Adjacency Matrix help Thanks so much for your quick reply. it seems to work.

[R] compiling r from source on Windows 7 (64 bit)

2011-08-13 Thread Erin Hodgess
Dear R People: Hope you're having a nice Saturday. I'm trying to compile R-2.13.1 from source on Windows 7 (64 bit). I've been able to compile on a 32 bit without any problems. I changed my BINPREF64, WIN, DEFS_W64 in MkRules.local and did the usual stuff with the jpeg, etc. But things are

Re: [R] compiling r from source on Windows 7 (64 bit)

2011-08-13 Thread Joshua Wiley
Erin, You can build packages without compiling from source (what made you think you couldn't?). Did you make sure when you installed the Rtools (I am assuming you are using those rather than going out and getting everything you need on your own) that you included everything for 64 bit builds?

Re: [R] seeking advice about rounding error and %%

2011-08-13 Thread Joshua Wiley
Hi Paul, What about using: x[x != as.integer(x)] - NA I cannot think of a situation off hand where this would fail to turn every non integer to missing. I wonder if there is really a point to this? Can the client proceed with data analysis with any degree of confidence when an unknown

[R] [R\ Compiling R from source on Windows 7 (64 bit) solved

2011-08-13 Thread Erin Hodgess
Hello again. Due to the excellent help from Josh Wiley, I ran back in the C:/R directory with only changing WIN = 64 in the MkRules.local file (other than the JPEG, etc). All was well. Thanks, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences

[R] Using get() or similar function to access more than one element in a vector

2011-08-13 Thread Joseph Sorell
Dear R-users, I've written a script that produces a frequency table for a group of texts. The table has a total frequency for each word type and individual frequency counts for each of the files. (I have not included the code for creating the column headers.) Below is a sample: Word Total

Re: [R] seeking advice about rounding error and %%

2011-08-13 Thread Ken
How about something like: If(round(x)!=x){zap} not exactly working code but might help Ken On Aug 13, 2554 BE, at 3:42 PM, Paul Johnson pauljoh...@gmail.com wrote: A client came into our consulting center with some data that had been damaged by somebody who opened it in MS Excel. The

Re: [R] Using get() or similar function to access more than one element in a vector

2011-08-13 Thread Joshua Wiley
Hi Joseph, Without a reproducible example, you probably will not get the precise code for a solution but look at ?list Rather than doing what you are doing now, put everything into a list, and then you will not need to use get() at all. You will just work with the whole list. It can take a bit