Re: [R] How to define new matrix based on an elementary row operation in a single step?

Is this sufficiently single steppish for you? D <- A <- matrix(1:16, 4) D[3, ] <- 2 * D[1, ] + D[3, ] # Alternately, you could do this # but it is much messier, and I do not see how # two steps is really an issue # you want to end up with two matrices anyways # so it's not like you save memory by

Re: [R] Looking for an image (R 64-bit on Linux 64-bit) on Amazon EC2

Thanks! But how could I find out their names on Amazon EC2? -- View this message in context: http://r.789695.n4.nabble.com/Looking-for-an-image-R-64-bit-on-Linux-64-bit-on-Amazon-EC2-tp2338938p2357731.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] How to define new matrix based on an elementary row operation in a single step?

Sorry for possible misunderstanding: I want to define a matrix (B) based on an existing matrix (A) in a single step and keep A unchanged: > #Existing matrix > A=matrix(1:16,ncol=4) > A [,1] [,2] [,3] [,4] [1,]159 13 [2,]26 10 14 [3,]37 11 15 [4,]4

[R] Calculating p and q values with R

Hi, I have a huge dataset (53 million records). I have to calculate the p and q values of my data. How can I do it in R or perl? I have downloaded R (I'm completely new to R). and the package qvalue but I don't understand how can I call/use qvalue package with R. When I type library(qvalue), it gi

Re: [R] expression() and plot title

Try this: > LinePlot(1,1) > LinePlot=function(a,b){ + # a = slope + # b = y intercept + x=seq(-10,10,0.4) + y=a*x+b + plot(x,y, type="l") + title(paste("a=",a,"b=",b)) + } > > #test > LinePlot(a=-2,b=9) HTH -- View this message in context: http://r.789695.n4.nabble.com/expression-and-plot-t

Re: [R] expression() and plot title

Here is a suggestion you may care to develop func <- function(a, b) { plot(1:10) title(main = bquote(a == .(a)*','~ b == .(b))) invisible() } try with func(1,2) func(36, 2^10) &c -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-proj

Re: [R] expression() and plot title

On Aug 28, 2010, at 12:53 AM, Sancar Adali wrote: What I want to do is put the arguments I supply to a function into the title of a plot Say I'm calling func.1 func.1(a=4,b=4) plot(,..., title("a=4, b=4")) If I'm calling func.1 with different arguments, I want the plot title to refl

Re: [R] About plot graphs

Hi Gavin, Thanks for your advice and the examples explaining plotting settings. The steps on your examples work on my test. > 2) Don't attach() things, you are asking for trouble > If a function has a formula method (which plot does) then use it like > this: plot(Draft_No. ~ Day_of_year, data

[R] expression() and plot title

What I want to do is put the arguments I supply to a function into the title of a plot Say I'm calling func.1 func.1(a=4,b=4) plot(,..., title("a=4, b=4")) If I'm calling func.1 with different arguments, I want the plot title to reflect that. A small detail is that func.1 might have an argumen

Re: [R] How to define new matrix based on an elementary row operation in a single step?

On Aug 28, 2010, at 12:29 AM, Cuckovic Paik wrote: Thanks for respose. you still used two steps to get the new matrix tst2: For Pete"s sake! Can't you see that I didn't _need_ to use tst2. The same operations would have worked on tst. I was just using tst2 because I wanted to be able to

Re: [R] How to define new matrix based on an elementary row operation in a single step?

Thanks for respose. you still used two steps to get the new matrix tst2: step 1: tst2 = tst step 2: perform the row operation in tst2. Can you do this in a single step? A similar example: > tst [,1] [,2] [,3] [,4] [1,]159 13 [2,]26 10 14 [3,]37 11

Re: [R] How to define new matrix based on an elementary row operation in a single step?

On Aug 27, 2010, at 9:21 PM, Cuckovic Paik wrote: I'm not absolutely sure I know what you mean by elementary row operations ... you are supposed to offer test cases and specify the desired output on r-help to support the aging faculties of the helpeRs in this case ... having only the v

Re: [R] Scatterplot question

Is this what you want? x=rnorm(10) y1=0.5*rnorm(10)+0.5*rexp(10) y2=0.5*rnorm(10)+0.5*runif(10) xlimit=range(x) ylimit=range(c(y1,y2)) plot(x,y1,pch="1",col=2,xlim=xlimit, ylim=ylimit, ylab="Y", xlab="X", cex=0.7) points(x,y2,pch="2", col=4,cex=0.7) -- View this message in context: http://r.7

Re: [R] R.matlab package help

Henrik, The line before that is: Received cmd: 2 save C:\Users\FAN'S~1\AppData\Local\Temp\tpe2b4012b_f9ed_402d_af0f_f21ebd8116a6.mat -V6 B Thanks, Michael [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

Re: [R] R.matlab package help

On Fri, Aug 27, 2010 at 6:16 PM, michael wrote: > Henrik, >          When I set remote=FALSE, I got this error: > > ??? Error: A MATLAB string constant is not terminated properly. > Error in ==> MatlabServer at 197 >      eval(expr); Also the lines before. /Henrik > > In R: > >   Retrieving var

[R] How to define new matrix based on an elementary row operation in a single step?

-- View this message in context: http://r.789695.n4.nabble.com/How-to-define-new-matrix-based-on-an-elementary-row-operation-in-a-single-step-tp2341768p2341768.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mail

Re: [R] R.matlab package help

Henrik, When I set remote=FALSE, I got this error: ??? Error: A MATLAB string constant is not terminated properly. Error in ==> MatlabServer at 197 eval(expr); In R: Retrieving variables from the Matlab server: 'B' Sending expression on the Matlab server to be evaluated...:

Re: [R] Problems loading rjags in R 2.11.0 under OS X 10.5.8

dear Professor Ripley, thank you very much for your prompt reply. You are absolutely right. I didn't have the latest version installed although I used the command: install.packages("rjags"). Apparently, you only get the latest version of packages if you have the most up to date version of R, whi

Re: [R] predict.loess and NA/NaN values

From: Philipp Pagel > > In a current project, I am fitting loess models to subsets of data in > order to use the loess predicitons for normalization (similar to what > is done in many microarray analyses). While working on this I ran into > a problem when I tried to predict from the loess models a

Re: [R] R.matlab package help

Hi, let's stick with the remote=FALSE case first. What output does Matlab produce now? That is really key information in order to solve this. Also, make sure that you really restart the MatlabServer. /H On Fri, Aug 27, 2010 at 3:56 PM, michael wrote: > Henrik, >          I updated the MatlabS

Re: [R] R-help

I ran some comparisons of `lp.assign' and `solve_LSAP' on a couple of moderately large problems. The `solve_LSAP' appears to be substantially faster than `lp.assign'. Best, Ravi. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Ravi

Re: [R] R.matlab package help

Henrik, I updated the MatlabServer.m from the link you sent to me, but it still doen't work. I have the same problems as before, when remote=FALSE, getVariable doesn't work; when remote=TRUE, setVariable says: Received cmd: 5 Will read MAT file structure of length: 984 bytes. ??? Undefin

Re: [R] R.matlab package help

://www.braju.com/R/patches/R.matlab/20100827/MatlabServer.m Let me know if it works. If so, I'll submit a new version of R.matlab to CRAN with the fix. /Henrik On Fri, Aug 27, 2010 at 12:27 PM, michael wrote: > Henrik, >         Here is the matlab response: > >>> Matlab

Re: [R] LSAP and the Hungarian algorithm [was: R-help]

Ravi Varadhan jhmi.edu> writes: > > However, The "clue" package has the solve_LSAP() function (as pointed out by > Hans Werner) that solves the linear sum assignment problem, but it uses the > Hungarian algorithm that was requested by you. > > The lp.assign() function in "lpSolve" package (as p

Re: [R] make.rm

On Aug 27, 2010, at 2:38 PM, Bruce Johnson wrote: Where do I find this function? Search? (How else?) rhelpSearch <- function(string, restrict = c("Rhelp10", "Rhelp08", "Rhelp02", "functions" ), matchesPerPage = 100, ...) { RSiteSearch(string=string, restrict = res

Re: [R] make.rm

On 08/27/2010 08:38 PM, Bruce Johnson wrote: > Where do I find this function? On Google, first hit... Now really, If you can't find it, how would you expect anyone else to, without any indication of what the function is supposed to do nor from where you got the idea to look for it? -- Peter Dal

Re: [R] Scatterplot question

On Aug 27, 2010, at 2:16 PM, ashz wrote: Hi, I have two XY datasets (e.g., longitude and concentrations) who share the same X scale. How can I make a simple scatterplot which will combine them both with different colors for the two Y groups? (plot and xyplot solutions are fine with me

Re: [R] interpreting date-related error message

Thanks to everyone for the explanations. Toby On 27 Aug 2010 at 12:46, Phil Spector wrote: > Toby - > Since dat$doy is just a number, the default S3 method > for format is used, where the second argument is the trim > parameter. I suspect you are confusing format (which is for > output) wi

Re: [R] 05 - Re: R and MySQL

Hello De-Jian Thanks. That helped. It works now. On Fri, Aug 27, 2010 at 2:30 AM, De-Jian,ZHAO [via R] < ml-node+2340662-1760624908-138...@n4.nabble.com > wrote: > Waitlist reason: > subscription-bounces+138...@n4.nabble.comis > not on your > Guest List |

[R] make.rm

Where do I find this function? [[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 provide comme

[R] Scatterplot question

Hi, I have two XY datasets (e.g., longitude and concentrations) who share the same X scale. How can I make a simple scatterplot which will combine them both with different colors for the two Y groups? (plot and xyplot solutions are fine with me). Thanks in advance. Cheers -- View this messa

[R] how to use glm for a time course data set

Dear list, I used glm to analyze the effect of treatment, measurement-methods etc. on lesion size in a time course experiments. Could anyone points me out how to extract detailed comparisons of interest? For example: Measurement.T2.M6 vs. Measurement.T2.Baseline, Or even Measurement.T2.M6.Trea

Re: [R] How to plot an expression-label with variable text

On Aug 27, 2010, at 2:51 PM, Dieter Menne wrote: Thanks to both of you. I noted that my example was over-simplified. Looks like I need to correct the environment when nested in a function, but I have to catch the last bus now. Dieter plotExp = function(what) { plot.new() lab =express

Re: [R] Grouping sets of data, performing function and re-assigning values

Hi Johnny, Something like this rbind(NA, dat.med)[as.numeric(dat$image.group), ] should do the trick (with the data you provided and Ista's code). The key is that dat.med has a different row for each level of the factor image.group (and in the same order). The idea is to convert the factor cre

Re: [R] How to calc ratios base on current and previous row?

For this particular case there is a nice shortcut (there is still looping , but it is internal and quick): a <- c(2,2,4, sample(1:10, 97, TRUE) ) b <- cumsum( (1/2)^(99:0)*a )/( (1/2)^(99:0) ) ## compare bb <- numeric(100) bb[1] <- a[1] for( i in 2:100 ) { bb[i] <- 1/2 * bb[i-1] + a[i

Re: [R] interpreting date-related error message

Toby - Since dat$doy is just a number, the default S3 method for format is used, where the second argument is the trim parameter. I suspect you are confusing format (which is for output) with strptime (which is for input). For example, strptime(dat$doy,'%j') will assume that the dates

Re: [R] interpreting date-related error message

On Aug 27, 2010, at 3:29 PM, Toby Gass wrote: Hello, helpeRs, I have a vector of numbers from 1-365 (days of the year) that I would like to convert to a date. There are no NA's and no missing values. I did not insert leading zero's for numbers less than 100. Using the syntax: dat$doy.1 <- a

Re: [R] interpreting date-related error message

Toby, What is it that you're trying to accomplish? There seems to be several ideas confused in your post. It sounds like you have input x <- 1:365 and you want to call some function to return dates? Which date should be returned for input 1? January 1, 2010? Your error is because if you giv

[R] interpreting date-related error message

Hello, helpeRs, I have a vector of numbers from 1-365 (days of the year) that I would like to convert to a date. There are no NA's and no missing values. I did not insert leading zero's for numbers less than 100. Using the syntax: dat$doy.1 <- as.numeric(format(dat$doy, "%j" )) I get the fol

Re: [R] R.matlab package help

Henrik, Here is the matlab response: >> MatlabServer Matlab v7.x or higher detected. Saving with option -V6. Added InputStreamByteWrapper to dynamic Java CLASSPATH. -- Matlab server started! -- Trying to open server socket (port )...done. so I

Re: [R] Grouping sets of data, performing function and re-assigning values

HI Ista, Thanks for the help. The 'cut' function seems to do the trick . I'm not sure why you suggested this line of code: > ddply(dat, .(image.group), transform, measure.median = median(Measurement)) I think I might have confused the issue by putting a 'Measurement' column in my example in th

Re: [R] How to calc ratios base on current and previous row?

Hi Marcus, I am guessing you are thinking in terms of Excel, where in column B, you could enter the formula =(A2 + 0.5*B1) and just drag it down however many cells you wanted. In R, it would be expressed a bit differently. If k indexes your rows, I believe you want: B[k] = A[k] + 0.5 * B[k - 1]

Re: [R] R.matlab package help

Hi, could you send me what Matlab is outputting. When I send the following from R: > evaluate(matlab, "A=2;"); Sending expression on the Matlab server to be evaluated...: 'A=2;' and Matlab should print something like: Received cmd: 1 "eval" string: "A=2;" Received an 'OK' reply (0) from the Ma

Re: [R] non-linear plot parameters

Just a small fix to my solution; inserted below. On 2010-08-27 3:51, Peter Ehlers wrote: On 2010-08-26 15:52, Marlin Keith Cox wrote: I agree. I typically do not use non-linear functions, so am seeing the "art" in describing functions of non-linear plots. One last thing. I tried to use a self-s

Re: [R] How to plot an expression-label with variable text

On Aug 27, 2010, at 2:28 PM, baptiste auguie wrote: hi, try this lab =bquote(paste("Estimated ", t[50]," from ",.(what))) bquote doesn't need the paste() if you use plotmath separators: lab =bquote(Estimated ~t[50]~from~.(what)) text(0.5,0.2,lab) HTH, baptiste On 27 August 2010 20:19

Re: [R] Parsing a XML file

Ok. Pardon me. Its is not a bug. My XML is malformed. The program that generates writes things like: >5yr Which make the XML unsable. Sorry everyone Caveman On Wed, Aug 25, 2010 at 5:31 AM, Duncan Temple Lang wrote: > > xmlDoc() is not the function to use to parse a file. > > Use > > doc

Re: [R] How to plot an expression-label with variable text

Thanks to both of you. I noted that my example was over-simplified. Looks like I need to correct the environment when nested in a function, but I have to catch the last bus now. Dieter plotExp = function(what) { plot.new() lab =expression(paste("Estimated ", t[50]," from tgv")) text(0.5,0

Re: [R] How to plot an expression-label with variable text

On Aug 27, 2010, at 2:19 PM, Dieter Menne wrote: Disclaimer: I have read plotmath, but maybe it's too late today: How do I get the two labels to be the same: plot.new() lab =expression(paste("Estimated ", t[50]," from tgv")) text(0.5,0.5,lab) # Should look the same as above. I could not get t

Re: [R] Nestad ANOVA with random Factors

Dear Dennis, I sampled five different plots inside six predetermined areas in three diferent days. Am I doing the analysis right? Thanks, Pablo -- View this message in context: http://r.789695.n4.nabble.com/Nestad-ANOVA-with-random-Factors-tp2340725p2341484.html Sent from the R help mailing

Re: [R] How to plot an expression-label with variable text

hi, try this lab =bquote(paste("Estimated ", t[50]," from ",.(what))) HTH, baptiste On 27 August 2010 20:19, Dieter Menne wrote: > > plot.new() > lab =expression(paste("Estimated ", t[50]," from tgv")) > text(0.5,0.5,lab) > # Should look the same as above. I could not get the substitute righ

Re: [R] step

On Fri, 2010-08-27 at 11:03 -0300, Silvano wrote: > Hi, > > how can I change the significance level in test F to select > variable in step command? > > I used > > step(model0, ~x1+x2+x3+x4, direction=c("forward"), test='F', > alpha=.05) > > but it does't work. Well, there is no 'test' argume

Re: [R] How to obtain seed after generating random number?

It is nearly impossible to go back after the fact and figure out the seed you started with. So, you need to be careful to record the seed first. If you are doing your simulations with a function then it is a good idea to always start in the function by saving a record of the seed and returning it w

[R] How to plot an expression-label with variable text

Disclaimer: I have read plotmath, but maybe it's too late today: How do I get the two labels to be the same: plot.new() lab =expression(paste("Estimated ", t[50]," from tgv")) text(0.5,0.5,lab) # Should look the same as above. I could not get the substitute right: what = "tgv" lab =expression(pas

Re: [R] Importance of levels in a factor variable

Thanks Greg. Actually, we have 5000 levels and it is not an import problem. I looked into combine.levels in the Hmisc package. The problem with this approach is that it takes the frequency of levels, then combines infrequent levels into one level called "Others". If you apply this to the complete d

Re: [R] Grouping sets of data, performing function and re-assigning values

Hi Johnny, If I understand correctly, I think you can use cut() to create a grouping variable, and then calculate your summaries based on that. Something like dat <- read.csv("~/Downloads/exampledata.csv") dat$image.group <- cut(dat$a.ImageNumber, breaks = seq(0, max(dat$a.ImageNumber), by = 3))

Re: [R] R.matlab package help

Hi, what is the version of Matlab you are running? Could you show a minimum verbatim example - starting with a fresh R session and library("R.matlab") - that gives you the error? Have you tried the example of help(Matlab)? Then, for troubleshooting it could be useful to ask Matlab to display va

Re: [R] pairwise correlations for large dataset

Actually, I've answered my own question. It turns out that transposing the expression matrix first, outside of the loop, significantly improves the speed. It now looks like the entire matrix should be calculated in a day or two. So I think this solution should be fine. I now have this: exp

Re: [R] How to calc ratios base on current and previous row?

On 27-Aug-10 17:17:50, Marcus Drescher wrote: > Hi all, > I want to calculate in each row a ratio based on number in the current > row and the previous row. > Is there a way to do this without for-loops because that is extremely > slow. > > A B > [1] 2 2 > [2] 2

Re: [R] Grouping sets of data, performing function and re-assigning values

Hi all, Since I could not attach a file to my original e-mail request, for those who want to look at an example of a data file I am working with, please use this link: http://dl.dropbox.com/u/4637975/exampledata.csv Thanks again, Johnny. __ R-help

Re: [R] How to calc ratios base on current and previous row?

Your question is unclear. If the calculation of B[3] now depends on the newly calculated value of B[2], then the answer is no: you must loop. If B[3] is based on the original B[2] value, then the answer is yes, and the solution is just a matter of simple indexing, which I leave as an exercise. --

Re: [R] multivariate distributions

Hi: One form of a multivariate beta distribution is the Dirichlet, so using package sos, library(sos) # install if necessary u <- findFn('Dirichlet distribution') grepFn('Dirichlet distribution', u, column = 'Description', ignore.case = TRUE) reveals about 25 matches, about half of which l

[R] How to calc ratios base on current and previous row?

Hi all, I want to calculate in each row a ratio based on number in the current row and the previous row. Is there a way to do this without for-loops because that is extremely slow. A B [1] 2 2 [2] 2 3 [3] 4 5,5 ... B2 = A2 + 0.5*B

[R] pairwise correlations for large dataset

Hello, I need to calculate the correlation for all pairwise combinations in a very large matrix. I have 25,000 elements and need to calculate the pairwise correlation with a different set of 5,000 elements. I have written code that works, but it is extremely slow. At the current rate, it wil

[R] OT: Re: Fwd: R SOFTWARE

On 18/08/10 10:46, Nokuzola Simakuhle wrote: > NB: This email and its contents are subject to the Eskom Holdings Limited > EMAIL LEGAL NOTICE > > which can be viewed at http://www.eskom.co.za/email_legalnotice > > Since your email policy allows for no copying or distribution of your message,

Re: [R] R.matlab package help

Henrik, Thanks for your reply. I am using Matlab version 7.9.0(R2009b). This is my first try using R.matlab. I started matlab in the same machine, and I run the MatlabServer.m in matlab, here is my code in R: library(R.matlab) library(MASS) matlab <- Matlab() open(matlab) setVerbose(mat

Re: [R] Error: package/namespace load failed for 'IlluminaHumanMethylation27k.db'

Hi Kristel -- On 08/27/2010 05:48 AM, Kristel van Eijk wrote: > Hello everyone, > > I have a problem when loading the library: IlluminaHumanMethylation27k.db > > I installed the package from zip (from the Bioconductor website), and also > installed all the packages it asked for. It went all succ

[R] step

Hi, how can I change the significance level in test F to select variable in step command? I used step(model0, ~x1+x2+x3+x4, direction=c("forward"), test='F', alpha=.05) but it does't work. -- Silvano Cesar da Costa Departamento de Estatística Universida

[R] Grouping sets of data, performing function and re-assigning values

Hi there, I hope you have time to read this question and offer a suggestion or two. My basic question is this: I have data in sets of three. I would like to combine the data from each set, perform a function (probably just taking the median and MAD), then re-assign these values to each of t

Re: [R] AIC using nls function

John: 1. As always, and as requested (see posting guide), a small reproducible example might help. 2. What is CLi in your model? 3. In general, AIC may not be particularly meaningful as a measure of fit quality penalized for model complexity in NON-linear models unless the different models are

[R] Error: package/namespace load failed for 'IlluminaHumanMethylation27k.db'

Hello everyone, I have a problem when loading the library: IlluminaHumanMethylation27k.db I installed the package from zip (from the Bioconductor website), and also installed all the packages it asked for. It went all succesfull: "package 'IlluminaHumanMethylation27k.db' successfully unpacked an

[R] [R-pkgs] New package: mpt

Dear all, I have submitted a new package called mpt to CRAN. It contains functions for fitting and testing multinomial processing tree (MPT) models, a class of statistical models for categorical data that involve latent parameters. These parameters are often interpreted as psychological processi

[R] [R-pkgs] TraMineR version 1.6-1

Dear R users, The TraMineR package for mining, describing and visualizing sequences of states or events, and more generally discrete sequential data has been updated to version 1.6-1. Its primary aim is the analysis of biographical longitudinal data in the social sciences, such as data descri

[R] How to maintain class signature in splom

All, I was having trouble trying to create a new class of data and pass it on to splom (in the lattice library). I mentioned this to Martin Morgan after a talk he gave. Following is not so much a question, but rather an answer from Morgan that might be useful to others. Here is the edited

Re: [R] multivariate distributions

On Aug 27, 2010, at 10:42 AM, Rofizah Mohammad wrote: Hi, How can I generate data from multivariate gamma distribution & multivariate beta distribution? I only found command for multivariate normal only. You need to improve your searching skills: rhelpSearch function(string, res

Re: [R] LFA package

Hi: To find functions in R for particular tasks, package sos is an invaluable resource: library(sos) findFn('latent factor analysis') produces 72 hits on my system, not all of which are relevant, but at least it gets you in the vicinity of what you're looking for. In addition to sem, it appears

[R] AIC using nls function

Using the nls function I fit the following model (and some others) to my data. mod1=nls(CLr ~ A-(A-CLi)*exp(-k*d), start = list(A=60,k=0.005)) I would like to rank a set of models using AIC. I calculated AIC as AIC(mod1) However, it appears to use an incorrect number of parameters (3 instead of 2

[R] multivariate distributions

Hi, How can I generate data from multivariate gamma distribution & multivariate beta distribution? I only found command for multivariate normal only. Many thanks in advance :) Regards Rofizah [[alternative HTML version deleted]] __ R-help@r-pr

Re: [R] Band-wise Sum

On Aug 27, 2010, at 9:49 AM, Vincy Pyne wrote: Hi I have a large credit portfolio (exceeding 5 borrowers). For particular process I need to add up the exposures based on the bands. I am giving a small test data below. I would think that cut() would be the accepted method for defining

Re: [R] R-help

However, The "clue" package has the solve_LSAP() function (as pointed out by Hans Werner) that solves the linear sum assignment problem, but it uses the Hungarian algorithm that was requested by you. The lp.assign() function in "lpSolve" package (as pointed out by Girish) also solves the linear s

Re: [R] Nestad ANOVA with random Factors

Hi: I'm kind of wondering why Area and Plot are nested within Day. Are you measuring different areas and different plots on different days? Dennis On Fri, Aug 27, 2010 at 12:43 AM, Pablo R wrote: > > Hi, > > I need a help. I am new in R and I need to run a nested anova with fixed > and > ran

Re: [R] R.matlab package help

I updated R and still have the problem. When I set remote = true, it reports: ??? Undefined variable "java" or class "java.io.InputStreamByteWrapper". Error in ==> MatlabServer at 262 reader = java.io.InputStreamByteWrapper(4096); and I changed remote = false, same problem pops up: ??? Erro

Re: [R] plot data by class tag

Hi, This should do it. I let plot() choose whatever colors it wanted, but if you really wanted you could set them manually. Many thanks to David for the data. plot(temp$V1, temp$V2, pch = as.character(temp$Class), col = temp$Class) # As an alternative using ggplot2 # I do not use "A", "B", "C

[R] Band-wise Sum

Hi I have a large credit portfolio (exceeding 5 borrowers). For particular process I need to add up the exposures based on the bands. I am giving a small test data below. rating <- c("A", "AAA", "A", "BBB","AA","A","BB", "BBB", "AA", "AA", "AA", "A", "A", "AA","BB","BBB","AA", "A", "AAA","

Re: [R] plot data by class tag

On Aug 27, 2010, at 9:09 AM, arvin mer wrote: Hello to all,I have a data file as This came to the list mangled because of your failure to follow Posting Guide advice to use plain text. POST IN PLAIN TEXT. > dput(temp) structure(list(Class = structure(c(1L, 1L, 2L, 3L, 3L, 2L), .Label =

Re: [R] export 4D data as povray density files

Thanks – figure 3 is actually pretty close to what I had in mind. I had forgotten about this package. I still hope to work out the povray route as the end result look really good (and that's probably the main point of such graphics). Best, baptiste On Aug 27, 2010, at 12:37 PM, Deepayan Sark

Re: [R] sqldf syntax

I had checked those references before posting, actually. SQLite has a very limited implementation of the standard. To do a single table update I would not go to sql. It's easy enough to do in R. The problem is when I need to do an update from a left outer join, which I had to do with sqlSave (t

Re: [R] Sorting groups in bwplot chart

On Aug 27, 2010, at 9:03 AM, Jan Hornych wrote: Hi all, I am just curious how to sort the groups in a categorical box plot chart bwplot here is the example: d<-data.frame(sample(rep(month.abb,20), 100), runif(100,1,10)); colnames(d) <- c("Month", "Value"); bwplot(d$Month ~ d$Value); as y

Re: [R] calculate the elasticities by linear.hypothesi commander

Dear Jinghua Xie, This is, I assume, the linear.hypothesis function in the car package, which is deprecated in favour of linearHypothesis. Please see below: > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of Jinghua Xie > Sent:

Re: [R] Sorting groups in bwplot chart

On Aug 27, 2010, at 3:03 PM, Jan Hornych wrote: > Hi all, > > I am just curious how to sort the groups in a categorical box plot chart > bwplot > > here is the example: > > d<-data.frame(sample(rep(month.abb,20), 100), runif(100,1,10)); > colnames(d) <- c("Month", "Value"); > bwplot(d$Month ~

[R] plot data by class tag

Hello to all,I have a data file as Class V1V2A -2.00.0A 0.90.7B 0.10.6C 4.1 0.4C 1.01.9B 1.10.5 I am plotting this data in R as V1 verses V2> temp<-read.table('temp.dat', header=T)> attach(temp)> plot (V1,V2, col='red')> text(x=V1, y=V2, lab

[R] Odp: Sorting groups in bwplot chart

Hi make Month factor, sort its levels and do bwplot as you did ?factor ?levels Regards Petr r-help-boun...@r-project.org napsal dne 27.08.2010 15:03:06: > Hi all, > > I am just curious how to sort the groups in a categorical box plot chart > bwplot > > here is the example: > > d<-data.fram

Re: [R] LFA package

The sem package may be what you are looking for, hth, Ingmar On Thu, Aug 26, 2010 at 9:24 PM, David Joubert wrote: > > Hello all- > > Just wondering if anyone has heard of a package performing latent factor > analysis, similar to what is done with LatentGold. I know that PoLCA does > latent clas

[R] Sorting groups in bwplot chart

Hi all, I am just curious how to sort the groups in a categorical box plot chart bwplot here is the example: d<-data.frame(sample(rep(month.abb,20), 100), runif(100,1,10)); colnames(d) <- c("Month", "Value"); bwplot(d$Month ~ d$Value); as you can see Months are not sorted alphabetically. Does a

Re: [R] R.matlab package help

David, Thanks for your reply, I compile it under the old version of R, I'll update R and run it again. Michael On Fri, Aug 27, 2010 at 7:59 AM, michael wrote: > Hi,all >                      I have a problem running R.matlab package > (under 2.10.1 version). I can set up the matlab ser

[R] calculate the elasticities by linear.hypothesi commander

Dear all If I run the model and get the estimated parameter a11. Then I want to use the estimated parameter to calculate the elasticities by using the formula e11=a11/mw1-1. What I have done is using the command of linear. Hypothesis. > formulas1=dWfresh~dlnPfresh+dlnPfrozen+dlnPsmoke+d

[R] R.matlab package help

Hi,all I have a problem running R.matlab package (under 2.10.1 version). I can set up the matlab server under local machine(run the MatlabServer.m), " And I can use setVariable and evaluate matlab functions in R. But when I ask Matlab to send the value back to R using getVar

Re: [R] Problems when Apply a script to a list

Where exactly did you put the sink() statement? I tried it with a 1000 dataframes and I have no problem whatsoever. Cheers Joris On Fri, Aug 27, 2010 at 6:56 AM, wrote: > Joris, > thank you very much for your help. > It is very helpful for me. > I still have a problem with sink stack although I

Re: [R] Again: wget parameters

On Fri, 2010-08-27 at 13:11 +0200, [Ricardo Rodriguez] Your XEN ICT Team wrote: > OK! > > Gavin Simpson wrote: > > > > > Sorry, I meant: > > > > install.packages("RCurl", type = "source") > > > > G > > > > Here the output. Now, it seems to install without a glitch. I'm afraid I > won't le

Re: [R] Again: wget parameters

OK! Gavin Simpson wrote: Sorry, I meant: install.packages("RCurl", type = "source") G Here the output. Now, it seems to install without a glitch. I'm afraid I won't learn this time the origin of that message! http://xen.net/txt/installSourceRCurl20100827.txt Thanks for your help, Ri

Re: [R] export 4D data as povray density files

On 08/27/2010 08:11 PM, baptiste auguie wrote: Dear list, I wish to visualise some 4D data as a kind of colour / translucent cloud in 3D. I haven't seen such plots in R (but perhaps I missed a feature of rgl). The easiest option I found would be to export the data in povray's df3 (density file)

  1   2   >