Re: [R] [ how can sample from f(x)~x^(a-1)

khazaei said: > how can sample from f(x)~x^(a-1)*ind(0,min(b,-log(u)) in R? > where a and b is positive constand and 0http://n4.nabble.com/how-can-sample-from-f-x-x-a-1-tp978822p979342.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Is SEM package of R suitable for sem analysis

Besides sem package, OpenMx (http://openmx.psyc.virginia.edu/) is another option for R users. Best 2009/12/27 Joe King : > I am going to take SEM this next quarter in my doctoral program. My > suggestion is to use the program your professor suggests and try to > re-create your models in R using t

Re: [R] How to manipulate tables

On Dec 26, 2009, at 4:56 PM, James Rome wrote: Thanks David. I wanted to calculate the Poisson distribution from my histograms to see how closely they match it. So the formula is something like pprob=((lambda**cnts)/factorial(cnts))*exp(lambda) (Please don't reply privately.) You have tw

Re: [R] Is SEM package of R suitable for sem analysis

Your welcome, I wanted to say I agree with Bruno that the comparative software is very expensive, even more prohibitively so for students and even though MPlus is good I think R will catch up rapidly and even overtake those as people who use these modeling techniques become more integrated into the

Re: [R] Multiple CHOLMOD errors when attempting poisson glmm

On Sat, Dec 26, 2009 at 1:21 PM, Marielle Postava-Davig wrote: > I was actually curious if there was another version of lme4 I could use. You could try glmmPQL in the MASS package. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

Re: [R] Is SEM package of R suitable for sem analysis

Dear Bruno and Joe, Thanks for advising! Reeyarn On Sat, Dec 26, 2009 at 9:32 PM, Bruno Falissard wrote: > A few years ago it could have been true, but now the package has improved > (especially with the bootstrap procedure). > At the moment there is no argument to recommend AMOS. On Sun, De

Re: [R] Is there a problem with lattice axes?

OK, I take it all back" (See original below). I had overlooked A: That the Y Label for the histograms I drew was "Percent of Total", and that there is a parameter 'type for which type="percent" is the default. B: That the histograms I used as an example had N=100 for the case with larges

Re: [R] Is there a problem with lattice axes?

Thanks, Dennis, that does it! And, now that you have pointed it out, I finally found that particular case lurking deep in "?xyplot". This still leaves the question: When the (default) "same" is used, why does histogram() plot the bar-heights as if "free" had been used? The result is that all but o

Re: [R] Question regarding if statement in while loop

Sounds like what you want is either try() or trycatch() -- check out the documentation. Basically, putting your eblest() function inside a try() call allows you to return an error message if eblest crashes without crashing your main loop. BTW, for debugging purposes, I'd recommend creat

Re: [R] [BioC] How to do RMA without summary to probeset level?

You'll need to use calls to bg.correct and the other functions in sequence and save the probe values along the way. See the example in ?normalize.AffyBatch.normalize2Reference in the caret package. Also, don't email R help about this. Keep it on the bioC mailing list. Max On Sat, Dec 26, 2009 a

Re: [R] [GGPLOT] Legends at different layers

Hi Hadley, recently I run this code however got following error (this error seems not be there at 1st time) : > dat <- data.frame(x = rnorm(100)) > dat1 <- data.frame( + x = c(0,0), + y = c(1,0), + Label = c("Point1", "Point2") + ) > > > ggplot(dat, aes(x)) + + geom_histogra

Re: [R] something similar to %include() in sas?

thanks, all, i think source() is the right one i am looking for. On Sat, Dec 26, 2009 at 3:13 PM, David Winsemius wrote: > > On Dec 26, 2009, at 3:10 PM, Wensui Liu wrote: > > i am just wondering if there is an effective way to include other external >> codes into the program. >> >> > ?source >

Re: [R] something similar to %include() in sas?

source("abc.R") if you want to run the code in the abc.R file at that point. Or if you are looking for a macro facility see defmacro in gtools which is based on Thomas Lumley's R News 1/3 article where defmacro is presented. On Sat, Dec 26, 2009 at 3:10 PM, Wensui Liu wrote: > i am just wonderin

Re: [R] something similar to %include() in sas?

?source What is the problem you are trying to solve? Sent from my iPhone. On Dec 26, 2009, at 15:10, Wensui Liu wrote: i am just wondering if there is an effective way to include other external codes into the program. thanks. [[alternative HTML version deleted]] _

Re: [R] something similar to %include() in sas?

On Dec 26, 2009, at 3:10 PM, Wensui Liu wrote: i am just wondering if there is an effective way to include other external codes into the program. ?source -- David Winsemius, MD Heritage Laboratories West Hartford, CT __ R-help@r-project.org

[R] something similar to %include() in sas?

i am just wondering if there is an effective way to include other external codes into the program. thanks. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] Is there a problem with lattice axes?

After answering a previous post by James Rome: "Re: [R] Why do histogram bars vary their width?" I noticed that the lateral axis scales in the lattice histogram which I used for illustration are inconsistent with the frequencies. The example was: set.seed(54321) X <- c(rnorm(100),rnorm(50,1.0,0.

Re: [R] How to manipulate tables

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of James Rome > Sent: Saturday, December 26, 2009 11:03 AM > To: r-help@r-project.org > Subject: [R] How to manipulate tables > > I am sorry to be bothering the list so much. > > I

Re: [R] [BioC] How to do RMA without summary to probeset level?

I think that you misunderstood me. As far as I know, RMA does three things: background correction, quantile normalization, and summary from probes to probesets. I want the probe values after background correction and quantile normalization but before the summary. On Sat, Dec 26, 2009 at 12:07 PM,

Re: [R] How to manipulate tables

On Dec 26, 2009, at 2:02 PM, James Rome wrote: I am sorry to be bothering the list so much. I made a table of counts of flight arrivals by hour: No, you made a list of tables, which is different. cnts=tapply(Arrival4,list(Hour),table). There are up to 15 arrivals in a bin. Why not work

[R] How to manipulate tables

I am sorry to be bothering the list so much. I made a table of counts of flight arrivals by hour: cnts=tapply(Arrival4,list(Hour),table). There are up to 15 arrivals in a bin. > cnts $`0` 1 2 3 4 5 6 7 8 9 10 13 1 2 5 9 2 7 5 4 2 4 1 $`1` 1 2 3 4 3 2 2 1 $`2` 1 3 2 2 . .

Re: [R] if else does not return right value

Thanks David and Patrick, I need to use two statements for the if condition specified. Since the first statement feeds into the second, the ifelse function may just work for the dataset. I'm not sure if there will be issues later but I've tested it and it works. x <- as.matrix(read.table("te

Re: [R] if else does not return right value

See below: On Dec 26, 2009, at 12:58 PM, Muhammad Rahiz wrote: Hi all, I'm not getting the right results for values that are >99 using the if else function. The following illustrates the problem > x <- as.matrix(read.table("test.txt")) > x V1 V2 V3 [1,] 47 1 43 [2,] 83 2 42 [3,]

Re: [R] Question regarding if statement in while loop

Hi Stephanie, it sounds like R's exception handling may help, something like this: foo <- try(eblest(i, dir5, sterr5, weight5, aux5)) if ( class(foo) == "try-error" ) next Take a look at ?try. HTH, Stephan Stephanie Coffey schrieb: Hi all, I'm running R version 2.9.2 on a PC. I'm having a

Re: [R] Why do histogram bars vary their width?

On 26-Dec-09 15:33:45, James Rome wrote: > histogram(~(Arrival4) | as.factor(Hour), type="count", > breaks=16,ylab="Arrival Count", > xlab="Arrival Rate/4",main="Friday EWR A22R D22L Configiration", > layout=c(6,4), par.strip.text=list(cex=0.7)) > > Why do I get plots with different bar widths

[R] Question regarding if statement in while loop

Hi all, I'm running R version 2.9.2 on a PC. I'm having a problem with a loop, and have tried using an if statement within to fix it, but to no avail. Any advice would be appreciated. Here is my code: * eblest <- function(i,dir, sterr, weight, a

Re: [R] Multiple CHOLMOD errors when attempting poisson glmm

Thank you so much for replying. The data set I used for this particular output has 600 entries. The original had 1657 lines so I've already reduced it considerably. Here is a copy of the entire output with lme4 package information: > substrate=read.csv(file.choose(),header=T) > attach(substrate

[R] if else does not return right value

Hi all, I'm not getting the right results for values that are >99 using the if else function. The following illustrates the problem > x <- as.matrix(read.table("test.txt")) > x V1 V2 V3 [1,] 47 1 43 [2,] 83 2 42 [3,] 1 3 41 [4,] 39 4 40 [5,] 23 5 39 [6,] 23 6 38 [7,] 39 7

Re: [R] Is SEM package of R suitable for sem analysis

I am going to take SEM this next quarter in my doctoral program. My suggestion is to use the program your professor suggests and try to re-create your models in R using the SEM package. We are not going to use AMOS though we use EQS, and another prof on campus who teaches it in a different departme

Re: [R] Fwd: Why do histogram bars vary their width?

On Dec 26, 2009, at 10:38 AM, James Rome wrote: I did histogram(~(Arrival4) | as.factor(Hour), type="count", breaks=16,ylab="Arrival Count", xlab="Arrival Rate/4",main="Friday EWR A22R D22L Configiration", layout=c(6,4), par.strip.text=list(cex=0.7)) Why do I get plots with different bar wi

Re: [R] Reading Input file

Use newrate[[i]] <- read.csv(rates[i]) better yet, use a list and lapply newrate <- lapply(paste('rate', 1:4, '.csv', sep=''), read.csv) On Sat, Dec 26, 2009 at 8:09 AM, Maithili Shiva wrote: > Dear R helpers / Johannes Sir, > > Unfortunately it didn't work. > > I am getting following warning

Re: [R] Help with SVM package Kernlab

Hi David, Thank you so much for the pointer. I get it now. I did try the str(testSeq_df) and since it gave me more than 2 factors for each column, I believed that it was fine. I get the point clearly now. Thanks again for all your help. I really appreciate it. Sincerely, vishal On Sat, Dec 26,

[R] Fwd: Why do histogram bars vary their width?

I did histogram(~(Arrival4) | as.factor(Hour), type="count", breaks=16,ylab="Arrival Count", xlab="Arrival Rate/4",main="Friday EWR A22R D22L Configiration", layout=c(6,4), par.strip.text=list(cex=0.7)) Why do I get plots with different bar widths? See http://dl.dropbox.com/u/537118/Friday.p

Re: [R] [ how can sample from f(x)~x^(a-1)

On Fri, Dec 25, 2009 at 8:07 AM, wrote: > Hello all > how can sample from f(x)~x^(a-1)*ind(0,min(b,-log(u)) in R? > where a and b is positive constand and 0 If the idea is that X is a random variable, then you need to decide what kind of random variable it is. For example, if you wanted to as

Re: [R] Help with SVM package Kernlab

On Dec 26, 2009, at 3:53 AM, Vishal Thapar wrote: Hi All, Thank you for your replies so far. I was hoping I could get some more input from you on this issue. It seems to me that I have hit a dead end here and would really appreciate some feedback. I have followed all the suggestions you

Re: [R] Help with SVM package Kernlab

Perhaps my response is not showing up on the r-help list due to the large number of recipients generated by the reply-all option triggering some sort of spam filter. So I am trimming them. On Dec 26, 2009, at 3:53 AM, Vishal Thapar wrote: Hi All, Thank you for your replies so far. I was ho

Re: [R] Reading Input file

First of all there seems to be something wrong with your equation > rate_name[i] = (paste(`rate', i, ‘.csv`, sep = ‘’)) Try this rate_name[i] <- paste("rate",i,".csv", sep="") I am not --- On Sat, 12/26/09, Maithili Shiva wrote: > From: Maithili Shiva > Subject: Re: [R] Reading Input file >

Re: [R] Is SEM package of R suitable for sem analysis

A few years ago it could have been true, but now the package has improved (especially with the bootstrap procedure). At the moment there is no argument to recommend AMOS. Of course, some software like Mplus are actually better, but very specialized (and rather expensive). This is only a personal po

Re: [R] Reading Input file

Dear R helpers / Johannes Sir,   Unfortunately it didn't work. I am getting following warning messages.   1: In neweate[i] = read.csv(rates[i]) :   number of items to replace is not a multiple of replacement length. 2:  similar message 3:  similar message 4: similar message   Let me be speci

Re: [R] Help with SVM package Kernlab

Hi All, Thank you for your replies so far. I was hoping I could get some more input from you on this issue. It seems to me that I have hit a dead end here and would really appreciate some feedback. I have followed all the suggestions you have mentioned but they still this is stuck. Earlier I thoug

[R] Is SEM package of R suitable for sem analysis

Dears, I'm a college student and In doing my statistics homework. I use R with SEM package as my tool for sem analysis, but my teacher told me AMOS is more suitable for such analysis. Could someone help tell me whether it is true that some commercial software is better accepted in academic field

Re: [R] String question

Knut Krueger wrote: Will this do? temp <- paste("m", 1:3, sep="",collapse=",") Unfortunately not, because I explained the example not detailed enough. The string could have different Items, like November, December, Monday, Tuesday, Daylight and so on Therefore I must count the Items of t

Re: [R] Reading Input file

On Sat, Dec 26, 2009 at 10:21 AM, Maithili Shiva wrote: > > > Dear R helpers > > I have some files in my say 'WORK' directory and the file names are say > rate1.csv, rate2.csv, rate3.csv, rate4.csv > > Because of some other requirement, I need to run the following commands > > n = 4 > > rates = NU

[R] Reading Input file

Dear R helpers   I have some files in my say 'WORK' directory and the file names are say rate1.csv, rate2.csv, rate3.csv, rate4.csv   Because of some other requirement, I need to run the following commands   n = 4    rates = NULL   for (i in 1:n) rates[i] = (paste(`rate', i, ‘.csv`, se