search for: daalgard

...tapply I know I cannot use length to get a count where there are NAs. summaryBy appears to work the same way. I do know how to get a count of the entire column using sum: sum(!is.na(lh.forgotten$RLWTEST)) However, this does not give me a count split up by my factor (STATUS). I have looked through Daalgard (2002) and Verzani (2005), and have searched the help files, but with no luck. Thank you in advance for your help. I love R and am interested in making it more accessible to social scientist types like me. I know it can do everything SPSS can and more, but sometimes the very simplest things seem t...

.... News (with respect to rw0613) -------------------------------- (a) I rewrote the console using a portable toolkit (GraphApp, see http://www.cs.usyd.edu.au/~loki/graphapp). Now, output is buffered until some input is required. This means that it is faster (for a good test, I learned it from Peter Daalgard, try `ls("package:base")' ). Of course, buffering can be disabled (you need to disable it ONLY if your command require a lot of time and you want to look to the output line by line). In addition, there are some command history support. For a mini help on the new console, use the help/...

.... News (with respect to rw0613) -------------------------------- (a) I rewrote the console using a portable toolkit (GraphApp, see http://www.cs.usyd.edu.au/~loki/graphapp). Now, output is buffered until some input is required. This means that it is faster (for a good test, I learned it from Peter Daalgard, try `ls("package:base")' ). Of course, buffering can be disabled (you need to disable it ONLY if your command require a lot of time and you want to look to the output line by line). In addition, there are some command history support. For a mini help on the new console, use the help/...

.... News (with respect to rw0613) -------------------------------- (a) I rewrote the console using a portable toolkit (GraphApp, see http://www.cs.usyd.edu.au/~loki/graphapp). Now, output is buffered until some input is required. This means that it is faster (for a good test, I learned it from Peter Daalgard, try `ls("package:base")' ). Of course, buffering can be disabled (you need to disable it ONLY if your command require a lot of time and you want to look to the output line by line). In addition, there are some command history support. For a mini help on the new console, use the help/...

Dears, I saw in the R-project site some R-books. However, I'm new in this community and I didn't figure out what are the best books. Can you suggest me some "reference" books? My intentions with R is concerned to Artificial Intelligence simulations, Classification and general Statistics (e.g., regression, multivariate regression, etc). Thanks for any help, Anderson

...h to get a > count where there are NAs. summaryBy appears to work the same > way. I do know how to get a count of the entire column using sum: > sum(!is.na(lh.forgotten$RLWTEST)) > > However, this does not give me a count split up by my factor > (STATUS). I have looked through Daalgard (2002) and Verzani > (2005), and have searched the help files, but with no luck. > > Thank you in advance for your help. I love R and am > interested in making it more accessible to social scientist > types like me. I know it can do everything SPSS can and more, > but sometim...

Hello I have to solve a multiple linear regression. Most programs like Excel or Mathlab only support 5-10 dimensions. Now I have installed CRAN and I have no clue what to do next. At the moment I am entering my data into an excelsheet (for quick copy- paste). The Y-array will be 20 columns (=dimensions) and 128 rows (=variables). The X-array may also be 128 rows in Excel. Which package do I need

Greetings, I would like to use a cubic spline or smoother to model the fixed effects within nlme. So far the only smoother I have been able to get to run successfully in nlme is smooth(). I tried smooth.spline: fixed=list(lKa~1,lCL~smooth.spline(BSA, df=3)) the error I got was the following. Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid

Dear R'ers, when I e.g. install the integrate2.2-2 package under R1.1.1 (debian2.2) atlas is linked into the library: Installing package `integrate' ... libs gcc -shared -lc -shared -o /usr/lib/R/library/integrate/libs/integrate.so adapt.o adapt_callback.o bsrl.o flfm.o nxprt.o symrl.o wht.o -L/usr/local/lib -lf77blas -lcblas -latlas -lg2c -lm -L/usr/lib/gcc-lib/i386-linux/2.95.2

Hi, Can someone help me out to create a (for?) loop for the following procedure: x=rnorm(250,0,0.02) library(timeSeries) x=timeSeries(x) P=1000 loss=-P*x loss

Dear all, I have something about function outer() that I can't understand. Just see the following example. The two NaNs are due to 0/0, but I can't figure out the cause of the last two errors. I wonder if some one can explain this for me. ___________________________________________________________________ > sx=rbinom(10,1,0.5);ot=rbinom(10,1,0.5);ag <- rbinom(10,100,0.3);ho <-

Hi all, I have been struggling with ANOVAs on R. I am new to R, so I created a simple data frame, and I do some analyses on R just to learn R and then check them on SPSS to make sure that I am doing fine. Here is the problem that I've run into: when we use the aov function, it uses SS Type I as default (on SPSS it is Type III). Then I used the Anova function under cars package using the

Hi Again, First of all thank you for all the responses to my previous query. Your answers were very helpful and I did the job ;-). Now I hope you can answer as quick the following (sorry I am invading you with trivial questions): Let?s use again the following data.frame example: DF <- data.frame(x=rnorm(5), y=rnorm(5)) I want to obtain a new data.frame (or matrix) that contains only n

Hi, I noted a discrepancy between R and openepi when I ran a fisher test with the same matrix. In R: > a=matrix(c(1,2,6,17), nrow=2) > a [,1] [,2] [1,] 1 6 [2,] 2 17 > fisher.test(a, conf.int=T) Fisher's Exact Test for Count Data data: a p-value = 1 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 0.02061498

...> Subject: RE: [R] cubic spline/smoother with nlme To: "Suzette Blanchard" <suzette at sdac.harvard.edu> Cc: r-help at stat.math.ethz.ch Message-ID: <CFBD404F5E0C9547B4E10B7BDC3DFA2F0B16072D at usrymx18.merck.com> Content-Type: text/plain Suzette, In addition to Professor Daalgard's suggestions of ns() and bs(), you could also try out rcs() from Frank Harrell Design package (you may need his Hmisc package as well). This function helps to fit natural (restricted cubic splines), and have been very useful for me in practice to use in tandem with lme() in modeling longitudin...