result <- multinom(Y ~ X+age, data = tdata)
z <- summary(result) $ coefficients/summary(result) $ standard.errors
p <- (1-pnorm(abs(z),0,1))*2
out <- data.frame(cbind(names(tdata)[q], t(summary(result) $ coefficients[,2]), t(p[,2])))
timo <- rbind(timo,out)
z <- summary(result)¥coefficients/summary(result)¥standard.errors
p <- (1-pnorm(abs(z),0,1))*2
out <- data.frame(cbind(names(tdata)[q], t(summary(result)¥coefficients[,2]), t(p[,2])))
timo <- rbind(timo,out)