similar to: lm for matrix of response...

Dera R users, I have written a function which computes variance, sd, r^2, R^2adj etc. But i am not able to return all of them in return statement. So how to return more than one variable from function. In C i used to return by pointers etc. is there any way like that. Thanks in advance.

I am interested in Neural network models in R. Is there any reference material/tutorial which i can use. Regards, ____________________________________________________________________________________ TV dinner still cooling? Check out "Tonight's Picks" on Yahoo! TV.

dear R users, I have fit the lm() on a mtrix of responses. i.e M1 = lm(cbind(R1,R2)~ X+Y+0). When i use summary(M1), it shows details for R1 and R2 separately. Now i want to use stepAIC on these models. But when i use stepAIC(M1) an error message comes saying that dropterm.mlm is not implemented. What is the way out to use stepAIC in such cases. regards,

I want to plot variation of more than one variable in single plot with different point types and points sizes . Can someone help me to do that, Thanks in advance. vinod ____________________________________________________________________________________ Now that's room service! Choose from over 150,000 hotels

Dear All, Can I use lm() to fit more than one response in single expression. e.g data is a matrix of these variables R1 R2 R3 X1 X2 X3 1 2 1 1 2 3 Now i wnat to fit R1~X1 R2~X2 R3~X3 in turn, and I don't want to do it use loops of couse it it easy to make it using loops,but the proceed is very slow since the data is very big Can't anybody give me some tips or help, my eamil:

Hi, I am trying to use optim to solve a heavy calibration problem. I supply the parameters in vector form. But before entering my target The call is simply: optim(par = parameters, fn = SumLSQ, method = "Nelder-Mead") the function SumLSQ is simply: SumLSQ<-function(parameters, data = timeseries){ print("sumLSQ") nbseries =

Hi all, I want to set up a firewall on CentOS 4.4. I wnat to know the diiffrence between port filtering and packet filtering ? Can iptables do both? Is there another pkg better than this? if so, pls let me know. The purpose of this is to setup a firewall for production use. -- Thank you Indunil Jayasooriya -------------- next part -------------- An HTML attachment was scrubbed... URL:

when i wnat to working with realtime and mysql for any context i have to insert (switch => Realtiem/context at extensions) statment into extensions.conf for example if i want to have 10 context, i have to insert these lines into extension.conf : [context1] switch => Realtiem/context1 at extensions [context2] switch => Realtiem/context2 at extensions [context3] switch =>

I am using Dovecot IMAP and added an IMAP folder to an Outlook client which is also popping another account. Both account folders exist in the same client. I am wondering if I can drag the pop-ed mail to the IMAP folder, thereby transferring old mail to the IMAP folder? Outlook allows me to click and drag folders into the IMAP folder. You see where I am going with this, I wnat to

Dear listers, I have two questions: (1) Is there a way in R to change the base-n of the calculations. I wnat to run some calculations either in binary (base-2) or base-4. Is there a way to specify that in R - to chnage from the decimal? (2) I also want to extract the digits from a larger number and store them as a vector. I could do it through converting top string, parsing the string and

okay. on /posts/show/13 i have show.rhtml, which does render_component(:controller=>''comments'', :action=>''display_comments'', :id=>post.id) on comments/display_comments i iterate through the comments and then, display_comments.rhtml render_component(:controller=>''comments'', :action=>''new'', :id=>@id)

Hi, I've made mysqlf a RPM and installed it. After that i did what i shouldn' do: remove the files it installed. Now, i wnat to install it agai and i can't. Already rebuilt the rpm database ( rpm -vv --rebuilddb), used --force, --replacepkgs, --replacefiles, --justdb, --erase, --initdb, but no good. Been at rpm.org, but found nothing that i already know. Any help would be

I have a router with about 300 clients connecting to it. It has htb with a class per client. I wnat to create a script to prioritise www trafic and ssh trafic over p2p trafic.... this is a sample of what i have now: /sbin/tc class add dev eth0 parent 1:5 classid 1:8012 htb rate 35Kbit ceil 281Kbit prio 6 /sbin/tc qdisc add dev eth0 parent 1:8012 handle 8012 sfq perturb 8 /sbin/tc filter add

Hi all, I am using polr(). The resulting coefficients of first levels are always 0. What to do if I wnat to get the coefficients of the last level 0. For example, suppose x has 3 levels, 1, 2, 3 probit <- plor(y ~ x, data1, method='probit') will get coefficients of level 2, 3 of x, but I want coefficients of level 1, 2 Thank you, Tian [[alternative HTML version deleted]]

Can this be done: Directory structure under ~/News News/agent/nntp/(various server names)/=> thousands of actual news messages/ But also a few files I wnat backed up. Example News/agent/nntp/some.news.server/agent.lib I want agent.lib directory but not all the other stuff under some.news.server ls News/agent/nntp/some.news.server/ agent.lib comp gnu alt The other directories

Hi all, I'm trying to perform an MDS of some data that I have. When I use cmdscale everything is fine and I get some interesting results however, the tends to be low. What I wnat to do is compare this with the Non-Metric MDS using isoMDS or sammon. However, when I try using these I get the following message. Error in isoMDS(x.dist) : zero or negative distance between objects 2 and 4

Dear R experts, it would be very kind if you could help me with two wireplot problems. First, when I make a wireplot and transform it into an .eps using the postscript function the eps-file leaves always a lot of space below the plot, as if it would leave space for a legend or something like that. How can i get the plot into the bottom corner without the space below? The space is not there when

Hey, All In principal component analysis (PCA), we want to know how many percentage the first principal component explain the total variances among the data. Assume the data matrix X is zero-meaned, and I used the following procedures: C = covriance(X) %% calculate the covariance matrix; [EVector,EValues]=eig(C) %% L = diag(EValues) %%L is a column vector with eigenvalues as the elements percent

This is a wrapper around gnm by Turner and Firth to make fitting log-multiplicative models as convenient as possible: it provides simple functions, good starting values, jackknife or bootstrap standard errors, and direct plotting of the results. In addition, it makes it possible to identify scores from RC(M) association models, which gnm does not allow without computing the SVD yourself. It will

This is a wrapper around gnm by Turner and Firth to make fitting log-multiplicative models as convenient as possible: it provides simple functions, good starting values, jackknife or bootstrap standard errors, and direct plotting of the results. In addition, it makes it possible to identify scores from RC(M) association models, which gnm does not allow without computing the SVD yourself. It will